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1 Introduction 



It is well known that the equations of motion for a mechanical system with nonholonomic constraints do 
not arise from a variational principle in the usual sense. As a consequence, they cannot be formulated 
as a classical Hamiltonian system. Instead, they are written with respect to an almost Poisson bracket 
that fails to satisfy the Jacobi identity. This formulation has its origins in [3^ [251 EI] and others. 

On the other hand, after a symmetry reduction, the equations of motion of a number of exam- 
ples allow a Hamiltonian formulation (sometimes after a time reparametrization) and one talks about 
Hamiltonization^ (see [IIl[IlllI5l[20l[I7l|19l[22l|28j and others). 

In this paper we employ recent developments of Poisson geometry to study this phenomenon from a 
geometric perspective. We use gauge transformations by 2-forms as introduced by Severa and Weinstein 
in [29] to construct different almost Poisson brackets describing the dynamics of the same nonholonomic 
system. Although our interest is in almost Poisson geometry, we consider more general objects known 
as almost Dirac structures |12j . as they provide the most natural setting for the definition and study 
of gauge transformations. 

We illustrate the need for our methods by working out the Hamiltonization of the motion of rigid 
bodies that are subject to generalized rolling constraints. These are nonholonomic constraints that 
relate the angular velocity uj of the body and the linear velocity x of its center of mass in a linear way 
(i.e. X = Au for a 3 x 3 matrix A). This type of constraints contain the celebrated Chaplygin sphere 
problem as a special case. The latter concerns the motion of an inhomogeneous sphere whose center of 
mass coincides with its geometric center, that rolls without slipping on the plane. As a consequence of 

difFerent meaning to Hamiltonization is given in where the authors study the unreduced system in connection 
with the inverse problem of the calculus of variations. 
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the Hamiltonization, we are able to show complete Liouville integr ability of the reduced dynamics of 
any rigid body subject to generalized rolling constraints. 

Incidentally, during our discussion, we discover that twisted Poisson brackets [29J appear in the 
study of nonholonomic systems. In particular, we show that in the original physical time (before the 
time reparametrization), the reduced dynamics of the Chaplygin sphere are formulated in terms of a 
twisted-Poisson bracket. Although these structures do not in general satisfy the Jacobi identity, they 
possess a fair amount of properties, including foliations, that might imply an interesting interplay with 
dynamical features. To date, the interest in these brackets has been mainly geometrical. 

1.1 Hamiltonization 

Perhaps the most interesting example of Hamiltonization concerns the Chaplygin sphere. Even though 
the formulation and integration of the equations of motion by Chaplygin dates back to 1903 [TU], the 
Hamiltonian structure of the reduced equations (after a time reparametrization) was only discovered 
in 2001 by Borisov and Mamaev [6]. 

Recently, Jovanovic [22] proved that the multidimensional version of the Chaplygin sphere problem 
introduced in ^16j is also integrable and Hamiltonizable when the vertical angular momentum is zero. 
This gives a partial solution to a problem that remained open for many years. His approach to prove 
integrability involves in a crucial way the Hamiltonization of the problem. Another important example 
where the integration of a nonholonomic system follows from its Hamiltonization is the multidimensional 
Veselova system treated by Fedorov and Jovanovic in [15j . We also mention the recent work of Ohsawa, 
Fernandez, Bloch and Zenkov [28j in connection with Hamilton-Jacobi theory. 

The relationship between Hamiltonization and integrability may have been the original motivation 
for Chaplygin to consider the problem of Hamiltonization back in 1911 |11] . In this work, Chaplygin 
proved the famous Chaplygin reducing multiplier Theorem that applies to the so-called G- Chaplygin 
systems. These are nonholonomic systems with the property that the tangent space to the orbits of 
a symmetry group G exactly complements the constraint distribution on the tangent space TQ of the 
configuration manifold Q. Stated in modern geometric terms, the Theorem says that if the shape space 
Q/G is two-dimensional, and the reduced equations have an invariant measure, then they can be put 
in Hamiltonian form in the new time r defined by dr = ^dt. The positive function ^ : Q/G ^ M 
is known as the reducing multiplier'^ . There is a very neat interpretation of the multiplier ^ in terms 
of the invariant measure and as a conformal factor for an almost symplectic form that describes the 
dynamics, see [HI |T5l [20l [28] . This interpretation suggests that geometric methods may be useful to 
understand Hamiltonization in more general scenarios of nonholonomic systems with symmetry. 

Recently, Fernandez, Mestdag and Bloch [17j, derived a set of coupled first order partial differential 
equations for the multiplier ^ for G-Chaplygin systems whose shape space has arbitrary dimension. 
Even more, the set of equations found by the authors applies to general nonholonomic systems with 
symmetry that are not necessarily G-Chaplygin. This is done by writing the reduced equations of 
motion in Hamilton- Poincare-D'Alembert form as described in ^ i4j. The issue of Hamiltonization 
is thus reformulated as a problem of existence of a solution for the aforementioned system of partial 

^We denote the reducing multiplier by ^ instead of to be consistent with our exposition which takes the Poisson 
rather than the symplectic perspective. 
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differential equations. 

Our approach to Hamiltonization contains the same degree of generahty but is more intrinsic. 
Denote by Xfi the vector field describing the dynamics on the reduced space TZ and by T-L-ji the reduced 
Hamiltonian. We formulate the reduced equations of motion in almost Poisson form with respect to 
a collection of bivector fields t^^^^b. Each member in this collection describes the reduced dynamics 
in (almost) Hamiltonian form (i.e. {'K^^^B)\d'H'ji) = and arises as the reduction of a bivector 

field TTj^ associated to a bracket {•, -j^. Such bracket {•, ■}^^^ is obtained through what we define as a 
dynamical gauge transformation by a 2- form B of the noholonomic bracket {•, -jnh defined in [301 [251 EI] 
(see discussion in subsection 11.21 below) . 

Having said this, we reformulate the issue of Hamiltonization by requiring that one of the bivector 
fields TT^^^B in the collection described above is conformally Poisson, i.e. 

[^■^rcdB,^1T,,^B]=0, (1.1) 

for a positive function cp,^ and where [•, ■] is the Schouten bracket. The scaling of tt^^^b by if is 
dynamically interpreted as the time reparametrization dr = ^dt (see Section [4. 4p . Note that, for each 
^rcdS ! equation p.l|) locally defines a set of coupled first order partial differential equations for (/?. This 
seems to be in agreement with the results of Fernandez, Mestdag and Bloch |17j . 

If the symmetries are of G-Chaplygin type, the bivector fields tt^^^b are everywhere non-degenerate 
and the equations of motion can be written with respect to the associated almost symplectic form 
^^^^b. If a multiplier if satisfying (jl.ip exists, then i^^^^B is conformally closed, (d(^^^rods) = 0), and 
one speaks of Chaplygin Hamiltonization [14J. 

The term Poissonization was introduced by Fernandez, Mestdag and Bloch in [T7| to refer to the 
case where the bivector field vr^^^s satisfying (jl.ip is degenerate. Their motivation to distinguish this 
case is to study the relationship between Hamiltonization and the existence of invariant measures for 
the reduced equations. To simplify our exposition and to treat the problem in a unified manner, we will 
not use their terminology and simply talk about Hamiltonization whenever there exists a solution to 
(jl.ip . We give a discussion on the existence of invariant measures for nonholonomic systems admitting 
a Hamiltonization in this generality in subsection 14. 4[ 

In general terms, the main contributions of this paper to the problem of Hamiltonization are the 
clear geometric formulation of the problem using recent developments on the field of Poisson geometry 
(mainly those in Severa, Weinstein [2Q] ) and the illustration of the usefulness of these techniques in the 
study of rigid bodies subject to generalized rolling constraints. 

1.2 Gauge transformations in nonholonomic mechanics 

The Hamiltonization of the Chaplygin sphere can be obtained in two ways: 

In the first one, described in Borisov and Mamaev [7], one performs the reduction in two stages. 
In the first step one reduces the translational symmetries of the problem corresponding to the abelian 
action of M^. On the second stage one reduces the internal symmetry corresponding to rotations about 
the vertical axis described by the action of S^. The second stage is performed in Borisov and Mamaev [7] 

■^In fact the function if : TZ ^ is not arbitrary. It is required to be basic with respect to the fibered structure of 
the reduced space TZ. Just as for G-Chaplygin systems, we think of </? : Q/G — >■ R^. 
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using Routh's reduction method. The Hamiltonian interpretation of the second reduction is dehcate 
and is studied in Hochgerner, Garcia-Naranjo [20]. Effectively, it is shown that in order to further 
reduce the system in a Marsden-Weinstein fashion, the geometric data needs to be modified and the 
authors propose a method for this cahed truncation. 

A second manner to achieve the Hamihonization of the Chaplygin sphere is described in Garcia- 
Naranjo |19j . In this approach one first formulates the equations of motion in terms of an affine almost 
Poisson bracket that differs from the usual nonholonomic bracket {•, considered in |30 1 1251 [2T] and 
others. In this framework, the Hamiltonization can be obtained as a single (one step) reduction by the 
symmetry group of rigid transformations on the plane, SE(2). 

In this paper we further elaborate on the latter approach. Our main tool for this is to incorporate 
gauge transformations as introduced in Severa and Weinstein |29j . 

Recall that the nonholonomic bracket {•, considered in [30 1 1251 [2T] is defined on functions on the 
constraint phase space M.. The vector field X-^h on M that describes the nonholonomic dynamics is 
(almost) Hamiltonian with respect to this bracket and with respect to the constrained Hamiltonian T-Lm- 
Now, let be a 2-form on A^. li B satisfies a certain technical condition, then a gauge transformation 
of the nonholonomic bracket by B defines a new bracket {•, -j^ on M with the same (non-integrahle) 
characteristic distribution. As a consequence, the (almost) Hamiltonian vector fields associated to 
{•, -j^ satisfy the nonholonomic constraints. If in addition, the 2-form B satisfies ij^ ^5 = 0, then we 
say that B defines a dynamical gauge transformation and the vector field is (almost) Hamiltonian 
with respect to the gauged bracket {-j-}^ and the constrained Hamiltonian I-Lm- Iii this way, we 
distinguish a family 5^ of almost Poisson structures that describe the dynamics of our nonholonomic 
system corresponding to different dynamical gauge transformations. We show (Remark 14. 4|) that the 
affine almost Poisson brackets defined in Garcia-Naranjo [19] are particular members in 5^. 

The main motivation to consider the large family 5^ of almost Poisson brackets for our nonholonomic 
system is to have a larger choice of structures to describe the reduced dynamics and hope to find one 
amongst them that Hamiltonizes the problem. In particular, the Hamiltonization of the Chaplygin 
sphere arises as the reduction of a dynamically gauged bracket {■, -j^ that differs from the standard 
nonholonomic bracket {•,•}„!,• 

In fact, the members of the family 5^ behave quite differently after reduction. Our examples show 
that while some of them might yield a true Poisson structure after reduction, others yield an almost 
Poisson bracket with a non-integrable characteristic distribution. 

In order to study the gauge transformations of almost Poisson brackets, we formulate the dynamics 
of nonholonomic systems on almost Dirac structures [12j. These are more general geometric objects 
that provide the framework in which gauge transformations are more natural. These structures had 
been already considered in connection to nonholonomic mechanics by Yoshimura and Marsden [321 133j. 
and by Jotz and Ratiu |23j. However, the issue of Hamiltonization and the incorporation of gauge 
transformations are not treated in these works. 

1.3 Twisted Poisson brackets in nonholonomic mechanics 

It is well known that nonholonomic systems are formulated in terms of almost Poisson brackets that 
fail to satisfy the Jacobi identity. However, very little research, if any, has been done in understanding 
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how far away these brackets are from being true Poisson. 

A very strong property of a Poisson manifold is that the characteristic distribution is integrable 
and defines a fohation by even dimensional leaves. This property is also shared by conformally Poisson 
brackets whose bivector field tt satisfies [(/?7r, ipir] = 0, for a certain positive function ip called the 
conformal factor. These brackets have been considered in the study of Hamiltonization of nonholonomic 
systems, and, as mentioned before, the conformal factor ip defines the time reparametrization dr = ^ dt. 

Another example of almost Poisson brackets that possess a foliation by even dimensional leaves, 
is given by twisted Poisson brackets that were introduced by Klimcik and Strobl in [23] and later in 
Severa and Weinstein [29j from a more geometric point of view. Twisted Poisson brackets correspond 
to almost Poisson structures whose associated bivector field tt satisfies 

^[7r,7r] = vr" ((/)), 

for a certain closed 3-form (f). To our knowledge, the present paper is the first one to explore the 
connection between this type of structures and nonholonomic mechanics. 

In this paper we show that an almost Poisson structure with a regular (constant rank), integrable 
characteristic distribution is twisted (Corollary 13. 7p . As a consequence we show (Remark 13. 8p that 
the reduced equations of the classical Veselova problem [3T] can be formulated in terms of a twisted 
Poisson bracket in the original physical time (prior to any time reparametrization). 

We also show (Theorems 15.71 and 15. 8p that the reduced equations of some examples of rigid body 
motion with generalized rolling constraints, that contain the Chaplygin sphere as a special case, are 
described by a twisted Poisson bracket in the original physical time. Moreover, we give an explicit 
formula for the twisting closed 3-form (/>. 

1.4 Outline and main results of the paper 

The paper is organized as follows. In Section [2] we introduce our motivating examples, rigid bodies 
subject to generalized rolling constraints. As mentioned before, these are nonholonomic constraints 
that relate the linear velocity of the body to its angular velocity via a 3 x 3 matrix A. After writing 
down the reduced equations of motion, we define two different (almost) Poisson structures for the 
reduced equations according to the rank of A, that varies from to 3. For each value of the rank 
of A, we show that only one of the brackets is Poisson (conformally Poisson if rank A = 1,2) while 
the other one possesses a non-integrable characteristic distribution. The geometric interpretation and 
construction of these brackets is one of the main goals of the paper and is postponed to Section [5l after 
the necessary tools are developed in Sections [3] and [H 

In Section [3] we develop the geometric background needed for our purposes. We focus on almost 
Dirac structures and their gauge transformations, introduced respectively in [12] and [29], and we 
collect some new results that are important in our study of nonholonomic systems. Proposition 13.11 
gives a characterization of regular almost Dirac structures that is used in Corollarv 13.41 to describe the 
structure of almost Poisson brackets having a regular characteristic distribution. Corollary 13.71 shows 
that an almost Poisson bracket possessing a regular, integrable, characteristic distribution is twisted. 
In fact this result is proved in the more general setting of almost Dirac structures in Theorem 13.51 We 



6 



also mention Theorem 13.111 that asserts that any two regular almost Dirac structures defining the same 
distribution are gauge related. 

In Section [J] we make the connection between the geometric methods developed in Section [3] and 
nonholonomic mechanics. In particular, we construct the nonholonomic bracket of [301 l25t [2T] using 
Corollary 13.41 and the framework for nonholonomic mechanics described in Bates and Sniatycki [2]. In 
Proposition 14.31 we show that the dynamics associated with this bracket coincide with the formulation 
of nonholonomic mechanics on almost Dirac structures considered in \32 \ \33 [ [23]. Next, we define the 
notion of dynamical gauge transformations for a nonholonomic system, and define a family 5 of almost 
Poisson brackets, possessing the same characteristic distribution, and that describe our nonholonomic 
system. Finally, we discuss the reduction of these brackets in the presence of symmetries and introduce 
our working definition of Hamiltonization. 

In Section O we resume the study of rigid bodies subject to generalized rolling constraints. In 
subsection 15 . 1 1 we show that the brackets given in Section [2] to describe the reduced dynamics, arise as 
a reduction of different members of the family ^ (Theorems 15.31 and 15 . 4p . In subsection 15 . 21 we establish 
the Hamiltonization of the reduced equations in detail and we conclude their integrability. Finally, in 
subsection 15.31 we focus on the twisted nature of the brackets that Hamiltonize the problem for the 
cases Rank A = 1,2, prior to the time reparametrization. 



2 Motivating Examples: Rigid bodies with Generalized Rolling Con- 
straints 



Consider the motion of a rigid body in space that evolves under its own inertia and is subject to the 
constraint that enforces the linear velocity of the center of mass, x, to be linearly related to the angular 
velocity of the body u, i.e., 

X = rAu. (2.2) 

Both vectors x and u> belong to and are written with respect to an inertial frame. The constant 
scalar r has dimensions of length and is a natural length scale of the system. The dimensionless constant 
3x3 matrix A is given and satisfies certain conditions that are made precise in the following Definition. 

Definition 1. The matrix A is said to define a generalized rolling constraint if it satisfies one of the 
following conditions according to its rank: 



(i) A= (^^ J ^ , with C G S0(2), if rank^ = 3. 
{ii) A = ( M , with C G S0(2), if rank^ = 2. 



.0 . 

[Hi) A = e^, if rank^ = 1, where eg is the third canonical vector in M^, and T denotes transpose. 

(iv) A = if rank^ = 0. 

The above conditions on A can be relaxed (see Remark 12.11 ahead) . However, for simplicity, we will 
assume that A has the form given by one of the items of the above Definition. If A satisfies any of the 
the conditions of the above Definition we say that (|2.2p is a generalized rolling constraint. 
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Our terminology is motivated by a particular example: the Chaplygin sphere. The problem, intro- 
duced by Chaplygin in 1903 [10], concerns the motion of a ball whose center of mass coincides with its 
geometric center that rolls on the plane without slipping. In this case, the matrix A is given by 



A 



and r is the radius of the sphere. 

The motion of the Chaplygin ball has been the subject of much research to our days. An important 
property is that (after a time reparametrization) the reduced equations can be given a Hamiltonian 
structure [6l[T9]. The geometry of the Hamiltonization of the problem is intricate. In order to study 
this phenomenon in a mathematically systematic fashion, we consider more general possibilities for the 
matrix A. 

The crucial property of A that determines many of the dynamical and geometrical features of the 
problem is its rank. The Chaplygin sphere corresponds to the case rank A = 2. Another familiar case 
occurs when rank A = 0. In this case the constraint (j2.2p becomes x = which can be interpreted as 
a conservation law for the free system that states that the center of mass of the body is at rest in the 
inertial frame. The motion of the system reduces to that of the classical free rigid body. 

We will also consider the cases where the rank of A equals 3 and 1 which, to our knowledge, have 
not yet been considered in the literature. 



2.1 Generalities 

The configuration space for the system is Q = S0(3) x W^. Elements in Q are of the form q = (5',x) € 
SO (3) X W^. The vector x € is the position of the center of mass in space and the orthogonal matrix 
g specifies the orientation of the ball by relating two orthogonal frames, one attached to the body and 
one that is fixed in space. We will assume that the body frame has its origin at the center of mass and 
is aligned with the principal axes of inertia of the body. These frames define the so-called space and 
body coordinates respectively. 

Recall that the Lie algebra so (3) can be identified with equipped with the vector product via 
the hat map: 

■n = {vi,m,V3) ^ v = { m o -vi \ ■ (2-3) 






-m 


m 







-m 


-m 


m 






Given a motion {g{t),x{t)) € Q, the angular velocity vector in space coordinates, cj € M^, and the 
angular velocity vector in body coordinates, E M^, are respectively given by 

u{t) = g{t)9-\t), n{t) = g~\t)g{t), 



and satisfy ft = g ^u. It will be useful to write the constraint (j2.2p in terms of the body angular 
velocity as 

± = rAgn. (2.4) 
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The kinetic energy of the rigid body defines the Lagrangian L : TQ ^ M by 

£(g,5,x,i) = i(m).n + ^||i||2, (2.5) 

where "•" denotes the Euchdean scalar product on M'^, m is the mass of the body and the 3x3 diagonal 
matrix I is the inertia tensor with positive entries Ii,l2,l3- 

Remark 2.1. It is not hard to see that if the space axes are rotated by an element h € S0(3), the 
Lagrangian C is invariant and the constraint (j2.2p is rewritten as 

X = rh^^ Ahu . 

Therefore, the conditions for A given in Definition [1] can be relaxed by allowing conjugation by matrices 
h G S0(3). o 



2.2 The equations of motion 

Let p = mx be the linear momentum of the body. In accordance with the Lagrange-D'Alembert 
principle, the constraint forces must annihilate any velocity pair (x,ri) satisfying (|2.4p . Therefore, the 
equations of motion are given by 

p = /X, in = m X n - rg-'^A^n, (2.6) 

where "x" denotes the vector product in and the multiplier G is determined uniquely from 
the constraint (|2.4p . 

Differentiating ()2.4p and using gCl = we find fj, = mrAgft. Thus, the second equation in ()2.6p 
decouples from the first to give 

m = m x 17 - mr'^g-^A^Ag tl. (2.7) 

In principle, this equation should be complemented with the reconstruction equation g = gCl. It will be 
shown ahead that it suffices to consider the evolution of the Poisson vector 'y := g^^ea that represents 
the vector 63 written in body coordinates. A direct calculation gives 

7 = 7 X fi. 

The decoupling in ()2.6p is due to the presence of symmetries that will be discussed in detail in Section 
12.41 Once this equation is solved for [g, fl), we obtain p = mrAgft that follows from ()2.4p . 

We introduce the kinetic momentum K € by 

K ■.= m + mr^g-'^A'^Ag^l. (2.8) 
This definition of the kinetic momentum allows us to define the (reduced) Hamiltonian 

nn = \{K-n), (2.9) 



9 



which coincides with the kinetic energy on the constraint space A4. 
A direct calculation using ()2.7p gives our final set of equations 

K = Kxn, 7 = 7xn. (2.10) 

To understand why the above equations define a closed system for (K, 7) G x M^, and to understand 
their structure, it is useful to perform a separate study for different values of the rank of the matrix A. 
This will also show that the Hamiltonian Ti-ji can be considered as a function of K and 7. 



2.3 A pair of (almost) Poisson brackets for the equations of motion 

For each value of the rank of A we will give two different brackets that define the equations of motion 
(|2.10|) . with respect to the reduced Hamiltonian T-Ln. In general, these brackets are almost Poisson, 
i.e. they do not satisfy the Jacobi identity but we will argue that one of them is more convenient than 
the other. They will be denoted by {•,-}Rankj and {^-l^ankj where j denotes the rank of the matrix 
A. Both brackets define the equations of motion (|2.10p in the sense that the directional derivative of 
any function / = /(7,K) G C°°(R3 x M^) along the fiow is given by / = {/, ^^T^laank, = {f,nn}Ln^r 
The geometric interpretation of these brackets is the subject of the subsequent sections. Concretely, 
in section [5] (Theorems 15.31 and 15. 4p . we will show that the bracket {•,-}Rankj arises as the reduction 
of the nonholonomic bracket introduced in |30j, and that {•, -jRankj arises as the reduction of a gauge 
transformation of the nonholonomic bracket. 

The following definitions will be useful in our discussion of the utility of the two brackets: 
Definition 2. Let P be a manifold equipped with an almost Poisson bracket {•, •}. 

1. The (almost) Hamiltonian vector field Xf of a function / G C°°{P) is the vector field on P 
defined as the usual derivation Xf{g) = {g, /} for all g G C°°{P). 

2. The characteristic distribution of {•,•} is the distribution on the manifold P whose fibers are 
spanned by the (almost) Hamiltonian vector fields. 

3. Due to Leibniz condition of {•, •}, there is a bivector field vr G T{/\^{TP)) such that for f,g& 
C°°{P) we have iT{df,dg) = {f,g}- We say that vr is the bivector field associated to {•,•} and 
we denote by vr^ : T*P — )• TP the map such that /3(7r^(Q)) = 7r(a,/3). We will occasionally refer 
to bivector fields simply as bivectors. Note that the characteristic distribution is the image of 
TT^ and the Hamiltonian vector field Xf = — 7r^((i/). The 3-vector field ['7r,7r], where [•,■] is the 
Schouten bracket, may be different from zero, and it measures the failure of the Jacobi identity 
through the relation 

{/, {9, h}} + {g, {h, /}} + {h, {/, g}} = ^i[^,,](<i/, dg, dh), (2.11) 
for /,5,/iGC-(P). 

4. The bracket is called conformally Poisson if there exists a strictly positive function ip G C°°{P) 
such that the bracket ip{-, •} satisfies the Jacobi identity, i.e. [(pir, ipn] = 0. 
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The well known symplectic stratification theorem states that the characteristic distribution of a 
Poisson bracket is integrable and its leaves are symplectic manifolds. Since multiplication of an almost 
Poisson bracket by a positive function does not change the characteristic distribution, a necessary 
condition for an almost Poisson bracket to be conformally Poisson is that its characteristic distribution 
be integrable. 

We now come back to the discussion of our example for the different values of the rank of A. 



If A has rank 3. In this case A~^ = AJ- and K = (I + mr'^E)Q, where E denotes the 3x3 identity 
matrix. It follows form ()2.10p that the rotational motion of the body is the same as that of a free rigid 
body whose total inertia tensor is I + mr'^E. It is trivial to write = (I + mr'^E)~^'K. and it is clear 
that equations ()2.10p define a closed system in M? X M?. 

The two brackets for the system for functions f,g £ C°°(M^ x M^) are given by 



(2.12) 



The above brackets are quite different. On the one hand, the bracket {•, -laj^nka satisfies the Jacobi 
identity. It in fact coincides with the Lie- Poisson bracket on the dual Lie algebra se(3)*. On the other 
hand, the bracket {•, - jRanka is not even conformally Poisson as the following Proposition shows. 

Proposition 2.2. The characteristic distribution of the almost Poisson bracket {•, -iRanks defined in 
()2.12p is not integrable. 



Proof. The (almost) Hamiltonian vector field Xf of a function / G C°°(M x M ) corresponding to the 
bracket {■, -jRanks is given by 

2r.x 9f df\ d ( df\ d 

and it is annihilated by the non-closed one-form 

X = 7 • dK + (K + mr'^n) ■ d-y. 



We have 



2r> \ ^ rr^ , 2f^ ^ 9 , d d _ d d 



0K2 0K3 072 073 0X2 oKs 



and thus 



x{[X^„Xk,]) = -dx{X^,,XK,) = -mr^ f , ^ + , J' 2 ) ^ 0- 



This shows that the commutator [X^-^ , X^i ] does not belong to the characteristic distribution which is 
therefore not integrable. □ 

Therefore, to obtain a true Hamiltonian formulation of the reduced equations of motion in the case 
where the rank of A is 3, one needs to work with the bracket {■, -iR^nks- 
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If A has rank 2. As mentioned before, this case has the Chaplygin sphere as a particular example. 
The analysis of the two brackets has been done in [19]. We include it here for completeness and to link 
it with clarity to other results of the present work. 

In view of the form of A given in item (ii) of Definition [H we can write A^A = E — esej and thus, 
according to (|2.8p . we get 

K = (I + mr'^E)n - mr^(ri • 7)7, 

which is precisely the expression for the angular momentum about the contact point for the Chaplygin 
sphere. 

The angular velocity Q, can be written in terms of K and 7 as 

n = {l + mr'Ey'K + mr' ( „„ ^ " ";;:'^^''7^, ^ ) (I + mr'E)-'j, 

\ ||7||^ — mr^7 • (11 + mr^£/) -^7/ 

so both the equations ()2.10p and the Hamiltonian T-L-ji are well defined on M'^ x R'^. 

In this case, the two brackets for the system for functions f,gG C°°(M^ x M^) are given by 

None of the above brackets satisfies the Jacobi identity but it is preferable to consider {•, •}^j^ank2- The 
reason is that this bracket is conformally Poisson with conformal factor 

(/?(7) = I7I |2 - mr-2 (7 • (I + mr2S)-i7). (2.14) 

This important observation was first made in j6]. The characteristic distribution of {•,-}aank2 thus 
integrable. The generic leaves are the level sets of the Casimir functions Ci(K,7) = K • 7 and 
C2(7) = IItIP- Another important feature of this bracket is that it is twisted Poisson (in the sense of 
\24\ I29j) as will be shown in Section [5.31 (Theorem [5T7|) . 

On the other hand, similar to Proposition 12.21 we have 

Proposition 2.3 (|19|). The characteristic distribution of the almost Poisson bracket {•, •}Rank2 defined 
in ()2.13p is not integrable. 

This can be shown exactly as we did for Proposition 12.21 Therefore if the rank of A is 2, just as in 
the case of rank 3, a Hamiltonian formulation of the reduced equations can only be obtained if we work 
with the bracket {•,-}^ank2- However, in this case one needs to multiply the bracket by a conformal 
factor. This can be interpreted as a time reparametrization, see the discussion in Section [4.41 

If A has rank 1. Taking into account the form of A given in item (Hi) of Definition [1] we have 
A^A = 6363^ and thus, in view of (|2.8p . we get 

K = m + mr^(n -7)7. 
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The expression for the angular velocity fi in terms of K and 7 is 




so again, both the equations ()2.10p and the Hamiltonian T-L-ji are well defined on R'^ x R'^. 
This time, the two brackets for the system are given by 




(2.15) 



for functions /,c/ E C~(R3 x R3). 

The properties of the brackets above are very similar to those obtained in the case where the rank 
of ^ is 2 except that the roles of {•, -jRanki and {•, -iH^nki are reversed. 

This time one can show that 



Proof. We have to prove that the scaled bracket on IZ defined as ip{-, -jRanki satisfies the Jacobi identity, 
i.e., 

Rankl ) h} Rankl Rankl ) /l} Rankl Rankl ) /2}Rankl = 

for all /i , /2 , /s G C"^ (M^ x M^) . In view of the derivation properties of the bracket, is enough to show the 
identity for the coordinate functions Kj, 7^. In our case, since {7i, 7j}Ranki = 0, it is immediate to check 
that the identity holds if two of the three functions are 7i's. A long but straightforward computation 
shows that the identity holds for the following three choices of functions /i = i^i,/2 = K2, = 7i; 
/i = Ki,f2 = K2,f3 = 73 and /i = i^i,/2 = K2, fs = K3. Since the definition of the bracket is 
symmetric with respect to the coordinate functions Ki,'ji, and since the Jacobi identity trivially holds 
if two of the three functions /i, /2, /s are equal, all of the other cases are either trivial or analogous. □ 

Hence, the characteristic distribution of {•, -jRanki is integrable and the generic leaves are again the 
level sets of the Casimir functions Ci(K,7) = K-7 and 6*2(7) = IItIP- If ^ill also be shown in Section 
15.31 that {•, - jRanki is twisted Poisson. 

On the other hand, analogous to Propositions 12.21 and 12.3) we have 

Proposition 2.5. The characteristic distribution of the almost Poisson bracket {•, - jRamji defined in 
()2.15p is not integrable. 

The proof is again similar. 

Thus, this time the Hamiltonian structure of the reduced equations can only be obtained with the 
bracket {•, - jRanki j again through the multiplication by a conformal factor that is interpreted as a time 
reparametrization. 



Proposition 2.4. The bracket {•, •} 



Rankl 



defined in ()2.15p is conformally Poisson with conformal factor 



ifii) = VIItIP + '^f'^ (7 ■ 
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If A has rank 0. In this case A is the zero matrix and the constraints are holonomic and can be 
seen as a conservation law for the standard free rigid body. We have K = Ifl and clearly the equations 
(I2.10p and the Hamiltonian T-L-ji are well defined on x R^. The two brackets are given by 

, _ J. fdf dg\ fdf dg dg df\ 

If V (T^ 2o^ fdf ^ dg\ f df dg dg df^ ^^'^^^ 
{/. SiRanko = -(K-mrn)- — X— -7- — X x — 



dK dKj ' \dK d-f dK dj J ' 

The situation is analogous to that of the case when the rank of vl is 3 but, once more, the roles of the 
brackets are reversed. While {•, -iRanko coincides with the Lie-Poisson bracket in the dual Lie algebra 
se(3)* (and hence satisfies the Jacobi identity), we have 

Proposition 2.6. The characteristic distribution of the almost Poisson bracket {^-jRanko defined in 
()2.16p is not integrable. 



The proof is identical to that of Proposition 12.21 

So in this case, the Hamiltonian structure of the reduced equations (|2.10p can only be seen by 
working with the bracket {•, -jRanko- 



2.4 Symmetries 

The reduced equations ()2.10p can be interpreted as the output of a reduction process that we now 
explain. We begin by noticing that the configuration space Q = SO (3) x can be endowed with the 
Lie group structure of the three dimensional euclidean transformations SE(3). The group multiplication 
is given by 

(5l,Xi)(5r2,X2) = (5152,51X2 +Xi). 

Let H be the Lie subgroup of SE(3) defined by 

H = {{h,y) e SE{3) : teg = eg}. 

For matrices A satisfying any of the conditions of Definition [H it follows that HA = Ah whenever 
(h, y) G H. We consider the left action of on Q by left multiplication. The tangent lift of the action 
to TQ maps 

(h, y) : {g, X, uj, x) {hg, /ix + y, huj, h±) or {h, y) : {g, x, Q, x) ^ [kg, /ix + y, Q, h±), 

depending on the trivialization of S0(3) that one is working with. Notice that the Lagrangian C given 
by ()2.5p is invariant under the lifted action. Moreover, since h commutes with A for any {h,y) € H, 
the constraint ()2.2p is also invariant. 

The momenta (K, p) are geometrically interpreted as coordinates on the fibers of the (trivial) 
cotangent bundle T*Q. The constraint space A4 G T*Q is determined by the condition p = mrAgfl, 
so the triple (5,x, K) € S0(3) x x specifies a unique point in A4. Reciprocally, any point in M. 
can be represented by a triple (5,x, K). 
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By invariance of the Lagrangian and the constraints, the Ufted action of H to T*Q leaves M. 
invariant and, therefore, restricts to M.. The restricted action is free and proper so the orbit space 
TZ := M./H is a smooth manifold. The reduced space TZ can be identified with xR^; the projection 
p : — 7> 7^ is given by 

p(5,x,K) = (7,K), (2.17) 

where 7 = ^^^63 E S^, and is a surjective submersion. The conditions /ie3 = 63 and hA = Ah, that 
are satisfied for (/i, y) G ensure that the above mapping is well defined (in particular notice that K 
is invariant). The reduced equations on Tt are precisely ()2.10p when restricted to the level set ||7|| = 1. 
Notice that TZ inherits the (trivial) vector bundle structure S'^ xM^ — >■ S'^ from M.. 

In this sense, the entries of 7 should be considered as redundant coordinates for the sphere and 
the entries of K as coordinates on the fibers of IZ. Notice that, for any j = 0, ... ,3, both brackets 
'iRankj and {•, •Inj^nkj rcstrlct to the level set |[7|[ = 1 since 6*2(7) = IItIP is a Casimir function. 

2.5 Kinematics and integrability of the constraint distribution 

The constraint distribution on Q defined by equation (|2.2p has fundamentally different properties 
according to the rank of the matrix A satisfying the conditions of Definition [TJ On one extreme we 
have the case where ^ = and the distribution is integrable (the 3-dimensional integral leaves are 
given by S0(3) x {x} for x € M^). As mentioned before, in this case the constraints are holonomic and 
the problem reduces to the classical free rigid body problem (the center of mass of the body x remains 
constant in our inertial frame). 

The extreme opposite case occurs when rank A = 3. In this case the corresponding distribution 
is completely nonholonomic or bracket- generating, see e.g. [27]. By Chow's theorem, any two points 
in the configuration space Q can be joined by a curve {g{t),x{t)) satisfying the constraints. Thus, at 
least at the kinematical level, there are no restrictions on the values of x. 

The cases where the rank of A is 1 or 2 lie in between the situations described above. If the rank 
of A = 2, the third component X3 of x remains constant during the motion. This is in agreement with 
our observation that the Chaplygin sphere problem is a particular case of this type of constraints - the 
sphere rolls on a horizontal plane X3 = const. In this case, the constraint distribution is non-integrable 
but is nevertheless tangent to the foliation of Q by 5-dimensional leaves defined by constant values of 
X3. 

Finally, for the case where the rank of A equals 1, the first two components xi, X2, of x remain 
constant during the motion. The body goes up or down along the X3 axis at a speed that is proportional 
to its angular velocity about this axis. This time the constraint distribution is non-integrable but 
tangent to the 4-dimensional leaves given by constant values of xi and X2- 

Without going into technical definitions, we simply state that the degree of non-integrability of the 
constraint distribution increases with the rank of A, passing from an integrable distribution if A = 
to a completely nonholonomic distribution if rank A = 3. It is interesting to see how this correlates 
with the need of a gauge-transformation to Hamiltonize the problem (Remark 15. 5p . 
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3 Geometric Setting 



This section is concerned with the basic features of (almost) Dirac structures [12], with focus on the 
regular case, as well as their gauge transformations As we will see, these geometric structures 

provide the setup that gives rise to the different brackets introduced in Section [2j Although we will 
be mostly interested in the geometry of bivector fields, our discussion is presented at the general level 
of (almost) Dirac structures, as they provide the framework in which gauge transformations are most 
natural. 

3.1 Dirac and almost Dirac structures 

A Dirac structure on a manifold P is a subbundle L of the Whitney sum TP T*P such that 
(i) L is a maximal isotropic subbundle of TP © with respect to the pairing (•, •) given by 
{{X,a),iY,(3)) =a(y)+/3(X), for (X,a),(y,/3) eTP®T*P. 

{a) r(L) is closed with respect to the Courant bracket defined on T{TP © T*P) given by 

l{X,a),{Y,P)\ = {[X,Y],£xP -Wda), 
for {X,a),{Y,p) eT{TP®T*P), i.e., lr(L), r(L)l C r(L). 

The underlying manifold P is sometimes referred to as a Dirac manifold. 

Let pri : TP © j'*p be the projection onto the first factor of TP © T*P. A Dirac structure 

L on the manifold P carries a Lie algebroid structure with anchor pri\i : L — )• TP and bracket 
given by the Courant bracket |-,-] restricted to r(L). It follows that, for a Dirac structure L, the 
distribution pri{L) C TP is integrable, i.e., P can be decomposed into leaves O such that for each 
X (z P, TxO = pr\{Lx)- If pri(L) has constant rank (i.e., pri{Lx) C T^P has the same dimension for 
all X & P), then we say that L is a regular Dirac structure, and pri{L) defines a regular foliation. Just 
as a Poisson manifold P is the disjoint union of its symplectic leaves, each leaf of a Dirac manifold P 
carries a presymplectic form. 

Examples Let be a closed 2-form and vr be a Poisson bivector field, and consider the maps fi^ : 
TP T*P given by n^X) = n{-,X) = -ix^ = for X G TP and vr" : T*P TP as in Definition [2i 
Then 

Ln := graph{n^) = {{X, a) eTP® T*P : ix^ = -a} 

and 

:= graph{7r^) = {{X, a) £TP® T*P : 7r«(a) = X} 

are Dirac structures. Note that pri identifies with TP as Lie algebroids. Similarly, L-j^ can be 
naturally identified with T*P, and the Lie algebroid structure induced on T*P by L^r has anchor 
vrt* : T*P TP, and bracket 

[a,/3]^ = i?^tt(„)(/3) - £^^^)ia) - (i(7r(a,/3)) = £^tt(a)/3 - 

This bracket is uniquely characterized by [df,dg]Tr = d{f,g} and the Leibniz identity. 
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If a subbundle L of TP®T*P satisfies (z) above, but not necessarily (ii), then L is called an almost 
Dirac structure. Condition (ii) is called the integrability condition. We say that L is a regular almost 
Dirac structure when the distribution pri(L) on P has constant rank. Notice that this distribution 
might not be integrable in the almost Dirac case. 

Examples If is an arbitrary 2-form or vr an arbitrary bivector field, then their graphs Lq and 
Ltt are almost Dirac structures. The failure of the integrability with respect to the Courant bracket 
of Lq and L^^ is measured by dO, and ^[7r,7r], respectively. For = graph{7r'^), the distribution 
Tr^{T*P) = pri(L^) is generally non-integrable. If it has constant rank we call the almost Poisson 
structure regular. Note that the bracket [•, -j^ defined as in (j3.ip is R— bilinear, skew-symmetric and 
satisfies the Leibniz identity. However, in general, vr" does not necessarily preserve the bracket; instead, 
(see e.g. [8]), 

7r«([a,/3]^) = [7r«(«),7r«(/3)] -iw/3[vr,^], for a, /? G ^^(P). (3.18) 

Note that an almost Dirac structure L on P is of the form = graph{T:^) for a bivector vr if and 
only if 

rPnL = {0}, (3.19) 

and L is of the form = graph{Q^) for a 2-form (7 if and only if T*P f] L = {0}, see |12] . Another 
example of an almost Dirac structure that will be very useful for our purposes is given by L C TP©T*P 
defined as 

L := {{X, a) eTP® T*P : X e F, ixftlp = -alp}, (3.20) 

where F C TP is a subbundle, O is a 2-form on P and • \f denotes the point-wise restriction to F. If 
the subbundle F is an integrable distribution and is closed, then L is a Dirac structure. 

Proposition 3.1. The following statements hold: 

(i) There is a one-to-one correspondence between regular almost Dirac structures L C TP (BT* P and 
pairs {F,Qp), where F is a regular distribution on P and £ r(/\^P*). 

(ii) Let F C TP be a regular distribution on P. Given a section F*), there exists a 2-form 
Q on P such that Q\f = Qp. 

Proof, (i) Let L C TP (BT*P be a regular almost Dirac structure with distribution F := pri[L) C TP 
on P (not necessarily integrable). Consider the section ^2^? in T{/\^ F*) given, at each x G P, by 

^f(x)(X(^), Y'(^)) = -a(^.)(y(^)), for X,Y e T{F) such that (X(^),Q(^)) G L^^. 

It is a straightforward computation to see that il^? is well defined, i.e., it is independent of the choice of 
a. Conversely, given a regular distribution P on P and ili? G r(/\^ P*), we may define the subbundle 
L C TP © j'*p as the pairs (X, a) such that X G P and \x^F = —0(\f- 

{ii) Let W C TP be a regular smooth distribution such that it is a complement of P on P, i.e., 
r2;P = Fx ® Wx for each x G P (e.g., can be chosen to be the orthogonal complement of Fx with 
respect to a Riemmanian metric). The 2-form on P can be defined by 

a{X,Y) = ^lF{XF,YF), 
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for X,Y £ X{P) such that X = Xp + Xw and Y = Yf + Y-^, where Xf,Yf G r(F) and Xv^,Yw G 
r(VV). Differentiabihty of $7 follows from its definition and the smoothness of F and W. 

□ 

Corollary 3.2. Given a regular almost Dirac structure L, there exists a 2-form Q on P and a regular 
distribution F C TP such that L is written in the form ()3.20p . 

Remark 3.3. Note that the 2-form Q € i7^(P) is not uniquely defined and, in general, there is no 
canonical choice for it. o 

Given a subbundle F C TP, we say that a section in r(/\^ F*) is nondegenerate if it is nonde- 
generate as a bilinear form on F at each point. It follows from (j3.20p that Ker(r2ir) = L n TP, and as 
a consequence of condition ()3.19p we obtain 

Corollary 3.4. Let L he a regular almost Dirac structure and {F,Qf) the pair associated to it in the 
sense of Proposition HOI Then ilp is nondegenerate if and only if L is the graph of a hivector field vr. 
Explicitly, the relation between {F^VLp) and vr is 

7r''(a) = —X if and only if ix^F = 0-\f^ 

where X G r(F) and a G Q}{P). 

Following notation of Definition [21 if {•,•} is the bracket associated to the bivector field vr in the 
above Corollary, then {f,g} = nF{Xf,Xg), for all f,ge C°°{P). 

3.2 Twisted Poisson and twisted Dirac structures 

Poisson structures may be viewed as encoding integrability in two levels: first, the characteristic dis- 
tribution 7r^(T*P) C TP is integrable, i.e., tangent to leaves; second each leaf carries a nondegenerate 
2-form that is closed (and this leads to the Jacobi identity). Twisted Poisson structures are special 
types of almost Poisson structures that retain the integrability of 7r^{T*P) but allow the leafwise 2-form 
to be non closed. These objects turn out to be related to Hamiltonization. We start with the more 
general notion of twisted Dirac structures. 

Consider a closed 3-form <j) on P, and define the <j)-twisted Courant bracket [29] as follows: 

[(X,a),(y,/3)1^ = {[X,Y],£x/3-iYda + ixAY^), (3.21) 

for {X,a) and (Y, (3) in T{TP®T*P). Now, a subbundle L oiTP(BT*P is a 4>-twisted Dirac structure 
|29j if L is maximal isotropic with respect to (•, •) and the integrability condition 

[r(L),r(L)i^cr(L) 

is satisfied. 

As in the ordinary case, a twisted Dirac structure L on P induces a Lie algebroid on L given by 
the anchor map prili and the bracket [•, •]0|r(L)- Therefore pri{L) is an integrable distribution on P. 
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Examples If is any 2-form on P, then = graph{Cl ) is ((ir2)-twisted Dirac. One may check that 
a bivector field vr on P such that L^^ = graph^-ir'^) is </>- twisted Dirac verifies (see [29]) 



^[7r,7r] =7r« (3.22) 

The fohowing result gives more examples: 

Theorem 3.5. Let L C TP T*P be a regular almost Dirac structure such that pri{L) C TP is an 
integrable distribution on P. Then, there exists an exact 3-form (p with respect to which L is a (p-twisted 
Dirac structure. 

Proof. Let F := pri{L) C TP and Op be the section in T{/\^ F*) associated to L given by Proposition 
13. 1[ Since F is integrable, defines a 2-form f^ci on each leaf O where Fx = T^O at each x G P. By 
Corollary 13.21 there exists a 2-form 0, on P such that LqO, = ^lo where lq : O ^ P is the inclusion. 
We assert that L is a ((ir2)-twisted Dirac structure. In fact, for {X,a) and {Y, f3) in r(L), 

1{X, a), {Y, /?)l(do) = {[X, Y], £xf3 - iyda + ixArdn) e T{L) 

if and only if 

^[x,y]^\f = -{£xP - Wda + lx^Yd9.)\F■ 
Since F is an integrable distribution we obtain that, 

-{£xli -Wda + ixAYd^)\F = -£x{!i\F) + Wd{a\F) -'iXAYd{Vt\F) 

= £x^Y^F ~ iydixf^F ~ ^x/\Yd^F — \x,y]^f- 

which completes the proof. □ 
Remark 3.6. 

(i) Note that if L is a twisted Dirac structure then L is also twisted with respect to any closed 
3-form 6' such that id) — 6') vanishes on the leaves. 



(ii) There is no canonical choice for the 3-form (p given in Theorem 13.5 



Tvi^isted Poisson bivectors. 

Bivector fields tt such that = graph{Ti'^) is a (/)-twisted Dirac structure are called (p-twisted Poisson 
bivectors |24[ I29j. i.e., tt verifies condition ()3.22p . We are especially interested in these kind of bivector 
fields since, as we will see, they appear naturally in the examples of nonholonomic systems introduced 
in Section [2j If {•,•} is the bracket given by the (/>- twisted Poisson structure vr, then relation ()2.1ip 
becomes 

{/, {g, h}} + {g, {h, /}} + {/i, {/, g}} + <p{Xf,Xg,Xh) = 0, 

for f,g,h E C°°{P) and Xf = {•,/}. So the failure of the Jacobi identity is controlled by the closed 
3-form (p. 
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Example If vr'' : T*P — t- TP is an isomorphism, then vr is twisted and (j) = dVL where is the 
nondegenerate 2-form associated with vr, i.e., Q!' o ti'^ = Id where as usual Q!'{X) = —\x^- Other, less 
trivial examples, will be presented ahead in Corollary 13.71 

Let TT be a </)-twisted Poisson structure on the manifold P and [•, -JyT be the bracket on T*P given 
by ()3.ip . Note that vr^ does not preserve this bracket. However, using (jS.lSp and ()3.22p we obtain 

[7rS(a),7r''(/3)] = vr" ([a, /3]^ + , 

for 1-forms a, /? on P. The i;^-twisted Courant bracket induces a modification of the bracket (j3.ip via 
the identification of T*P and Ljr, 

[a,/3]0 = £^tt(a)/3 - M(/3)da + i7rtt(a)A7rtt(/3)'A, 

such that (r*P, [•,-]0,7rf) is a Lie algebroid [29j (see also [8j). The characteristic distribution 'k\T*P) 
defines an integrable distribution on P (that may be singular). Each leaf O of the corresponding 
foliation of P is endowed with a non-degenerate 2-form Q,o that is not necessarily closed. If vr is 
(/)- twisted, then <Klo = ''o'^^ where to ■ O ^ P is the inclusion. 

Important examples of twisted Poisson structures are contained in the following Corollary of The- 
orem [33) 

Corollary 3.7. Let n be a bivector field on P with an integrable regular characteristic distribution. 
Then, there exists an exact 3-form (p on P with respect to which vr is (p-twisted. 

Remark 3.8. Mechanical Example. It is shown in [18] that the (semi-direct) product reduction of 
the Veselova system yields a regular conformally Poisson bracket on the reduced space. It follows 
that its characteristic distribution is integrable and thus, by Corollary \2>.1\ it is also twisted-Poisson. 
This is a first example of a nonholonomic system whose reduced equations are formulated in terms 
of a twisted-Poisson bracket. Other examples (related to the motion of a rigid body with generalized 
rolling constraints) are made explicit in Section [5.3i o 

Remark 3.9. An interesting question, that remains to be answered, is to give a characterization of 
almost Poisson brackets possessing an integrable characteristic distribution that is non-regular. o 

Regular conformally Poisson bivectors. 

An interesting class of almost Poisson structures admitting an integrable characteristic distribution is 
given by conformally Poisson structures. Recall from Section [5] that they are bivector fields vr for which 
exists a strictly positive function ip G C°^{P), such that (p-rr is Poisson. A conformally Poisson manifold 
(P, vr) is the disjoint union of conformally symplectic leaves. 

Note that this property is stronger than asking for (P, vr) to be a Jacobi manifold since a conformal 
factor implies the global existence of a function such that ipn is Poisson, while in Jacobi manifolds the 
factor ip may be only locally defined. 

From Theorem 13.51 we observe that any regular bivector admitting a conformal factor is also a 
twisted Poisson bivector. The following Proposition explains the relation between these two properties. 
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Proposition 3.10. Let n be a regular conformally Poisson bivector field on P with conformal factor 
ip. Let Q, € ri^(P) be as in Corollary Then any closed 3- form verifying ()3.22p for vr coincides with 

—dip Ail. on the leaves. 

Proof. Since vr admits a conformal factor ip S C°^{P), tiien [tTjVt] = ^^<^ A vr. On the other hand, 
if $7 is the 2-form on P associated to vr given by Corollary 13.21 then for 51,(72 G C^{P) we have 
n{Xg-^,Xg2) = TT{dgi,dg2). Thus, for 51, 32, 53 S C°°{P) we obtain 

^[iT,TT]{dgi,dg2,dg3) = A7r{dgi,dg2,dg3) = -^dip A n{Xg^, Xg^, Xg^). 

Then we conclude that any closed 3-form cj) satisfying ()3.22|) coincides with {-^dip A il) on the leaves. 

□ 

Let (P, vr) be a regular conformally Poisson manifold. A 2-form on P satisfying the conditions 
of Corollary 13.21 verifies that LqQ is conformally symplectic on each leaf O. However, may not 
necessarily be conformally closed. 

3.3 Gauge transformations. 

In this section we will consider a natural action of the abelian group of 2-forms on P on the almost 
Dirac structures on P. This action is given by gauge transformations of almost Dirac structures by 
2-forms, and it was introduced in [29j. 

More precisely, consider an almost Dirac structure L in TP (BT*P. A gauge transformation by the 
2-form P is a map tb : L ^ TP ® T*P, given by Ts((X,a)) = {X,a + ixB) for {X,a) G L. The 
subbundle tb{L) of TP T*P given by 

TBiL) = {{X,a + ixB) : (X,a)GL} 

is an almost Dirac structure. If the 2-form B is closed and L is Dirac then tb{L) is again Dirac. Thus, 
the 3-form dB is what determines the integrability with respect to the Courant bracket. It is a direct 
computation to see that if L is a (/>-twisted Dirac structure, then the gauge transformation of L by the 
2-form B is {(j) — (iP)-twisted Dirac (see e.g. [29]). 

If Li and L2 sue almost Dirac structures on P and there exist a 2-form B on P such that tb{Li) = 
L2, then we say that Li and L2 are gauge equivalent or gauge related. 

Note that a gauge transformation does not modify the distribution pri[L). So, for a Dirac struc- 
ture L, the foliation associated to L will be the same as the one associated to tb{L). However, the 
presymplectic form on each leaf is modified by the pullback of B to the leaf. If P is a regular almost 
Dirac structure determined by the pair (P, 0^) in the sense of item (i) of Proposition 13. 1| a gauge 
transformation by the 2-form B corresponds to the operation: 

TB : {F,nF) ^ {F,nF-B\F). (3.23) 

Theorem 3.11. Any two regular almost Dirac structures Li and L2 are gauge related if and only if 
pri{Li) = pri{L2). 
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Proof. It remains to prove the "only if" part of the statement. Let us denote F := pri{Li) = pri{L2) 
and let and f2|, be the 2-sections associated to Li and L2 respectively (Proposition 13. II (i)). Define 
the section Bp G T{f\^{F*)) by Bp := - nj, and let B G n'^{P) such that B\f = Bp (Prop. O 
(ii)). We claim that r^Li = L2. In fact, if {F,Q,p) is the pair associated to the almost Dirac structure 
tbLi, then by equation ()3.23p . 

= — Bp = i^'p. 

Since the sections associated to r^Li and L2 coincide, by Proposition 13. II (i) we conclude that r^Li = 
L2 which means that Li and L2 are gauge related. 

□ 

We are especially interested in gauge transformations of almost Poisson structures. Consider the 
almost Poisson manifold {P, vr) and a 2-form B on P. Then, the gauge transformation of L.^ := 
graph{'K^) is tb{Lt^) = {{X,a + ixB) G TP T*P : X = 7r''(a)} which does not necessarily 
correspond to the graph of a new bivector vr^. A necessary and sufficient condition for this to happen 
is that 

TB(L,)nrp = {0}, 

which is equivalent to the fact that the endomorphism Id + i?'' o vr" : T*P — ?> T*P is invertible [29]. 
Indeed, if such a bivector field vr^ exists, then, in view of (|3.3|) . for any 1-form a on P we have 

TB ((vr«(a) , a)) = (7r«(a) , a + = ([n^f (a + , a + . 

Thus, TT^ is characterized by the condition 

(7r«)«(a + i,«(„)i?) =^tt(a), (3.24) 

and we can write 

(^S)tt = 7r«o(Id-S^o^«)-i. (3.25) 
In particular, if vr and tt^ are non-degenerate bivector fields, equation (j3.25p is equivalent to 

((^^)tt)-i = (vr«)-i-i?^ 

Therefore, any two 2-forms on a manifold are gauge related. This is not necessarily the case with 
bivector fields. A necessary condition is that their characteristic distributions coincide. In view of 
Theorem I3.11[ this condition is also sufficient if such distributions are regular. 

Although gauge related bivectors have the same characteristic distribution, their Schouten brackets 
[tTj-tt] and [7r^,7r^] may not coincide. The following Proposition makes this precise. 

Proposition 3.12 ([29]). // a (p-twisted Poisson bivector vr is gauge related with another bivector tt^ 
via the 2-form B, then vr^ is {(j) — dB) -twisted. That is, 

l[n'',n^] = {n^)KcP-dB). 

In particular, if is Poisson, then vr^ is {—dB) -twisted. 
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4 Applications to Nonholonomic Systems and Hamiltonization 



In what follows, we will analyze the geometry of nonholonomic systems in the framework presented in 
the previous section. We introduce gauge transformations of the bracket describing the nonholonomic 
dynamics in order to study the process of Hamiltonization. 



4.1 Nonholonomic systems 

A nonholonomic system consists of an n-dimensional configuration manifold Q with local coordinates 
q G f/ C M", a Lagrangian £ : TQ — )■ R of the form £(q, q) = ^^(q)(q, q) — ^(q), where ^ is a kinetic 
energy metric on Q and y : Q ^ 1^ is a potential, and a regular non-integrable distribution V C TQ 
that describes the kinematic nonholonomic constraints. In coordinates, the distribution T> is defined 
by the equation 

e(q)q = 0, (4.26) 

where e(q) is a /c x n matrix of constant rank k where k < n\s the number of constraints. The entries 
of e(q) are the components of the M'^-valued constraint 1-form on Q, e := e(q) dq. 

The dynamics of the system are governed by the Lagrange-D'Alembert principle. This principle 
states that the forces of constraint annihilate any virtual displacement, so they perform no work during 
the motion. The equations of motion take the form 

d fdC\ dC r r ^ . . 

( — \ = /i'^e(q). (4.27) 



dt \dqj dq 

Here fi : TQ — )• M.^ is an M'^-valued function whose entries are referred to as Lagrange multipliers. 
Under our assumptions, it is uniquely determined by the condition that the constraints ()4.26p are 
satisfied. The equations (|4.27p together with the constraints (|4.26p define a vector field Y^^^ on T> whose 
integral curves describe the motion of the nonholonomic system. A short calculation shows that along 
the flow of Y^, the energy function Ec := ^ • q — £, is conserved. 

The above equations of motion can be written as a first order system on the cotangent bundle T*Q 
via the standard Legendre transform, Leg : TQ — >■ T*Q, that defines canonical coordinates (q, p) on 
T*Q by the rule Leg : (q, q) i— > (q, p = dC/di\). The Legendre transform is a global diffeomorphism 
by our assumption that Q is a. metric. 

The Hamiltonian function, %: T*Q ^ M., is defined in the usual way n := Eco Leg-^. The 
equations of motion (j4.27p are shown to be equivalent to 

and the constraint equations ()4.26p become 

e(q)^ = 0. (4.29) 

The above equation defines the constraint suhmanifold A4 = Leg(D) C T*Q. Since the Legendre 
transform is linear on the fibers, 7W is a vector sub-bundle of T*Q that for each q Q specifies an 
n — k vector subspace of T*Q. 
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Equations ()4.28p together with ()4.29p define the vector field on A4, that describes the motion 
of our nonholonomic system in the Hamiltonian side and is the push forward of the vector field by 
the Legendre transform. The vector field is defined uniquely in an intrinsic way by the equation 

= + /^^^*^)' (4-30) 

where Qq is the canonical symplectic form on r*Q, l : A4 ^ T*Q is the inclusion and t : T*Q ^ Q is 
the canonical projection. The constraints (|4.29p and their derivatives are intrinsically written as 

X„h € C := TM n T, (4.31) 

where T is the distribution on T*Q defined as := {v G T{T*Q) : {t*€,v) = 0}. Denote by Qm the 
pull-back of Qq to Ai, i.e. '■= '-*f^Q- The following Proposition is of great importance for our 

setup of the equations of motion as an almost Hamiltonian system. 

Proposition 4.1 (\34:\ [2]). The distribution C on AA defined by (|4.3ip is regular, non-integrable, and 
the point-wise restriction of to C, denoted by ilc; is non- degenerate. 

The non- integr ability of C is a direct consequence of the non- integr ability of T). One shows that the 
rank of C is 2 (n — k) and that along Ad we have the symplectic decomposition 

Tm{T*Q) = C®C^'^, (4.32) 

where denotes the symplectic orthogonal complement to C. 

Since r*e vanishes on C, by restricting (j4.30p to C and denoting Hm '■= ^ C°°(A^), the 
equations of motion can be written in the appealing format 

'^X^S^c = {d-HMk, (4.33) 

where {d'HM)c is the point- wise restriction of dV-M to C. The above equation uniquely defines the vector 
field and is central in our treatment; with this in mind, we collect the data of the nonholonomic 
system in the triple {A4 , > T~iM ) • 

Even though ()4.33p defines the vector field uniquely, and resembles a classical Hamiltonian 
system, notice that since the distribution C is non-integrable, then Qc is a section in /\^ C* ^ A4 (not 
a 2-form). 

Let {^A,^lc,^^M) be a nonholonomic system. For every / € C°°{AA), let Xf denote the unique 
vector field on M with values in C defined by the equation 

ix^nc = {df)c, (4.34) 

where {df)c denotes the point-wise restriction of df to C. The vector field Xf defined by equation 
()4.34p is called the the (almost) Hamiltonian vector field associated to f. 

Since Qc is nondegenerate, by Corollary 13.41 there is a unique bivector field vr^h on Ai associated 
to the pair {C,il,c)j that is ^T^^^{a) = —X if and only if ij^^c = ct|c- On exact forms we have, for 
/ G C^{M), 

ixMc = {df)c if and only if ^L(d/) = -^/, (4.35) 
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which is consistent with notation of Definition [2j The bracket {•, - jnh on functions on Ai associated to 
the bivector TTnh describes the dynamics in the sense that 

X„,(/)(m) = Xnj^if)im) = {f,nMUirn) for ah / G C^iM). 

It follows from ()4.35p that the characteristic distribution of the bracket {•, ■}nh is C. Since C is non- 
integrable then {■, is an almost Poisson bracket that does not satisfy the Jacobi identity. 

Remark 4.2. If the constraint distribution V were integrable, the same would be true for the distri- 
bution C. Let A/" C 7W be a leaf of the corresponding (regular) foliation of A4 (i.e. Cx = T^M for all 
X G N). In view of ()4.32p the submanifold M is symplectic. Therefore, in this case, our construction 
of {•, -jnii coincides with the usual construction of the Dirac bracket on each leaf N of the foliation of 
M (see |131 112j). Hence, in this case, the Jacobi identity holds. o 

Inspired by equation (j3.20p and our data, it is natural to define the almost Dirac structure L^^ on 
M by 

L^^:={{X,a) eTM®T*M : X € C, ix^mIc = -a\c}, (4.36) 
as considered (up to a minus sign) in [32l |33] , and subsequently in [23] . 

Proposition 4.3. Let 7r„h be the bivector field defined in (j4.35p and {^-Inh the corresponding bracket. 
The following statements hold: 

(i) The almost Dirac structures -/^tt^^ := graph{Tr^-^i^) and L^ii given in (j4.36p coincide. 

{ii) The almost Poisson bracket {•, -jnh coincides with the classical almost Poisson bracket for non- 
holonomic systems defined in |30[ [25] . 

Proof, (i) It is immediate since both almost Dirac structure are defined by the same pair (C, ^c)- 
{ii) The bracket on M for nonholonomic systems introduced in |30[I25] was shown in [9] to be given 

by 

{f.g} = ^M{VYj,VYg), lor f,geC°°{M), (4.37) 

where V : T_m(T*Q) — )■ C is the projector associated to the symplectic decomposition ()4.32p . and Yj 
the free Hamiltonian vector field on the symplectic manifold {T*Q, Qq) defined by i^Og = df, where 

/ is an arbitrary smooth extension of / to T*Q. 

It is easy to check that along M one has i-py-^c = {df)c, so VYj' coincides with the almost 

Hamiltonian vector field Xj defined by equation ()4.34p . Therefore, for any f,g(z C°°{A4), 

{f,g} = nM{VYf-,rYg) = nc{Xf,Xg) = -df{7rl{dg)) = {f,gU. 

□ 

The second item of the above Proposition should not be surprising since the expression ()4.37p is the 
nonholonomic version of the Dirac bracket (see discussion in [21\ lU]). Hence, its description naturally 
falls in the ambit of almost Dirac structures as described above. As a consequence of the above 
Proposition, we can equivalently describe our nonholonomic system with the triple (A^, TTnh, "Hx). 
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Definition 3. Let TTnh, ^A4) be a nonholonomic system. 

1. The bivector field vTnh on A4 given by ()4.35p is called the nonholonomic bivector field and the 
bracket {■, is called the nonholonomic bracket. 

2. We say that an almost Dirac structure L describes the dynamics of the nonholonomic system if 
the pair (— G ^(L). 



4.2 Gauge transformations of the nonholonomic bracket 

The main idea of using gauge transformations in our setting is that it opens the possibility to modify 
the geometric structure on Ai that describes the dynamics. 

Consider the nonholonomic system {^A,7^^^,,'HJ^) and continue to denote L,rnh = 9^(^pf^{'^ih)- The 
gauge transformation of 7r„ij associated to a 2-form i3 on gives 

TB{L^„J = {iX,a + ixB)eTM®T*M : 7rl{a) = X}. (4.38) 

First of all we are interested in knowing when the pair {—X^i-^dJiM) is a section of tb{L-,^^_^^). On 
the other hand, we would also like to know whether the almost Dirac structure tb{Lt^^^J corresponds 
to the graph of a bivector field or not. 

If the 2-form B on M verifies ij^ = 0, then from equation ()4.38p we see that the pair 
(— Xnh, dHj^) belongs to r(rB(L^^^)). Moreover, in view of (|3.23|) . the gauge transformation of Tr^h by 
the 2-form B has the form 

TBiL^J = {{X,a) eTM(BT*M : X e C, ixi^M - B)\c = -a\c}. (4.39) 

Thus the equations of motion (|4.33p are equivalently written as 

ix A^c - Be) = {dTiMk, 

where Bq is the point- wise restriction of to C. 

Therefore, as a particular case of Corollary 13.41 we observe that if the section 0,c ~ Be is non- 
degenerate then the gauge transformation of 7r„ij associated to the 2-form B is again a bivector field 
vTj^^ . It follows from ()3.25p that the non-degeneracy of — B on C is equivalent to the invertibility of 
the endomorphism (Id — o tTj^i,) on T*A4 . 

Definition 4. Let (P, vr) be an almost Poisson manifold with a distinguished Hamiltonian function 
H £ C°°(P). Given a 2-form B on P, the gauge transformation of vr associated to the 2-form B is said 

to be a dynamical gauge transformation if 

(i) 'ixhB = 0, where Xjj is the (almost) Hamiltonian vector field associated to H and 

(zi) TB{graph(Tr^)) corresponds to the graph of a new bivector vr^, i.e., the endomorphism (Id — S^'ovr") 
on T*P is invertible. 
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Of course we are interested in dynamical gauge transformations of the nonholonomic bracket where 
the distinguished Hamiltonian function is Hm ^-iid the corresponding (almost) Hamiltonian vector field 
is Xnh- 

Note that if vr is regular, by equation (|3.23p . the gauge transformation defined by B is determined 
by the restriction Bp oi B io F where F := 'k^{T*P). Then condition (i) of the above Definition is 
equivalent to ix^Bp = 0. 

Remark 4.4. The definition of an affine almost Poisson bracket for a nonholonomic system made 
in [19] corresponds to a dynamical gauge transformation of the nonholonomic bracket by a 2-form 
B = —l*Qq where is a semi-basic form 2-form on T*Q. The proof is analogous to that of item (ii) 
of Proposition 14.31 In this case, the hypothesis that is semi-basic implies that the condition (ii) of 
Definition H] is satisfied (see Proposition 14.61 below) . o 

After this discussion, we observe that it is more appropriate to describe a nonholonomic system 
by the triple {A4 , ^, Hm ) where ^ is the family of bivector fields that are related to vr^h through a 
dynamical gauge transformation. Notice that C is the characteristic distribution of any bivector field 
in 5- Thus, the (almost) Hamiltonian vector fields defined by the corresponding brackets satisfy the 
nonholonomic constraints. It follows that the bivector fields in the family 5 are almost Poisson bivectors 
in a "strong" sense since the non-integrability of C prevents them from being twisted or conformally 
Poisson. Our interest in considering this big family of brackets relies on reduction. In the presence 
of symmetries, the bracket that Hamiltonizes the reduced equations may arise as the reduction of a 
member of 5 that is not necessarily TTnh- 

Since the distribution C is regular we observe 

Corollary 4.5. [of Theorem 13. 1 1] . All bivectors with characteristic distribution equal to C are gauge 
related (in particular, gauge related to the nonholonomic bivector tt^i-,). 

We finish this section by discussing some cases for which the second condition in Definition |3] is 
satisfied. Recall that M. C T*Q is a vector bundle over Q. We have 

Proposition 4.6. If B is a semi-basic 2-form on Ai, then the gauge transformation of tt^y, associated 
to B corresponds again to a bivector field. 

Proof. The graph of vrfj^ is an almost Dirac structure corresponding to the pair (C, Qc) in the sense of 
Proposition 13.11 (see Proposition 14.31 («)). Thus, in view of (|3.23p . the gauge transformation of Tr^h is 
the almost Dirac structure corresponding to the pair (C, {^m ~ B)\c). It is shown in [19j that \i B a. 
is a semi-basic 2-form, then the point-wise restriction of {^m — B) to C is non-degenerate. Thus, by 
Corollarv 13.41 tb{Lt^^^) corresponds to the graph of a bivector. □ 

In fact this Proposition is a special case of the following result: 

Proposition 4.7. Let P ^ Q be a vector bundle and vr a regular almost Poisson bivector on P. If 
for all semi-basic 1-forms a on P the vector field ■k'^{q) is vertical, then the gauge transformation of t: 
associated to a semi-basic 2-form B corresponds again to a bivector. 



27 



Proof. Consider local bundle coordinates (q, p) € C/ x 1/ C x M™", on P such that q are local 
coordinates on the base manifold Q. Since 7r^(dq) is vertical and B is semi-basic we obtain 

B^oJ{dci) = 
B'oJ{dp) = 6(q,p)dq, 

where 6(q, p) denotes the m x n matrix with entries 6aj(q, p) = {B^ o n'^ (^dpa) , , for a = l,...,m, 
and j = 1, n. Thus, the matrix representation of the endomorphism (Id — B^ o vr^) on T*P is 

f Idnxn -fe(q,P)'^ \ 

xm J 

This matrix has full rank and hence (Id — B^ o vr") : T*P — )• T*P is invertible. □ 

Remark 4.8. It is interesting for future work to drop the condition (ii) in Definition H] that requires 
TB{graph{'K\^)) to define a bivector field. In this case, the family 5^ consists of all the almost Dirac 
structures that are gauge related to and that describe the nonholonomic dynamics. In this sense, the 
Hamiltonization of the problem is achieved if the reduction of a member of is a Dirac structure (but 
not necessarily a Poisson structure). This approach requires the consideration of a general reduction 
scheme for almost Dirac structures. However, we are unaware of any examples of nonholonomic systems 
that justify the need of such a general framework. o 



4.3 Reduction by a group of symmetries 

We now add symmetries to the problem and perform the reduction. Our interest from the point of 
view of Hamiltonization is to find a bivector field in the family whose reduction is either Poisson or 
conformally Poisson (see Section [4. 4p . 

Let G be a Lie group acting freely and properly on Q. We say that that G is a symmetry of the 
nonholonomic system if the lifted action to TQ is free and proper, and leaves the constraint distribution 
T> C TQ and the Lagrangian C : TQ M invariant. 

Denote hy ^ : G x T*Q ^ T*Q the cotangent lift of the action to T*Q. If G is a symmetry for 
our nonholonomic system, then ^ leaves both the constraint submanifold A4 and the Hamiltonian 
H : T*Q —?■ M invariant. We continue to denote by ^ the restricted action to A4. 

One can show that the tangent lift of ^ to TAi, preserves the distribution C and the section fi^. 
As a consequence, if G is a symmetry group for our nonholonomic system, then the action ^ preserves 
the standard nonholonomic bracket {•, -Inh- That is, for f,g(z C°°{A4), we have 

By freeness and properness of the action, the reduced space TZ := Ai/G is a smooth manifold and the 
orbit projection map p : M. ^ TZ \s a surjective submersion. Notice that TZ inherits a vector bundle 
structure from M over the shape space Q/G. Moreover, TZ is equipped with the reduced nonholonomic 
bracket {■, -jrod that is characterized by 

{/,5}rcd op(m) := {/op,5op}„i,(m) for m G and /, 5 G G°° (7^) . (4.40) 
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The corresponding bivector field will be denoted tt^^^. The reduced nonholonomic bracket describes 
the reduced dynamics in the sense that the nonholonomic vector field X^i, is p-related to the (almost) 
Hamiltonian vector field X^^^ = {•,7i'jz}^cd associated to the reduced Hamiltonian Ti-ji defined by the 
condition T-Lm = ^7?. ° P- 

Remark 4.9. The reduction of nonholonomic systems performed by Bates and Sniatycki in [2J shows 
that it is possible to define a 2-form cj^^d on IZ which is non-degenerate along a distribution C C TTZ. 
The definition of C is given by C := Tp{U) where [/ := C n (C n V)^^ and V is the distribution on 
M. tangent to the orbits of G. In fact, the pair {C^lo^^^) is just the pair associated to the almost Dirac 
structure given by the bivector field -k^^^ in the sense of Corollary 13.41 o 

The analysis of the reduction of the bivector fields vr^^ in the family ^ that are related to vTnh by a 
dynamical gauge B follows from: 

Proposition 4.10. Let (-P, vr) he an almost Poisson manifold and B a 2-form on P such that the 
endomorphism (Id — B^ o n^) : T*P j'*p {g invertible. If the Lie group G, acting freely and properly 
on P, preserves the almost Poisson structure vr and leaves B invariant, then G preserves the bivector 
field TT^ obtained by the gauge transformation of n associated to B. 

Proof. In view of equation ()3.25p . we see that (vr^)" is a composition of invariant maps and we conclude 
that TT^ is invariant as well. □ 

As a direct consequence of this Proposition we have: 

Proposition 4.11. IfG is a symmetry group of the nonholonomic system ( , Tr^h, 7^_a4 ) and B is aG- 

invariant dynamical gauge on A4, then is G invariant. In particular there is a reduced bivector 

field iT^.^^B on the reduced space TZ that determines a well defined bracket {•, ■}^.^^b on TZ satisfying 

{/,5LdB op(m) = {/op,5op}f,(m), for f,g E C^{n). 

Moreover, the reduced bracket {■,-}^^^b also describes the reduced dynamics in the sense that Xy^^ = 

There is a good reason why we did not denote the reduced bivector field tt^^^b by vr^^^ and that is 
that in general the reduced bivector fields Tr^ed and vr^^^s need not be gauged related. We shall see this 
explicitly with the analysis of our mechanical examples (Remark 15. 5p . 

Remark 4.12. Since the almost Dirac structure given in (|4.36|) is the graph of a bivector (see 
Proposition 14.31 (i)), then the reduction of 1/^1, as an almost Dirac structure is simply the reduction 
of the bivector n^i, in the classical way. The same observation is valid for the reduction of the almost 
Dirac structure rs(L^^^) defined in (|4.39p since tb{Lt^^^) is also the graph of a bivector (($7>i — B) is 
non-degenerate on C). o 
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The G-Chaplygin case. 



If the reduced bivector fields vr^cd and tt^^^b happen to be everywhere non-degenerate, then they are 
gauge-related. This is the scenario that one finds after reduction of external symmetries of G-Chaplygin 
systems (see [II])- These systems are characterized by the property that the tangent space to the orbits 
of the symmetry group exactly complements the constraint distribution on the tangent space TQ of 
the configuration manifold Q. 

In this case there exist unique non-degenerate 2-forms Q^^^ and 0,^^^b on TZ satisfying 

7rLol7L=Id, nl^Bonl^B=ld, 

where as usual, = — ix^^rcd and r2^^^s(X) = —ix^^^^B for all X E X(JZ). In particular, the 

reduced equations can be written as: 

ix^^^^rod = ixn^^rcdB = dHu, 

which gives different almost symplectic formulations of the reduced equations (compare with the results 
in [20]). The bivector fields Tr^cd and tt^^^iB are gauge related via the 2-form iJ^cd := ^rcd — ^vcd^ on TZ. 
Moreover, the 2-forms Q^^^ and 0,^^^b satisfy 

(p*aed)|c = and (p*aedB)lc = i^M - B)c, 

and thus iJ^ed verifies that (/9*-Brod)|c = Be- In other words, we have the following commutative diagram 

Ln u — > ^TT^ 

^71-rod ^ ^^r^dB ■ 

We stress that one does not have such a commutative diagram in more general situations, see Remark 
[531 



4.4 Hamiltonization 

We continue with the notation of the previous sections and suppose that G is a symmetry group for 
our nonholonomic system. The solutions of the reduced equations on the reduced space IZ = M./G are 
the integral curves of the reduced vector field Xy^^ and preserve the reduced Hamiltonian l-L-ji- 

The issue of Hamiltonization in our context concerns answering the question of whether the reduced 
vector field X-^^^ on TZ is Hamiltonian. Our candidates for the Hamiltonian structure come from the 
reduction of the (invariant) bivector fields that belong to the family 

It turns out that the above condition is too restrictive. We relax it by asking that the vector field 
X-y^ can be rescaled by a basic^ positive function 99 : 7^ ^ M in such a way that the resulting vector 
field ipX-fi^ is Hamiltonian. 

*We use the term basic with respect to the fibered structure of 71 inherited from M. That is (p — (p o r where 
r : 7?. — > Q/G is the bundle projection and if> : Q/G ^ M.. 
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In view of Proposition 14. 1 1| for any G-invariant bivector field vr^^ belonging to ^, the rescaled vector 
field ipXq-i^ satisfies 

{if7^,^,Bf{d'Hn) = -vX'H^- (4.41) 
Hence, we are interested in finding a bivector field vr^.^^s satisfying 

[(^^,,dfl,'/^vr^eds] =0, (4.42) 

that is, we want to find vr^.^^s conformally Poisson (see Definition . 

Definition 5. If there exist an invariant bivector field vr^^,, belonging to ^ and a strictly positive, 
basic function 99 : 7^ — )• M such that ()4.4ip and ()4.42p hold, we say that the nonholonomic system 
is Hamiltonizable. Moreover, we say that the reduced equations are Hamiltonian in the new time r 
defined by dr = ^dt (see discussion below). 

In Section [5] we will extend the discussion of Section [2] and show that all generalized rolling systems 
are Hamiltonizable. The table ()5.48p in Section [5] shows how different scenarios of the Hamiltonization 
scheme described above are realized according to the rank of the matrix A. 



Time reparametrizations 

It is common in the literature to interpret the rescaling of the vector field X-^^ by the basic positive 
function ip as a. nonlinear time reparametrization. One argues as follows, let c{t) S 7?- be a flow line of 
X-^i^ (i.e. ^{t) = X'^^{c{t))). Introduce the new time r by integrating the relation 

dr = — , , dt. 

(^(c(t)) 

Since (p > 0, the correspondence between t and r is one-to-one and one can express t as a function of 
T. The curve c(r) := c(t(r)) is checked to be a flow line of (fXy^^ (i.e. ^(t) = (/3(c(r))X^^(c(r))). 
This interpretation of the rescaling is quite subtle. The definition of r depends on the particular flow 
line c(t), so different initial conditions induce different reparametrizations. It is therefore not possible 
to interpret the time rescaling as a "global" operation. This contrasts with the natural procedure of 
multiplying the vector field X-}{^ by the positive function if. 

Remark 4.13. One might wonder why we only care about basic and not arbitrary functions ip : 
TZ M+. A detailed answer to this question would involve a careful study of the structure of the 
equations of motion that would gear us away from the main subject of this paper. We refer the reader 
to |15] where one can find a very good discussion on the Hamiltonization of G-Chaplygin systems. We 
simply mention that physically, the fact that (p is basic means that it is independent of the momentum 
variables and only depends on the (reduced) configuration variables. Thus, the time reparametrization 
changes the speed at which the trajectories are traversed depending on the position of the system but 
independently of the velocity itself. It is shown in [15] how in order to obtain Darboux coordinates 
for the reparametrized system, one can keep the same (reduced) configuration variables but should 
rescale the momenta by ^. Since (p is basic, the rescaled momenta continue to depend linearly on the 
velocities. o 



31 



Measure preservation 

Hamiltonization is strongly related to the existence of invariant measures. Suppose for simplicity that 
we are dealing with a G-Chaplygin system. As mentioned before, in this case the bivector field tt^^^b 
is everywhere non-degenerate and the reduced equations can be written as: 

where ^^^^b is the non-degenerate 2- form on TZ induced by tt^^^b . It follows that the scaled vector field 
^X-j-i^ satisfies 

Hence, Hamiltonization in this setting amounts to finding a positive function 99 such that the 2-form 
^^rodS is closed (which under our hypothesis, is of course equivalent to ()4.42p ). Suppose for a moment 
that this is the case so iJZ, ^^^^^b) is a symplectic manifold. It follows from Liouville's theorem that 

the vector field LpX-^i^ preserves the symplectic volume {^^^^^b) where m = 2 dim 7^. Therefore, 
the volume form (^)™'^^(f^i.cds)'^ is preserved by the vector field X-^^. 

The above argument shows that a Hamiltonizable G-Chapligyn system possesses an invariant mea- 
sure. One might wonder if the reciprocal statement is true, namely, if any G-Chaplygin system with 
an invariant measure is Hamiltonizable. The celebrated Chaplygin's reducing multiplier theorem 
demonstrates that the answer is positive if m = 2. For m > 2, a characterization of the systems for 
which this is true is an open problem. Interesting examples where this holds for arbitrary values of 
m have been found by Fedorov and Jovanovic in the study of the multidimensional Veselova problem 
|15j . See also the discussion in [H] where a candidate for the conformal factor ^ is given under the 
hypothesis that there exists a preserved measure, and [17j where a set of coupled first order partial 
differential equations for the multiplier ^ are given. 

In the case where the nonholonomic system is Hamiltonizable but the corresponding bivector field 
t^^^^b is degenerate at some points in IZ, one can repeat the above argument to conclude that the 
reduced system preserves a measure on every leaf of the symplectic foliation of TZ corresponding to the 
Poisson bivector field (fTr^^^B . However, this does not imply the existence of a smooth invariant measure 
on TZ (and this is the motivation for Fernandez, Mestdag, and Bloch [T7j to talk about Poissonization) . 
An example of this situation is given by the reduction of the Chaplygin sleigh (see the discussion in 
|18j ) that exhibits asymptotic dynamics that contravene the existence of a global invariant measure. 
The problem of the existence of a global invariant measure in this case is most naturally attacked by 
considering the modular class of the Poisson manifold (7^, (/JTr^^^s), see [35] . 

5 Back to the Examples: Hamiltonization and Integrability 

In this section, we analyze the generalized rolling systems presented in Section [2] using the geometric 
framework that was developed in the previous sections. In subsection 15.11 we provide the geometric 
interpretation for the brackets {•,-}Rankj and {^-jpankj presented in Section [2j In 15.21 we consider the 
Hamiltonization and integrability of generalized rolling systems in detail and finally, in l5.3l we explicitly 
show that the brackets {•, -jRanki and {•, ■}'^^n^2 are twisted Poisson. To our knowledge, this is the first 
time that the appearance of such structures is made explicit in the field of nonholonomic mechanics. 
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5.1 The geometry of the rigid bodies with generalized rolling constraints 

We begin by computing the nonholonomic bracket for the motion of a rigid body subject to generalized 
rolHng constraints as introduced in Section [2j 

Nonholonomic bracket via the non-degenerate 2-section 

Consider again, as in Section [21 the motion of a rigid body in space subject to a generaUzed rolhng 
constraint as in ()2.2p . That is, the constraint relates the linear and the angular velocities of the body 
X = rAu = rAgft, where the matrix A satisfies any of the conditions of Definition [1] and cj, is the 
angular velocity written in space and body coordinates, respectively. 

Recall that the configuration space for the system is Q = SO (3) x M^. Denote by A (respectively, p) 
the left (respectively, right) Maurer-Cartan form on S0(3). Upon the identification of the Lie algebra 
5o(3) with by the hat map ()2.3p we think of A and p as valued 1-forms on S0(3). For a tangent 
vector Vg € Tg S0(3) we have 

^ = p{9){vg), n = X{g){vg), 

where u (respectively, Ct) denotes the angular velocity vector written in space (respectively, body) 
coordinates as discussed in Section [2j 

The Maurer-Cartan forms A and p are related by X{g) = g^^p{g) and satisfy the well-known 
Maurer-Cartan equations 

dp=[P,P], d\=-[\,X], 

where [•, •] is the commutator in the Lie algebra. For the rest of the section we will use three dimensional 
vector algebra notation in our calculations with differential forms and vector fields. In our convention, 
the scalar product of differential forms should always be interpreted as a wedge product (and is thus 
anti-commutative!). The Maurer-Cartan equations take the form 

dp = X p, d\ = X (5.43) 

where "x" denotes the standard vector product in M'^. 

The constraint distribution D, defined by the generalized rolling constraints, can be expressed in 
the terminology of subsection 14.11 as the annihilator of the M^-valued 1-form e on Q given by 

e = dx — rAp = dx — rAgX. 

We consider the (global) moving co-frame {A, dx} for T*Q that defines fiber coordinates (M, p) 
in the following sense. A co-vector Uq G T*Q is written uniquely as = M • A -|- p • dx, for a 
certain (M,p) € M'^ x M^. The Legendre transform Leg : TQ — ?• T*Q associated to the kinetic energy 
Lagrangian ()2.5p is defined by the rule: 

M = in, p = mx. 

Physically, p is the linear momentum of the body while M is the angular momentum of the body about 
the center of mass written in body coordinates. 
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In order to deal with the constraints, it is more convenient to work with the global moving co- 
frame {A, e} for T*Q. We denote by (K, u) the fiber coordinates defined by this co-frame. Putting 
K-A-|-u-e = M- A-|-p-dx implies 

K = M rg^A^p, u = p. 

Along the constraint submanifold Ad = Leg(P) we have p = mrAg^l so 

K = m mr'^g^ Agn, 

which is the expression for the kinetic momentum obtained in (|2.8|) . Notice that is a vector bundle 
over Q and that K is a natural coordinate for the fibers of 7W. In what follows we will use the 
components of (7,x, and K, as redundant coordinates on M.. 

Denote by = {X\',X^,X^) the moving frame of S0(3) that is dual to A = (Ai,A2,A3). The 
components of are the left invariant vector fields on SO (3) obtained by left extension of the canonical 
basis of M^. Along the points of the constraint subbundle Ad, the non-integrable distribution C defined 
in ()4.3ip is given by 

C = span|x^ + r/A^|-, A|. (5.44) 

The canonical 2-form fig on T*Q is given by 

VLq = -d{m • A p • dx) = -d(K • A + p • e) 
= A • dK — K • dA — p • de -|- e • dp, 

where "•" denotes the usual scalar product in R^. 

To compute de we use the identity dg = g\ where " denotes the hat map (j2.3p . Using the 
Maurer-Cartan equations (|5.43p we get 

de = -rA{dg)X - rAgdX = -rAg{\ x A) - rAgdX = rAgdX. 

Therefore, 

JIq = A • dK - K • dA - p • {rAg dX) + e ■ dp 
= X-dK-{K + rg'^A^p) -dX + e-dp. 

Let L : M ^ T*Q denote the inclusion. Since p = mrAgCl along Ai, we have 

Qm ■■= t^*{^Q) = A • dK - (K + mr^g'^A^AgQ) • dA + t*(e • dp). 

Since e vanishes along the non-integrable distribution C, we get the following expression for the restric- 
tion Qc of i^*{^q) to C: 

nc = X-dK-(K + mr'^g'^A^AgQ,) ■ dX. 
With this expression for ilc we are ready to show: 
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Proposition 5.1. The nonholonomic bracket {•,-}nh on A4 for the generalized rolling system is given 
in the redundant coordinates {^jj, x^, A";}, i,j,k,l = 1,2,3, for M by 

{xi, = r{Ag)u, {ffij, ^Jm. = 9ik, {Kh Kj}^^ = -^(K + mr'^{g^ Agn))i , 

with all other combinations equal to zero. In the above formulas the Einstein convention of sum over 
repeated indices holds and e^j denotes the alternating tensor, that equals if two indices are equal, it 
equals 1 if {i,j,k) is a cyclic permutation of {1,2,3), and it is equal to —1 otherwise. The entries of 
{g,x) G Q are denoted by {gij,Xk) and those o/K by Ki. 

Proof. We will rely on the identity ()4.35p that characterizes the associated bivector field Tr^h in terms 
of Vtc ■ Contraction of by the elements in the basis of C given in (|5.44|) gives 

ifxL , r„T^T^)f^c = dK - (K + mr'^ g'^ A^ AgO,) x A , i^^c = -A, 
where we have again made use of the Maurer-Cartan equations ()5.43p . It follows that 

\x.^+r9TAT^^(K+mr^gTATAgn)x^)^C = dK, i_,,^g^J^C = rAgX = d^\c. 

Therefore, according to (|4.35p we get 



nlidK) = -X^ - rg^A^^ + (K + mr^g^A-^AgQ) x 

ax oK 

7rL(A) = vrJi,(dx) = rAg— . 



(5.45) 



In addition, for any canonical vector ej € M'^ we have d{g ej) = {g ej) x A, so 

d 

T^L{d{g~^ei)) = (g^^ei) x — . 

The proof follows using the above formulas and recalling that for any f,g& C°°{A4), we have {/, g}nh = 
-dfiTrlidg)). 

□ 



Finally, we state without a formal proof that the nonholonomic vector field on A4 is given by 

X^^ = n-X.'^ + rAgft- — + {Kxft)-—. (5.46) 

ax aK 

The above expression can be shown by taking into account the equations ()2.10p . the constraint ()2.2p . 
and the definition of fl, or, alternatively, by computing the the almost Hamiltonian vector field Xy^j^ = 
—T^ihidTiM) corresponding to the Hamiltonian that coincides with expression ()2.9p . The latter 
approach requires one to write ft in terms of K and g as was done in Section [2] for the different values 
of the rank of A. 
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The gauge transformation of the nonholonomic bracket 

We will now construct a gauge transformation of the nonholonomic bracket in the sense of subsection 
14.21 We are interested in describing the same dynamics so we look for a 2-form B that defines a dy- 
namical gauge transformation as introduced in Definitional In our case, the distinguished Hamiltonian 
is TiM that has Xnh as its associated almost Hamiltonian vector field. 

Following [19J , we consider the bi-invariant volume form on S0(3) oriented and scaled such that 
v{X^ , X2 , X^) = 1. We consider the natural extension of as a 3-form on Q = S0(3) x M'^. Denote 
by G Q'^{T*Q) the 3-form given by iy = v*^ where v : T*Q ^ Q is the canonical projection. We can 
write = iA • (A X A). 

Let be the 2-form on A4 given by 

where, as before, l : A4 ^ T*Q is the inclusion. Note that i? is a semi-basic 2-form on A4 that 
vanishes upon contraction with the nonholonomic (almost) Hamiltonian vector field X^^. Therefore, 
by Proposition 14. 6t we can perform a dynamical gauge transformation of the nonholonomic bivector 
field TTnh by the 2-form B to obtain another bivector field vr^ that also describes the dynamics of our 
problem. 

Using the Maurer-Cartan equations ()5.43p and the expression ()5.46p for X^ih, we obtain 

B = -mr'^fl ■ dX. (5.47) 

To compute the bivector field vr^^ associated to the gauge transformation we use equation (|3.24|) . For 
an arbitrary one-form a on we have 

Setting a equal to A and dx and using (j5.45p and (j5.47p we obtain 

i^^fW = (vrfj«(dx) = rAg^. 

Similarly, putting a = dK. and noticing that 

M (dK)^ = -ixL-B = mr'^n x A, 

nh^ 

we deduce 

(vrfj«(dK) = -X^ - r/A^|- + (K + mr\g^A^Ag - E)n) x A, 

where E denotes the 3x3 identity matrix. 
The above formulas imply 

Proposition 5.2. The gauged nonholonomic bracket {•, ■}^^^ on M, associated to the bivector field vr^^, 
is given in the redundant coordinates {gij,Xk,Ki), i,j,k,l = 1,2,3, for Ai by 

{xi,Ki}^^ = r{Ag)u, {9ij,Ki}^^ = -e% ff.fc, {K^, K,}f, = -4(K + mr'^^g' A^ Ag - E)n)i , 

with all other combinations equal to zero. 
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Reduction of the symmetries 

Recall that the Lie group introduced in section 12.41 acts on the configuration space Q and that its 
lift to TQ leaves both the Lagrangian and the constraints invariant. From the discussion in section 
14.31 (and the regularity of the action) it follows that the reduced space TZ := M./H \s equipped with a 
reduced bracket {■, - jrcd determined by condition ()4.40p and that describes the reduced dynamics. 

We are now ready to give the geometric interpretation of the bracket {•, -jRankj introduced in section 
[2] for the different values of the rank of A. 

Theorem 5.3. The reduced bracket {•, - jjcd on TZ is precisely the restriction of the bracket {■, ■}^^nkj 
(defined in section\^ to the Casimir level set ||7|| = 1, for the different values j = 0, 1, 2, 3, of the rank 
of A. 

Proof. Recall from section [23] that the reduced space TZ can be identified with S'^ xM^ with redundant 
coordinates (')',K). Therefore, it makes sense to compare the two brackets on the Casimir level set 
II7II = 1 of the space (7,K) G x M^. 

Moreover, from the expression of the projection p : M ^ TZ given by (|2.17p . and condition ()4.40p . 
it follows that the reduced bracket of the (redundant) coordinate functions (7, K) can be computed 
using the formulas obtained in Proposition 15.11 (notice that 7 = (531, (732, 533))- 

The proof is completed by considering the particular form of A for the different values of its rank 
given in Definition [U and by writing the bracket {/i,/2}rcd of arbitrary functions /i,/2 G C°°{TZ) in 
terms of the derivatives ^ and ^ using Leibniz rule. □ 

We now turn to the study of the reduction of the gauged nonholonomic bracket {•, -j^. First of 
all notice that the 2-form B that defines the gauge transformation is written in (j5.47p in terms of left 
invariant objects on S0(3). Since the symmetry group H acts by left multiplication on the S0(3) 
factor of Q, it follows that B is invariant under the cotangent lifted action. Therefore, in accordance 
with Proposition 14. IH the gauged bracket {■,-}nh drops to TZ where it defines the bracket {•, -j^^dS that 
determines the dynamics. As usual, the corresponding bivector field on TZ will be denoted by tt^^^b. 

In analogy with Theorem 15.31 we have 

Theorem 5.4. The reduced bracket {■,-}^^^b on TZ is precisely the restriction of the bracket {■ , -y^ii^nkj 
(defined in section\B(j to the Casimir level set ||7|| = 1, for the different values j = 0, 1, 2, 3, of the rank 
of A. 

The proof is identical to that of Theorem 15.31 except that one uses the formulas obtained in Propo- 
sition 15.21 

According to the discussion in Section [21 the following table summarizes the properties of the 
reduced bivector fields Tr^cd and tt^^^b for a generalized rolling system according to the different values 
of the rank of A: 
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Rank of A 





1 


2 


3 








Non-integrable 


Non-integrable 




Poisson 


Conformally Poisson 


characteristic 


characteristic 








distribution 


distribution 




Non-integrable 
characteristic 
distribution 


Non-integrable 
characteristic 
distribution 


Conformally Poisson 


Poisson 



Remark 5.5. Notice that the reduced bivector fields ir^cd and tt^^^b are not gauge related. Indeed, 
from the table above one sees that for any value of the rank of A, only one of the two bivector fields Tr^cd 
or vTj.^jjS has an integrable characteristic distribution. It follows from Theorem 13.111 that there cannot 
exist a gauge transformation between their graphs. o 

Remark 5.6. Recall from subsection 12.51 that as the rank of A increases, the constraint distribution is 
less integrable or "more nonholonomic" . The table (j5.48p above seems to suggest that it is appropriate 
to perform a gauge transformation by the 2-form B when the nonholonomic effects are more important, 
while the reduction of the standard nonholonomic bracket works better for weaker nonholonomic effects. 

o 



5.2 Hamiltonization and integrability of rigid bodies with generalized rolling con- 
straints 

According to the notion of Hamiltonization introduced in Section 14.41 (Definition [5]), and the table 
()5.48p . it immediately follows that the problem of the motion of a rigid body subject to a generalized 
rolling constraint is Hamiltonizable for any value of the rank of A. 

If the rank of A equals (respectively, 3) the reduced equations are Hamiltonian with respect to 
the bracket {•, ■}^^^ (respectively, {•, - jrcds)- Recall that in both cases the reduced dynamics correspond 
to classical rigid body motion (with modified inertia tensor I -|- mr'^E if the rank of A equals 3). 

If the rank of A equals 1 or 2, the analysis of the Hamiltonization is a bit more delicate but it also 
follows directly from Definition [S] and the table (|5.48|) . In the case rank A = 2, it follows that the 
reduced equations are Hamiltonian in the new time T2 defined by dT2 = -;^dt and with respect to the 
bracket (/32{t } rod s where 

(/?2(7) = \/l - mr2 (7 • (I -h mr'^E)-^j). (5.49) 

Note that (p2 is a basic function on TZ corresponding to the restriction of (j2.14p to the level set ||7|| = 1. 

Analogously, if the rank of A equals 1, the reduced equations are Hamiltonian in the new time ti 
defined by dri = -^^dt and with respect to the bracket ipi{-, - jrcd where 

(y9i(7) = y^l + mr2 (7-1-17). (5.50) 
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Integrability of the reduced equations 



In view of the Hamiltonization of the problem, the integrability of the reduced equations of motion 
()2.10p can be easily established using the celebrated Arnold-Liouville Theorem for classical Hamiltonian 
systems, see e.g. [Ij. 

Indeed, for any value of the rank of A, the reduced equations are Hamiltonian on TZ (after a time 
reparametrization if rank A = 1,2). Independently of the rank of A, the symplectic leaves Oa of the 
foliation of TZ correspond to the level sets Ci(K, 7) = K • 7 = a and can be shown to be diffeomorphic 
to the tangent bundle TS^ of the sphere (see the discussion in chapter 14 of [26] for the coadjoint orbits 
on sc(3)*). 

Once the value of a is fixed, the reduced equations ()2.10p can be seen as a two degree of freedom 
classical Hamiltonian system on Oa (again, after a time reparametrization if rank A = 1,2). These 
equations possess two independent integrals, the Hamiltonian H-jz, and F = K • K, whose joint level 
sets are compact in Oa- It follows from the Arnold-Liouville Theorem that these level sets are invariant 
two-tori and the dynamics are quasi-periodic on them (notice that the flow on the tori is rectilinear 
but not uniform if the rank of A is 1 or 2). 

The Arnold-Liouville Theorem also tells us that the reduced equations are integrable by quadratures 
(after the time reparametrization if the rank of ^ is 1 or 2). 

Finally, we state without proof that the reduced equations of motion ()2.10p preserve the measure 
ld{'y) a A dKi A dK2 A dK^ where a is the area form of the sphere S^, and the basic density /i : — > M 
is given by 

'1 if rank A = 0,3, 

— ^ if rank A = 1, 



where (pi,ip2 S C°°(S^) are defined in ()5.50p and (|5.49p respectively. 



— ^ if rank A 

¥'2(7) 



5.3 Twisted Poisson structures for rigid bodies with generalized rolling constraints 

In Section 13.21 we presented twisted Poisson structures which have been extensively studied in other 
contexts but not in mechanics. Now, we will show explicitly that twisted Poisson structures appear 
naturally in the study of nonholonomic systems. 



Rigid body with generalized rolling constraints of rank 2 

Here we show that the bracket {•, •}^j,„,j2! in addition to being conformally Poisson, is twisted Poisson. 
Note that this cannot be the case for the other bracket {•, •}Rank2 that describes the dynamics since, as 
shown in Section [2l its characteristic distribution is not integrable. 

Recall from the discussion in 12.41 that {■,-}aank2 should be considered as a bracket on the reduced 
space TZ = S'^ xM^ with redundant coordinates (7, K). The characteristic distribution of the bracket 
is integrable and the leaves Oa of the foliation are the level sets Ci(7,K) = 7 • K = a. By regularity 
and integrability of the characteristic distribution, it follows from Corollary 13.71 that the bracket is 
(ji-twisted. The value of the 3-form <j) is given in the following. 
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Theorem 5.7. The bracket {•, •}Hank2 defined in \2.13\) (that in particular describes the reduced dy- 
namics of the Chaplygin sphere for the appropriate choice of A), is a (j)-twisted Poisson bracket with 
(j) = —dB where 

B = mr'^{n--f)a, (5.51) 
and where a denotes the area form of the sphere \\^\\ = 1. 

Proof. The idea of this proof is to show that the bracket is gauge related to a Poisson bracket via the 
2-form —B. Thus, by Proposition 13.121 the bracket is (— d;B)-twisted Poisson. More precisely, we will 
show that {•, •}Rank2 is — ;S-gauge related with the bracket {•, -jRanko defined in ()2.16p and that coincides 
with the Lie- Poisson bracket on se(3)*. 

According to Theorem [231 we denote the bivector field associated to the bracket {•, ■y^^^^^2 '^rodS- 
Using ()2.13p one gets 

(vr,„,B)«(dK) = (K - mr\n ■ 7)7) >< ^ + ^ >< i^..,B)Hdl) = 7 x 

Next, notice that the 2-form B defined by ()5.5ip is written in the redundant coordinates (7,K) as 

B = ^mr^(i7 • 7) 7 • (^7 x ^7) . 

The bivector field tt^^^b and the 2-form B verify hypothesis of Proposition 14.71 and thus the gauge 
transformation of vr^.^^s associated to B is again a bivector field that we will denote it by t^^^^b ■ Relying 
on equation (|3.24p . one computes 

where, as usual, E denotes the 3x3 identity matrix. Using these expressions we deduce 

(dK) = K X A + 7 X A, (vrl.)«(d7) = 7 x A, 
that complete the proof. □ 



The conformal factor and the 3-form (p. In accordance with Proposition 13.10] since the bracket 
{■, •}^^jjk2 is both conformally Poisson and twisted Poisson, there is relationship between the conformal 
factor (p2 (given by (|5.49p ). and the twisting 3-form cp (defined in Theorem 15. 7p . 

We leave it to the reader to check that on the leaves Oa of the foliation of TZ corresponding to 
the bracket {•, •}^ank2i the 3-form (j) coincides with V' •= ^ dip2 A where 2-form i7 is given in the 
redundant coordinates (7, K) by 

n = ^{K- mr^{n ■ 7)7) • (d7 x ^7) - 7 • {dK x ^7) . 

This choice of Q, satisfies the conditions of Corollary 13.21 for the graph of the bivector field vr^.^^s 
corresponding to {•, •Inj.nka on TTZ © T*TZ. 
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Rigid body with generalized rolling constraints of rank 1 

A completely analogous analysis can be performed if the rank of the matrix A equals one. This time 
it is the bracket {•, -laanki that is both twisted and conformally Poisson. In analogy with Theorem 15.71 
we have 

Theorem 5.8. The bracket {•, - jaanki defined in i2.15\) . is a (j)-twisted Poisson bracket with (j) = dB 
with B given by expression ()5.5ip . 

The proof is the same to the Rank 2 case. The bracket {•, -jRanki is S-gauge related with the bracket 
{'i 'iRanko defined in ()2.16p and that coincides with the Lie-Poisson bracket on se(3)*. 
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1 Introduction 



It is well known that the equations of motion for a mechanical system with nonholonomic constraints 
do not arise from a variational principle in the usual sense (see e.g. [3]). As a consequence, they cannot 
be formulated as a classical Hamiltonian system. Instead, they are written with respect to an almost 
Poisson bracket that fails to satisfy the Jacobi identity. This formulation has its origins in [341 [29l |23] 
and others. 

On the other hand, after a symmetry reduction, the equations of motion of a number of examples 
allow a Hamiltonian formulation (sometimes after a time reparametrization) , in which case one talks 
about Hamiltonization^ (see [I2l[l6l[171[l9l[2ll[22l[2ll32]). 

In this paper we employ recent ideas from Poisson geometry to study the Hamiltonization phe- 
nomenon from a geometric perspective. The main tool in our appraoch is the concept of gauge transfor- 
mation of bivector fields (in the sense of Severa and Weinstein [33j), an operation that uses differential 
2-forms to modify bivector fields keeping their characteristic distribution (see Section [51 Definition [2} 
unchanged. 

Let us consider a nonholonomic system on a constraint phase space A4, equipped with the almost 
Poisson bracket {•,-}nh [34[ I29[ [23]. known as the nonholonomic bracket, and Hamiltonian function 
TiM- The (almost) Hamiltonian vector field X^h, defined by {-j-jnh and TImj governs the dynamics 
of the system. In this paper, we consider new brackets {•, -j^ obtained by gauge transformations of 
{■j 'Inh with respect to suitably chosen 2-forms B on Ai. Since gauge transformations do not modify 
the characteristic distribution, the (almost) Hamiltonian vector fields associated to {•, -j^ satisfy the 
nonholonomic constraints. If, in addition, B verifies ij^ ^ i? = 0, we say that it defines a dynamical 

difFerent meaning to Hamiltonization is given in where the authors study the unreduced system in connection 
with the inverse problem of the calculus of variations. 
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gauge transformation; the terminology reflects the fact that such gauge transformations do not affect 
the dynamics, in the sense that the (almost) Hamiltonian vector field defined by {•, -j^ and Hm is still 
Xnh- In this way, we distinguish a family J of almost Poisson structures that describe the dynamics of 
our nonholonomic system. 

The main motivation to consider the family 5 of almost Poisson brackets is to have a larger choice of 
structures to describe the reduced dynamics and hope to find one amongst them that Hamiltonizes the 
problem. More specifically, in the presence of a group of symmetries G, we consider the dynamics on the 
reduced space TZ := M./G. The reduction of the invariant brackets in 5 yields a collection of reduced 
brackets {•, -^^^^b on 7^, and any member of this collection describes the reduced dynamics in (almost) 
Hamiltonian form with respect to the reduced Hamiltonian T-Ltz- In particular, the reconstruction of 
the original dynamics on M is exactly the same regardless of the choice of bracket {•, -j^^dS on the 
reduced space. 

In our approach, the issue of Hamiltonization is formulated by the requirement that one of the 
reduced brackets on TZ, obtained as the reduction of an element in "^^ is conformally Poisson; in other 
words, there should exist a bracket {•, ■}^^^b and a positive function (p such that (p{-, •}rcdS is a Poisson 
bracket on TZ (i.e., it satisfies the Jacobi identity)^. The scaling of {-j-jredS by (p is dynamically 
interpreted as the time reparametrization dr = ^dt (an idea that goes back to [l2]). In particular, the 
Hamiltonization of the celebrated Chaplygin sphere problem [llj arises as the reduction of a dynamically 
gauged bracket {•, -j^ that differs from the standard nonholonomic bracket {•, -jnij. Our approach to 
Hamiltonization is summarized in the following diagram. 



Non- 
holonomic 
system 



mi 



almost Poisson 
description in terms of the 
nonholonomic bracket {■, -jnh 



dynamical gauge 
transformation 
by a 2- form B 



almost Poisson description 
in terms of a family ^ of 
almost Poisson brackets {•, -j^^ 

{witli the same cliaracteristic distribution) 



symmetry 
reduction 

almost Poisson description of the 
reduced dynamics in terms of the 
reduction of invariant members 
in g" denoted by {•, -j^^dS 

(tliese members may have different characte- 
ristic distribution and may not be gauge related) 



Hamiltonization 



find an element {•, -jredS 
in the reduced family 
that is (conformally) Poisson 

A key observation is that, although the brackets in ^ are all gauge related, they may have fundamentally 
different features after reduction. For example, depending on the choice of the 2-form i?, the charac- 
teristic distribution of the reduced almost Poisson bracket {•, -j^^dS rnay or may not be integrable. The 

■^In fact the function ip : TZ is not arbitrary. It is required to be basic with respect to the fibered structure of 

the reduced space Ti. 



symmetry 
reduction 



The reduction of {•, -Inh 
gives the known almost 
Poisson bracket {•, -jred 

(e.g. m) 



{■, 'Ij-Qfj is an clement 
of the reduced family 
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integrability of the characteristic distribution is of central importance, as it is a necessary condition for 
{•, - jrodS to be conformally Poisson (and, hence, for the system to be Hamiltonizable) . 

An important class of almost Poisson brackets possessing an integrable characteristic distribution 
is given by the so-called twisted Poisson structures, introduced in |26j and [33j. For these structures 
the failure of the Jacobi identity is controlled by a global closed 3- form (see Section I3.2|) . In this 
paper we show that any almost Poisson structure whose characteristic distribution is both integrable 
and regular is a twisted Poisson structure. As a consequence, we show that the reduced equations 
of certain nonholonomic systems (including the Veselova problem [35] and the Chaplygin sphere) can 
be formulated in terms of twisted Poisson brackets in the original physical time (prior to any time 
reparametrization) . To our knowledge, this establishes the first connection between twisted Poisson 
structures and nonholonomic mechanics. This should serve as a motivation to investigate the interplay 
between the rich geometrical properties of twisted Poisson brackets and the dynamical features of the 
corresponding (almost) Hamiltonian vector fields. 

To illustrate the method depicted in the diagram above, we introduce the problem of the motion 
of rigid bodies that are subject to generalized rolling constraints. These are nonholonomic constraints 
that relate the angular velocity u of the body and the linear velocity x of its center of mass in a linear 
way, i.e., x = Au: for a 3 x 3 matrix A satisfying certain properties. This type of constraints generalize 
the Chaplygin sphere problem. In fact, the constraints vary from completely nonholonomic if the rank 
of A equals 3, to (holonomic) classical free rigid body motion if A = 0. For the Chaplygin sphere, the 
rank of A equals 2, in which case we recover the Hamiltonization method of [21]. By allowing gauge 
transformations, we prove that the problem is Hamiltonizable independently of the value of the rank 
of A. Using the Hamiltonian structure of the reduced system, we also show its complete integrability 
in the sense of Liouville. 

For this explicit class of examples, the behavior of the reduced brackets, according to the rank of A, 
is illustrated in table (jl.ip below. In our notation, {•, ■}^^^ corresponds to the reduction of the classical 
nonholonomic bracket, while the bracket {•, ■},.^^b corresponds to the reduction of a bracket obtained 
by a dynamical gauge transformation by a specific 2-form B (defined in Section [5]) . 



Rank of A 




(free rigid body) 


1 


2 

(Chaplygin sphere) 


3 

(completely 
nonholonomic) 


{'i 'jrcd 


Poisson 


Conformally Poisson 
and Twisted Poisson 


Non-integrable 
characteristic 
distribution 


Non-integrable 
characteristic 
distribution 


{•7 •}red« 


Non-integrable 
characteristic 
distribution 


Non-integrable 
characteristic 
distribution 


Conformally Poisson 
and Twisted Poisson 


Poisson 



(1.1) 



Hamiltonization in context: The most interesting example of Hamiltonization concerns the Chap- 
lygin sphere. Even though the formulation and integration of the equations of motion by Chaplygin 
dates back to 1903 [TT|, the Hamiltonian structure of the reduced equations (after a time reparametriza- 
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tion) was only discovered in 2001 by Borisov and Mamaev [6]. Their result is all the more remarkable 
as it is in apparent contradiction to the assertions of \15\ [T6] . 

Recently, Jovanovic [25] proved that the multidimensional version of the Chaplygin sphere problem 
introduced in [181 ^^^^ integrable and Hamiltonizable when the vertical angular momentum is zero. 
This gives a partial solution to a problem that remained open for many years. His approach to prove 
integrability involves in a crucial way the Hamiltonization of the problem. Another important example 
where the integration of a nonholonomic system follows from its Hamiltonization is the multidimensional 
Veselova system treated by Fedorov and Jovanovic in [17j . We also mention the recent work of Ohsawa, 
Fernandez, Bloch and Zenkov |32j that explores the connection of Hamiltonization with Hamilton- 
Jacobi theory. 

The relationship between Hamiltonization and integrability may have been the original motivation 
for Chaplygin to consider the problem of Hamiltonization back in 1911 [T2]. In this work, Chaplygin 
proved the famous Chaplygin reducing multiplier Theorem that applies to the so-called G- Chaplygin 
systems. These are nonholonomic systems with the property that the tangent space to the orbits of 
a symmetry group G exactly complements the constraint distribution on the tangent space TQ of the 
configuration manifold Q. Stated in modern geometric terms, the theorem says that if the shape space 
Q/G is two-dimensional, and the reduced equations have an invariant measure, then they can be put 
in Hamiltonian form in the new time r defined by dr = ^dt. The positive function ^ : Q/G ^ M 
is known as the reducing multiplier^. There is a very neat interpretation of the multiplier ^ in terms 
of the invariant measure and as a conformal factor for an almost symplectic form that describes the 
dynamics, see [16l |T71 [22l [32] . This interpretation suggests that geometric methods may be useful to 
understand Hamiltonization in more general scenarios of nonholonomic systems with symmetry. The 
geometric approach to the study of nonholonomic systems has received a lot of interest in the last 
couple of decades and has its origins in the seminal paper of Koiller ^7] . 

Recently, Fernandez, Mestdag and Bloch [T9], derived a set of coupled first order partial differential 
equations for the multiplier ^ for G-Chaplygin systems whose shape space has arbitrary dimension. 
Even more, the set of equations found by the authors applies to general nonholonomic systems with 
symmetry that are not necessarily G-Chaplygin. This is done by writing the reduced equations of 
motion in Hamilton- Poincare-D'Alembert form as described in [31 H]. The issue of Hamiltonization 
is thus reformulated as a problem of existence of a solution for the aforementioned system of partial 
differential equations. 

Our approach to Hamiltonization contains the same degree of generality but is more intrinsic. 
Indeed, our requirement that one of the members {•, -j^^dS in the collection of brackets describing the 
reduced dynamics is conformally Poisson can be rewritten as 

[(/J7r^,dS,V'vr^eds] = 0, (1.2) 

that locally defines a set of coupled, partial differential equations for the multiplier ip. Here tt^^^b is the 
bivector field associated to {•, ■}^^^b and [•, •] is the Schouten bracket. Our approach is also broader than 
that of }19j since we have the freedom of selecting the bivector field tt^^^b in (|1.2p within the collection of 
brackets describing the reduced dynamics. We stress that this liberty is crucial to Hamiltonize certain 

■^We denote the reducing multiplier by ^ instead of ip to be consistent with our exposition which takes the Poisson 
rather than the symplectic perspective. 
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examples. As mentioned before, and as at is shown in table (jl.ip . the Hamiltonization of the Chaplygin 
sphere (by a "one-step reduction") depends on this freedom. The approach taken in O [19] involves a 
two-step reduction and the symplectic interpretation of the second reduction is delicate (see [22] ) . 

If the symmetries are of G-Chaplygin type, the bivector fields vr^^^s are everywhere non-degenerate 
and the equations of motion can be written with respect to the associated almost symplectic form 
ri^pjB. If a multiplier satisfying ()1.2p exists, then ^^^^b is conformally closed, (d(^^^rods) = 0), and 
one speaks of Chaplygin Hamiltonization [16]. The term Poissonization was introduced in [19j to refer 
to the case where the bivector field tt^^^b satisfying ()1.2p is degenerate. Their motivation to distinguish 
this case is to study the relationship between Hamiltonization and the existence of invariant measures 
for the reduced equations. To simplify our exposition and to treat the problem in a unified manner, 
we will not use their terminology. We give a discussion on the existence of invariant measures for 
nonholonomic systems admitting a Hamiltonization in this generality in subsection 14.41 

Outline of the paper: The paper is structured as follows. In Section [2] we introduce our motivating 
examples: rigid bodies subject to generalized rolling constraints. After writing down the reduced equa- 
tions of motion, we define two different almost Poisson structures for these equations according to the 
rank of the matrix A that defines the constraints. For each value of the rank, we show that only one of 
the brackets is Poisson or conformally Poisson, while the other one possesses a non-integrable charac- 
teristic distribution. The geometric interpretation and construction of these brackets is postponed to 
Section [5l after the necessary tools are developed in Sections [3] and HI 

Section [3] presents the geometric background needed for our purposes, with focus on the study of 
gauge transformations of bivector fields in the sense of [33j. Although our main interest is in almost 
Poisson brackets, we consider more general objects known as almost Dirac structures [IH], as they pro- 
vide the most natural setting for the definition and study of gauge transformations (Dirac structures 
have been also considered in connection with nonholonomic mechanics in [361 1371 [25] although the issue 
of Hamiltonization and the incorporation of gauge transformations are not treated in these works). 
Using a general description of regular almost Dirac structures in terms of their characteristic distri- 
butions (Proposition 13. 1|) . we show that any almost Poisson bracket possessing a regular, integrable, 
characteristic distribution is a twisted Poisson bracket. We also prove that any two regular almost 
Poisson (or, more generally, Dirac) structures defining the same characteristic distribution are gauge 
related (Theorem 13. lip . 

In Section [J] we make the connection between the geometric methods developed in Section [3] and 
nonholonomic mechanics. Using Proposition 13.11 and the framework of nonholonomic mechanics de- 
scribed in Bates and Sniatycki [2] we construct the nonholonomic bracket of [34 [ [29l I23| . Then we 
introduce the notion of dynamical gauge transformations for a nonholonomic system and define the 
family ^ of almost Poisson brackets describing the nonholonomic system (as in the diagram above). 
After discussing the reduction of these brackets in the presence of symmetries, we introduce our working 
definition of Hamiltonization. 

In Section[5]we resume the study of rigid bodies subject to generalized rolling constraints. We show 
that the brackets given in Section [2] to describe the reduced dynamics arise as a reduction of different 
members of the family 5 (Theorems 15.31 and 15. 4p . We then establish the Hamiltonization of the reduced 
equations in detail and we conclude their integrability. Finally, we consider the twisted nature of the 
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brackets that Hamiltonize the problem prior to time reparametrization when the rank of A is 1 or 2; 
in these cases, we also provide explicit expressions for the associated closed 3-form. 



2 Motivating Examples: Rigid bodies with Generalized Rolling Con- 
straints 

Consider the motion of a rigid body in space that evolves under its own inertia and is subject to the 
constraint that enforces the linear velocity of the center of mass, x, to be linearly related to the angular 
velocity of the body u, i.e., 

X = rAu. (2.3) 

Both vectors x and u> belong to and are written with respect to an inertial frame. The constant 
scalar r has dimensions of length and is a natural length scale of the system. The dimensionless, 
constant, 3x3 matrix A is given and satisfies certain conditions that are made precise in the following 
definition. 

Definition 1. The matrix A is said to define a generalized rolling constraint if it satisfies one of the 
following conditions according to its rank: 



(i) ^ = ^ J ^ , with C € S0(2), if rank^ = 3. 
{ii) A=l M , with C € S0(2), if rank^ = 2. 



.0 . 

[Hi) A = e^, if rank A = 1, where is the third canonical vector in M^, and T denotes transpose. 
{iv) A = if rank A = 0. 

The above conditions on A can be relaxed (see Remark 1 2 . 1 1 ahead) . However, for simplicity, we will 
assume that A has the form given by one of the items of the above definition. If A satisfies any of the 
the conditions of the above definition we say that (|2.3|) is a generalized rolling constraint. 

Our terminology is motivated by a particular example: the Chaplygin sphere. The problem, intro- 
duced by Chaplygin in 1903 [11], concerns the motion of a ball whose center of mass coincides with its 
geometric center that rolls on the plane without slipping. In this case, the matrix A is given by 

A 

and r is the radius of the sphere. 

The motion of the Chaplygin ball has been the subject of much research to our days. An important 
property is that (after a time reparametrization) the reduced equations can be given in a Hamiltonian 
structure [01 13 EI]- The geometry of the Hamiltonization of the problem is intricate. In order to study 
this phenomenon in a mathematically systematic fashion, we consider more general possibilities for the 
matrix A. 
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The crucial property of A that determines many of the dynamical and geometrical features of the 
problem is its rank. The Chaplygin sphere corresponds to the case rank A = 2. Another familiar case 
occurs when ranked = 0. In this case the constraint (|2.3p becomes x = which can be interpreted as 
a conservation law for the free system that states that the center of mass of the body is at rest in the 
inertial frame. The motion of the system reduces to that of the classical free rigid body. 

We will also consider the cases where the rank of A equals 3 and 1 which, to our knowledge, have 
not yet been considered in the literature. 

2.1 Generalities 

The configuration space for the system is Q = SO (3) x R^. Elements in Q are of the form q = (^jx) € 
SO (3) X M^. The vector x G is the position of the center of mass in space and the orthogonal matrix 
g specifies the orientation of the ball by relating two orthogonal frames, one attached to the body and 
one that is fixed in space. We will assume that the body frame has its origin at the center of mass and 
is aligned with the principal axes of inertia of the body. These frames define the so-called space and 
body coordinates respectively. 

Recall that the Lie algebra so (3) can be identified with equipped with the vector product via 
the hat map: 

■n = {rii^m^m) ^ "n = \ m o . (2.4) 

V -m m / 

Given a motion {g{t),:x.{t)) € Q, the angular velocity vector in space coordinates, u> € M^, and the 
angular velocity vector in body coordinates, € M^, are respectively given by 

u{t) = g{t)9'Ht), n{t) = g-Ht)g{t), 

and satisfy fi = g^^u. It will be useful to write the constraint (|2.3|) in terms of the body angular 
velocity as 

± = rAgn. (2.5) 
The kinetic energy of the rigid body defines the Lagrangian C : TQ — )• R by 

£(5,5,x,i) = ^(m)-fi + ^||i||2, (2.6) 

where "•" denotes the Euclidean scalar product on R'^, m is the mass of the body and the 3x3 diagonal 
matrix I is the inertia tensor with positive entries Ii,l2,l3- 

Remark 2.1. It is not hard to see that if the space axes are rotated by an element h € S0(3), the 
Lagrangian C is invariant and the constraint ()2.3p is rewritten as 

X = rh^^Ahu. 

Therefore, the conditions for A given in Definition [1] can be relaxed by allowing conjugation by matrices 
h G S0(3). o 
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2.2 The equations of motion 



Let p = mx be the linear momentum of the body. In accordance with the Lagrange-D'Alembert 
principle, the constraint forces must annihilate any velocity pair (x,ri) satisfying ()2.5p . Therefore, the 
equations of motion are given by 

p = /X, m = in X n - rg-^A^n, (2.7) 

where "x" denotes the vector product in M"^ and the multiplier /i € is determined uniquely from 
the constraint (|2.5p . 

Differentiating ()2.5p and using gQ = we find = mrAgfl. Thus, the second equation in ()2.7p 
decouples from the first to give 

m = m X n - mr'^g-^A^Ag ft. (2.8) 

In principle, this equation should be complemented with the reconstruction equation g = gCl. It will be 
shown ahead that it suffices to consider the evolution of the Poisson vector 7 := ^^^63 that represents 
the vector 63 written in body coordinates. A direct calculation gives 

7 = 7x0. 

The decoupling in ()2.7p is due to the presence of symmetries that will be discussed in detail in Section 
12.41 Once this equation is solved for (g, ft), we obtain p = mrAgfl that follows from (|2.5p . 

We introduce the kinetic momentum K G by 

K ■.= m + mr^g-'^A^Agit. (2.9) 
This definition of the kinetic momentum allows us to define the (reduced) Hamiltonian 

-H7^ = ^(K•f^), (2.10) 

which coincides with the kinetic energy on the constraint space A4. 
A direct calculation using (j2.8p gives our final set of equations 

K = Kxn, 7 = 7xr2. (2.11) 

To understand why the above equations define a closed system for (K,7) E x M^, and to understand 
their structure, it is useful to perform a separate study for different values of the rank of the matrix A. 
This will also show that the Hamiltonian "Hti can be considered as a function of K and 7. 

2.3 A pair of (almost) Poisson brackets for the equations of motion 

For each value of the rank of A we will give two different brackets that define the equations of motion 
(|2.1ip . with respect to the reduced Hamiltonian Tin. In general, these brackets are almost Poisson, 
i.e. they do not satisfy the Jacobi identity but we will argue that one of them is more convenient than 
the other. They will be denoted by {^-jaankj and {■rYna.nkj where j denotes the rank of the matrix 
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A. Both brackets define the equations of motion ()2.1ip in the sense that the directional derivative of 
any function / = /(7,K) G C°°{R^ x M^) along the flow is given by / = {/, T^TjIaank, = {/, ^T^ILnk,- 
The geometric interpretation of these brackets is the subject of the subsequent sections. Concretely, 
in section [5] (Theorems 15.31 and 15. 4p . we will show that the bracket {-j-jRankj arises as the reduction 
of the nonholonomic bracket introduced in [M], and that {•, -Ir^^j^j arises as the reduction of a gauge 
transformation of the nonholonomic bracket. 

The following definitions will be useful in our discussion of the utility of the two brackets: 
Definition 2. Let P be a manifold equipped with an almost Poisson bracket {•, •}. 

1. The (almost) Hamiltonian vector field Xj of a function / € C°°{P) is the vector field on P 
defined as the usual derivation Xf{g) = {g, /} for all g € C°°{P). 

2. The characteristic distribution of {•,•} is the distribution on the manifold P whose fibers are 
spanned by the (almost) Hamiltonian vector fields. 

3. Due to Leibniz condition of {•,•}, there is a bivector field vr S T{/\^{TP)) such that for f,g & 
C°°{P) we have Tr{df,dg) = {f,g}. We say that vr is the bivector field associated to {•,•} and 
we denote by vr^ : T*P — > TP the map such that /3(7r^(a)) = 7r(a,/3). We will occasionally refer 
to bivector fields simply as bivectors. Note that the characteristic distribution is the image of 
TT^ and the Hamiltonian vector field Xf = — 7r^((i/). The 3-vector field [vr, vr], where [•,•] is the 
Schouten bracket, may be different from zero, and it measures the failure of the Jacobi identity 
through the relation 

{/, {g, h}} + {g, {h, /}} + {h, {/, g}} = ^i[,,,](4f, dg, dh), (2.12) 
for f,g,hGC°-iP). 

4. The bracket is called conformally Poisson if there exists a strictly positive function ip G C°°{P) 
such that the bracket ip{-, •} satisfies the Jacobi identity, i.e. [ip7r,ipTr] = 0. 

The well known symplectic stratification theorem states that the characteristic distribution of a 
Poisson bracket is integrable and its leaves are symplectic manifolds. Since multiplication of an almost 
Poisson bracket by a positive function does not change the characteristic distribution, a necessary 
condition for an almost Poisson bracket to be conformally Poisson is that its characteristic distribution 
be integrable. 

We now come back to the discussion of our example for the different values of the rank of A. 

If A has rank 3. Li this case A-^ = A^ and K = (I + mr E)fl where E denotes the 3x3 identity 
matrix. It follows form (|2.1ip that the rotational motion of the body is the same as that of a free rigid 
body whose total inertia tensor is I + mr'^E. It is trivial to write $7 = (I + mr^i?)^"'^K and it is clear 
that equations ()2.1ip define a closed system in x M^. 
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The two brackets for the system for functions f,g €z C°°{M x M ) are given by 

rl-y I ' 

(2.13) 



{ /, 9h...s = -(K + mr^fi) ■ i ^ ^ ^) - 1 ■ ( ^ x - ^ ^ 



The above brackets are quite different. On the one hand, the bracket {•, - jHanka satisfies the Jacobi 
identity. It in fact coincides with the Lie-Poisson bracket on the dual Lie algebra se(3)*. On the other 
hand, the bracket {•, •jaanks is not even conformally Poisson as the following proposition shows. 

Proposition 2.2. The characteristic distribution of the almost Poisson bracket {•, - jRanks defined in 
(j2.13|) is not integrable. 

Proof. The (almost) Hamiltonian vector field Xj of a function / G C°°(1R^ x M^) corresponding to the 
bracket {■, - jRanks is given by 



2r^^ df df\ d f df\ d 

X,= ^(K + n^r^l2)x^ + 7X^j.^+(7X^j.^ 

and it is annihilated by the non-closed one-form 

X = 7 • dK + (K + mr'^n) ■ d-y. 

We have 

2r, \ ^ /.^ , 2n ^ 9 d d _ d d 



dK2 ' 'dKs '■'972 ^73' 8X2 ' dK^' 

and thus ^ ^ 

x{[X^,,Xk,]) = -dxiX^^XK,) = -mr^ (t^^ + , 2^ 2 ) + 0- 

This shows that the commutator \X^^ , X^i ] does not belong to the characteristic distribution which is 
therefore not integrable. □ 

Therefore, to obtain a true Hamiltonian formulation of the reduced equations of motion in the case 
where the rank of A is 3, one needs to work with the bracket {•, -jRa^ni^g- 



If A has rank 2. As mentioned before, this case has the Chaplygin sphere as a particular example. 
The analysis of the two brackets has been done in ^2i|. We include it here for completeness and to link 
it with clarity to other results of the present work. 

In view of the form of A given in item (ii) of Definition [H we can write A^A = E — esel' and thus, 
according to (j2.9p . we get 

K = (I + mr'^E)il - mr'^{n ■ 7)7, 

which is precisely the expression for the angular momentum about the contact point for the Chaplygin 
sphere. 
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The angular velocity Q can be written in terms of K and 7 as 

n = {I + mr'E)-'K + mr' ( „^ ^ ' ^^''j^^ ^ ) (I + mr^Ey'^, 

^ ' V||7||2 - mr27 • (1 + mr2£;)-i7y ^ ' 

so both the equations ()2.1ip and the Hamiltonian H-ji are well defined on x M^. 

In this case, the two brackets for the system for functions f,g£ C°^{W^ x M^) are given by 

^ -(K - ,n^in . . . ^) - , . X I - ^ . |i) . 

None of the above brackets satisfies the Jacobi identity but it is preferable to consider {•, ■Y^^.^i^i- The 
reason is that this bracket is conformally Poisson with conformal factor 

= y/\ \j\ |2 - mr2 (7 • (I + mr'^E)-^-/). (2.15) 

This important observation was first made in The characteristic distribution of {•,-}n!^nk2 thus 
integrable. The generic leaves are the level sets of the Casimir functions Ci(K,7) = K • 7 and 
C2(7) = IItIP- Another important feature of this bracket is that it is twisted Poisson (in the sense of 
[261 [33]) as will be shown in Section [5^31 (Theorem 15. 7p . 
On the other hand, similar to Proposition 12.21 we have 

Proposition 2.3 (|21|). The characteristic distribution of the almost Poisson bracket {•, •}Rank2 defined 
in ()2.14p is not integrable. 

This can be shown exactly as we did for Proposition 12.21 Therefore if the rank of A is 2, just as in 
the case of rank 3, a Hamiltonian formulation of the reduced equations can only be obtained if we work 
with the bracket {•,-}Rank2- However, in this case one needs to multiply the bracket by a conformal 
factor. This can be interpreted as a time reparametrization, see the discussion in Section 14. 4i 



If A has rank 1. Taking into account the form of A given in item (Hi) of Definition [T] we have 
A^A = 6363^ and thus, in view of (|2.9p . we get 

K = in + mr^(n -7)7. 

The expression for the angular velocity Ct in terms of K and 7 is 

n = r'K - mr^ f^^^iLpL^—) rS, 
V I ItI r + "T-r^ (7 • 11 7)/ 

so again, both the equations (|2.1ip and the Hamiltonian T-L-ji are well defined on x M'^. 
This time, the two brackets for the system are given by 

. -,K . ..Aa . . . ^) - , . . I - ^ . I) , 
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for functions f,g £ C°°{R^ x R^). 

The properties of the brackets above are very similar to those obtained in the case where the rank 
of ^ is 2 except that the roles of {•, -iRanki and {•, -j^anki reversed. 

This time one can show that 

Proposition 2.4. The bracket {•, - jRanki defined in ()2.16p is conformally Poisson with conformal factor 

ip{-y) = ■\/||7|p + mr^ (7 ■ I-I7). 

Proof. We have to prove that the scaled bracket on TZ defined as ip{-, -jRanki satisfies the Jacobi identity, 
i.e., 

Rankl ) Rankl + ^W{f2,f3} Rankl ) /l} Rankl + ^W{f3,fl} Rankl ) /2}Rankl = 

for all /i, /2, /s € C°°(M^x]R^). In view of the derivation properties of the bracket, is enough to show the 
identity for the coordinate functions Ki,ji. In our case, since {7i, 7j}Ranki = 0, it is immediate to check 
that the identity holds if two of the three functions are 7j's. A long but straightforward computation 
shows that the identity holds for the following three choices of functions /i = i^i,/2 = K2, fs = 71; 
/i = Ki,f2 = K2, fs = 73 and fi = Ki,f2 = K2, fs = K^. Since the definition of the bracket is 
symmetric with respect to the coordinate functions Ki,'ji, and since the Jacobi identity trivially holds 
if two of the three functions /i, /2, /a are equal, all of the other cases are either trivial or analogous. □ 

Hence, the characteristic distribution of {•, -jRanki is integrable and the generic leaves are again the 
level sets of the Casimir functions Ci(K,7) = K-7 and 6*2(7) = IItIP- It '^ill also be shown in Section 
15.31 that {^-jRanki is twisted Poisson. 

On the other hand, analogous to Propositions 12.2 1 and 12.3] we have 

Proposition 2.5. The characteristic distribution of the almost Poisson bracket {•, - jRa^ni,! defined in 
()2.16p is not integrable. 

The proof is again similar. 

Thus, this time the Hamiltonian structure of the reduced equations can only be obtained with the 
bracket {•, -iRanki) again through the multiplication by a conformal factor that is interpreted as a time 
reparametrization. 



If A has rank 0. In this case A is the zero matrix and the constraints are holonomic and can be 
seen as a conservation law for the standard free rigid body. We have K = Ifl and clearly the equations 
()2.1ip and the Hamiltonian Ti-ji are well defined on x M^. The two brackets are given by 

f f X _ -w f '^f ( ^9 ^9 

{/, 5)RankO = -^-^^^9|^j-^-^^Xa:^-a|^X^^ 

/f V (T^ 2o^ fdf dg dg df^ ^^'^^^ 

{/,5}Ranko = -(K-mrn)- — X— -7-(^x^-^x^ 



5K 9Ky ' ^7 5K ^7 

The situation is analogous to that of the case when the rank of ^ is 3 but, once more, the roles of the 
brackets are reversed. While {■, -jRanko coincides with the Lie- Poisson bracket in the dual Lie algebra 
se(3)* (and hence satisfies the Jacobi identity), we have 
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Proposition 2.6. The characteristic distribution of the almost Poisson bracket {•, •}^j,„i,o defined in 
(j2.17p is not integrable. 

The proof is identical to that of Proposition 12.21 

So in this case, the Hamiltonian structure of the reduced equations ()2.1ip can only be seen by 
working with the bracket {•, -jRanko- 

2.4 Symmetries 

The reduced equations ()2.1ip can be interpreted as the output of a reduction process that we now 
explain. We begin by noticing that the configuration space Q = SO (3) x can be endowed with the 
Lie group structure of the three dimensional euclidean transformations SE(3). The group multiplication 
is given by 

(5l,Xi)(5r2,X2) = (5152,51X2 +Xi). 

Let H be the Lie subgroup of SE(3) defined by 

H = {{h,y) e SE{3) : /leg = 63}. 

For matrices A satisfying any of the conditions of Definition [H it follows that HA = Ah whenever 
{h, y) € H. We consider the left action of H on Q by left multiplication. The tangent lift of the action 
to TQ maps 

(^,y) : (5,x,C4J,x) {hg,hyL + y,hu,h±) or (/i,y) : (5,x,f2,x) ^ /ix + y, fJ, /ix), 

depending on the trivialization of S0(3) that one is working with. Notice that the Lagrangian C given 
by ()2.6p is invariant under the lifted action. Moreover, since h commutes with A for any {h, y) E 
the constraint (|2.3p is also invariant. 

The momenta (K, p) are geometrically interpreted as coordinates on the fibers of the (trivial) 
cotangent bundle T*Q. The constraint space M <ZT*Q \s determined by the condition p = mrAgfl, 
so the triple (5,x, K) G S0(3) x x specifies a unique point in Ai. Reciprocally, any point in A4 
can be represented by a triple (5,x, K). 

By invariance of the Lagrangian and the constraints, the lifted action of H to T*Q leaves A4 
invariant and, therefore, restricts to A4. The restricted action is free and proper so the orbit space 
TZ := M./H is a smooth manifold. The reduced space TZ can be identified with S^ xM^; the projection 
p : 7W — )• 7?. is given by 

/5(5,x,K) = (7,K), (2.18) 

where 7 = ^"^63 € S^, and is a surjective submersion. The conditions he^ = 63 and hA = Ah, that 
are satisfied for (/i, y) € ensure that the above mapping is well defined (in particular notice that K 
is invariant). The reduced equations on TZ are precisely (|2.1ip when restricted to the level set ||7|| = 1. 
Notice that TZ inherits the (trivial) vector bundle structure S'^ xM^ — > S^ from Ad. 

In this sense, the entries of 7 should be considered as redundant coordinates for the sphere S^ and 
the entries of K as coordinates on the fibers of TZ. Notice that, for any j = 0, . . . , 3, both brackets 
•}Rankj and {■, •jaj.nkj Testrlct to the level set ||7|| = 1 since (^2(7) = ||7|P is a Casimir function. 
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2.5 Kinematics and integrability of the constraint distribution 

The constraint distribution on Q defined by equation (|2.3|) has fundamentally diff'erent properties 
according to the rank of the matrix A satisfying the conditions of Definition [TJ On one extreme we 
have the case where ^ = and the distribution is integrable (the 3-dimensionaI integral leaves are 
given by S0(3) x {x} for x € M^). As mentioned before, in this case the constraints are holonomic and 
the problem reduces to the classical free rigid body problem (the center of mass of the body x remains 
constant in our inertial frame). 

The extreme opposite case occurs when rank A = 3. In this case the corresponding distribution 
is completely nonholonomic or bracket- generating, see e.g. [31]. By Chow's theorem, any two points 
in the configuration space Q can be joined by a curve {g{t),x{t)) satisfying the constraints. Thus, at 
least at the kinematical level, there are no restrictions on the values of x. 

The cases where the rank of ^ is 1 or 2 lie in between the situations described above. If the rank 
of yl = 2, the third component X3 of x remains constant during the motion. This is in agreement with 
our observation that the Chaplygin sphere problem is a particular case of this type of constraints - the 
sphere rolls on a horizontal plane 2:3 = const. In this case, the constraint distribution is non-integrable 
but is nevertheless tangent to the foliation of Q by 5-dimensional leaves defined by constant values of 
X3. 

Finally, for the case where the rank of A equals 1, the first two components xi, X2, of x remain 
constant during the motion. The body goes up or down along the X3 axis at a speed that is proportional 
to its angular velocity about this axis. This time the constraint distribution is non-integrable but 
tangent to the 4-dimensional leaves given by constant values of xi and X2- 

Without going into technical definitions, we simply state that the degree of non-integrability of the 
constraint distribution increases with the rank of A, passing from an integrable distribution if ^ = 
to a completely nonholonomic distribution if rank ^4 = 3. It is interesting to see how this correlates 
with the need of a gauge-transformation to Hamiltonize the problem (Remark 15. 6p . 

3 Regular almost Dirac structures 

This section is concerned with the geometry of (almost) Dirac structures [13| . with focus on the regular 
case, as well as their gauge transformations [331. Our exposition begins with the basic definitions but 
it reaches out to present some original results. As we will see, these geometric structures provide the 
setup that gives rise to the different brackets introduced in Section [2l Although we will be mostly 
interested in the geometry of bivector fields, our discussion is presented at the general level of (almost) 
Dirac structures, as they provide the framework in which gauge transformations are most natural. We 
refer the reader who is unfamiliar with this material to the introductory exposition in section 2 of [9]. 

3.1 Dirac and almost Dirac structures 

A Dirac structure on a manifold P is a subbundle L of the Whitney sum TP T*P such that 



15 



(i) L is a maximal isotropic subbundle of TP © x*P with respect to the pairing (•, •) given by 
{{X,a),{Y,(3)) =aiY)+(3{X), for (X,a),(y,/3) eTP®T*P. 

{ii) r(L) is closed with respect to the Courant bracket defined on T{TP © T*P) given by 

[(X,a),(y,/3)1 =([X,y],£x/3-iyda), 

for Gr(rp©r*p), i.e., ir(L),r(L)i cr(L). 

The underlying manifold P is sometimes referred to as a Dirac manifold. 

Let pri : TP © T*P — ?> TP be the projection onto the first factor of TP © T*P. A Dirac structure 
L on the manifold P carries a Lie algebroid structure with anchor pri\i : L — )• TP and bracket 
given by the Courant bracket |-, •] restricted to r(L). It follows that, for a Dirac structure L, the 
distribution pri(L) C TP is integrable, i.e., P can be decomposed into leaves O such that for each 
X (z P, TxO = pri{Lx). If pri(L) has constant rank (i.e., pri{Lx) C T^P has the same dimension for 
all X € P), then we say that L is a regular Dirac structure, and pri{L) defines a regular foliation. Just 
as a Poisson manifold P is the disjoint union of its symplectic leaves, each leaf of a Dirac manifold P 
carries a presymplectic form. 

Examples Let be a closed 2-form and vr be a Poisson bivector field, and consider the maps Vi' : 
TP T*P given by n\X) = n{-,X) = -ixfl = for X G TP and vrtt : r*P ^ TP as in Definitions 
Then 

Ln := graph{n^) = {{X, a) eTP(B T*P : ix^ = -a} 

and 

:= graph{J) = {{X, a) £TP® T*P : 7r»(a) = X} 

are Dirac structures. Note that pri identifies Tq with TP as Lie algebroids. Similarly, can be 
naturally identified with T*P, and the Lie algebroid structure induced on T*P by has anchor 
J :T*P ^TP, and bracket 

[a, 13]^ = i^^tt(„)(/3) - i^^tt(/3)(a) - (i(7r(a,/3)) = £^tt(„)/3 - M(/3)C?a. 

This bracket is uniquely characterized by [df^dg],^ = d{f,g} and the Leibniz identity. 

If a subbundle L of TP®T*P satisfies (i) above, but not necessarily (ii), then L is called an almost 
Dirac structure. Condition (ii) is called the integrability condition. We say that L is a regular almost 
Dirac structure when the distribution pri{L) on P has constant rank. Notice that this distribution 
might not be integrable in the almost Dirac case. 

Examples If is an arbitrary 2-form or vr an arbitrary bivector field, then their graphs Lq and 
Ltt are almost Dirac structures. The failure of the integrability with respect to the Courant bracket 
of Lq and L^^ is measured by dQ and |[7r,7r], respectively. For L^r = graph{7r'^), the distribution 
7r''(T*P) = pri(L^) is generally non-integrable. If it has constant rank we call the almost Poisson 
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structure regular. Note that the bracket [-j-j-K defined as in ()3.ip is M— bilinear, skew-symmetric and 
satisfies the Leibniz identity. However, in general, vr" does not necessarily preserve the bracket; instead, 
(see e.g. [8]), 

7r«([a,/3]0 = [7r«(Q),7r«(/3)]-ii,^^[7r,7r], foi a, (3 e n^P). (3.19) 

Note that an almost Dirac structure L on P is of the form = graph{'K^) for a bivector vr if and 
only if 

TPnL = {0}, (3.20) 

and L is of the form = graph{^) for a 2-form Q. if and only if T*P f] L = {0}, see |13) . Another 
example of an almost Dirac structure that will be very useful for our purposes is given by L C TP®T*P 
defined as 

L := {{X, a) eTP® T*P : X e F, ixfilp = -oIf}, (3.21) 

where F C TP is a subbundle, is a 2-form on P and • \p denotes the point-wise restriction to F. If 
the subbundle F is an integrable distribution and is closed, then L is a Dirac structure. 

Proposition 3.1. The following statements hold: 

(i) There is a one-to-one correspondence between regular almost Dirac structures L C TP®T*P and 
pairs {F,^lp), where F is a regular distribution on P and G r(/\^F*). 

(a) Let F C TP be a regular distribution on P. Given a section Hp r(/\^ F*), there exists a 2-form 
Q on P such that Q\p = Qp. 

Proof (i) Let L C TP ®T*P be a regular almost Dirac structure with distribution F := pri{L) C TP 
on P (not necessarily integrable). Consider the section Vtp in T{/\^ F*) given, at each x G P, by 

Q.p{x){X(^^),Y(^^)) = -a(2.)(y(^)), for X,y G T{F) such that (X(^),a(^)) G L^;- 

It is a straightforward computation to see that Oi;' is well defined, i.e., it is independent of the choice of 
a. Conversely, given a regular distribution F on P and dp G T{/\'^ F*), we may define the subbundle 
L C TP © T*P as the pairs {X, a) such that X F and ix^F = —alp- 

(ii) Let W C TP be a regular smooth distribution such that it is a complement of F on P, i.e., 
T^P = Fx (B Wx for each x G P (e.g., can be chosen to be the orthogonal complement of Fx with 
respect to a Riemmanian metric). The 2-form on P can be defined by 

a{X,Y) = ^lp{Xp,Yp), 

for X, y G X(P) such that X = XiT' + Xw and y = Yf + Fw, where X^, G r(P) and Xw, G 
r(VV). Differentiability of follows from its definition and the smoothness of F and W. □ 

Corollary 3.2. Given a regular almost Dirac structure L, there exists a 2-form Q on P and a regular 
distribution F C TP such that L is written in the form ()3.2ip . 

Remark 3.3. Note that the 2-form Q G J7^(P) is not uniquely defined and, in general, there is no 
canonical choice for it. o 
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Given a subbundle F C TP, we say that a section in r(/\^ F*) is nondegenerate if it is nonde- 
generate as a bilinear form on F at each point. It fohows from (j3.2ip that Keic{il.F) = Lf] TP, and as 
a consequence of condition (|3.20p we obtain 

Corollary 3.4. Let L be a regular almost Dirac structure and (FjClp) the pair associated to it in the 
sense of Proposition COl Then ilp is nondegenerate if and only if L is the graph of a bivector field vr. 
Explicitly, the relation between {F,Qp) and vr is 

7r^(a) = —X if and only if ix^F = Oi\F, 

where X G r(F) and a G n^{P). 

Following notation of Definition [21 if {•,•} is the bracket associated to the bivector field vr in the 
above Corollary, then {f,g} = ^F{Xf,Xg), for all f,ge C°°{P). 

3.2 Twisted Poisson and twisted Dirac structures 

Poisson structures may be viewed as encoding integrability in two levels: first, the characteristic dis- 
tribution 7r''(T*P) C TP is integrable, i.e., tangent to leaves; second each leaf carries a nondegenerate 
2- form that is closed (and this leads to the Jacobi identity). Twisted Poisson structures are special 
types of almost Poisson structures that retain the integrability of ■k'^[T*P) but allow the leafwise 2- form 
to be non closed. These objects turn out to be related to nonholonomic systems. We start with the 
more general notion of twisted Dirac structures. 

Consider a closed 3- form (f) on P, and define the <j)-twisted Courant bracket [33] as follows: 

I(X,Q),(y,/3)l^ = {[X,Yi£xp-Wda + {x^Y(|>), (3.22) 

for (X, a) and iY,j3) in r(TP©T*P). Now, a subbundle L of TP(BT*P is a (p-twisted Dirac structure 
|33| if L is maximal isotropic with respect to (•, •) and the integrability condition 

[r(L),r(L)i^cr(L) 

is satisfied. 

As in the ordinary case, a twisted Dirac structure L on P induces a Lie algebroid on L given by 
the anchor map pri\L and the bracket |-, ■](/)|r(L)- Therefore pri{L) is an integrable distribution on P. 

Examples If Q. is any 2-form on P, then L^j = graph{Q^) is ((ir2)-twisted Dirac. One may check that 
a bivector field vr on P such that L^^ = graph(Tr'^) is twisted Dirac verifies (see |33j ) 

i[7r,7r] =7r« ((/<). (3.23) 

The following result gives more examples: 

Theorem 3.5. Let L C TP © X*P be a regular almost Dirac structure such that pri{L) C TP is an 
integrable distribution on P. Then, there exists an exact 3-form (f) with respect to which L is a (/)-twisted 
Dirac structure. 
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Proof. Let F := pri{L) C TP and 0^ be the section in r(/\^ F*) associated to L given by Proposition 
I3.1i Since F is integrable, VLp defines a 2-form on each leaf O where Fx = T^O at each x € P. By 
Corohary 13.21 there exists a 2-form Q on P such that i*Q^ = Q.o where lq '■ O ^ P \s the inclusion. 
We assert that L is a ((ir2)-twisted Dirac structure. In fact, for (X, a) and iY,fi) in r(L), 

I(X,a),(y,/3)](rfO) ={[X,Yl£x^i-Wda + \xAYdn)^T{L) 

if and only if 

i[x,y]^|F = -{£xP - iyda + ixAYdO,)\F. 
Since F is an integrable distribution we obtain that, 

-{£xf3 -iYda + ixAYdCl)\F = -£x{(3\f) + iYd{a\F) - ixAydi^lp) 

= £x^Y^F ~ lYdix^F ~ ^XAYd^F — \x,Y]^F ■ 

which completes the proof. □ 
Remark 3.6. 

(z) Note that if L is a 0— twisted Dirac structure then L is also twisted with respect to any closed 
3-form (p' such that {(j) — (j)') vanishes on the leaves. 

{ii) There is no canonical choice for the 3-form (p given in Theorem 13.51 

o 



Tvi^isted Poisson bivectors. 

Bivector fields tt such that = graph{7r'^) is a 0-twisted Dirac structure are called 4>-twisted Poisson 
bivectors [MIES], i.e., tt verifies condition ()3.23p . We are especially interested in these kind of bivector 
fields since, as we will see, they appear naturally in the examples of nonholonomic systems introduced 
in Section [2j If {•,•} is the bracket given by the i;^)- twisted Poisson structure vr, then relation ()2.12p 
becomes 

{/, {g, h}} + {g, {h, /}} + {h, {/, g}} + cl){Xf,Xg,Xh) = 0, 

for f,g,h E C°°{P) and Xf = {•,/}. So the failure of the Jacobi identity is controlled by the closed 
3-form (f). 

Example If vrtt : r*P ^ TP is an isomorphism, then vr is (/)-twisted and <f> = where Cl is the 
nondegenerate 2-form associated with tt, i.e., fi'' o tt^ = Id where as usual il!'{X) = —ix^- Other, less 
trivial examples, will be presented ahead in Corollary 13.71 

Let vr be a (/>-twisted Poisson structure on the manifold P and [•, be the bracket on T*P given 
by ()3.ip . Note that vr^ does not preserve this bracket. However, using ()3.19p and ()3.23p we obtain 
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for 1-forms q, /3 on P. The i;^-twisted Courant bracket induces a modification of the bracket ()3.ip via 
the identification of T*P and -L^r, 

[a,/3]0 = £^tt(Q,)/? - i^tt(/3)(ia + i^tt (a) attB 

such that {T*P, [■, -l^j^r") is a Lie algebroid [33] (see also [8]). The characteristic distribution 7r^(T*P) 
defines an integrable distribution on P (that may be singular). Each leaf O of the corresponding 
foliation of P is endowed with a non-degenerate 2-form that is not necessarily closed. If vr is 
twisted, then dOc) = '-c)';^, where to : O ^ P is the inclusion. 

Important examples of twisted Poisson structures are contained in the following Corollary of The- 
orem 13.51 

Corollary 3.7. Let tt be a bivector field on P with an integrable regular characteristic distribution. 
Then, there exists an exact 3-form (p on P with respect to which vr is (p-twisted. 

Remark 3.8. Mechanical Example. It is shown in [2D| that the (semi-direct) product reduction of 
the Veselova system yields a regular conformally Poisson bracket on the reduced space. It follows 
that its characteristic distribution is integrable and thus, by Corollary 13.71 it is also twisted-Poisson. 
This is a first example of a nonholonomic system whose reduced equations are formulated in terms 
of a twisted-Poisson bracket. Other examples (related to the motion of a rigid body with generalized 
rolling constraints) are made explicit in Section 15.31 o 

Remark 3.9. An interesting question, that remains to be answered, is to give a characterization of 
almost Poisson brackets possessing an integrable characteristic distribution that is non-regular. o 



Regular conformally Poisson bivectors. 

An interesting class of almost Poisson structures admitting an integrable characteristic distribution is 
given by conformally Poisson structures. Recall from Section [2] that they are bivector fields vr for which 
exists a strictly positive function (p € C°°{P), such that ipir is Poisson. A conformally Poisson manifold 
(P, vr) is the disjoint union of conformally symplectic leaves. 

Note that this property is stronger than asking for [P, vr) to be a Jacobi manifold since a conformal 
factor implies the global existence of a function such that (/jtt is Poisson, while in Jacobi manifolds the 
factor (p may be only locally defined. 

From Corollary 13.71 we observe that any regular bivector admitting a conformal factor is also a 
twisted Poisson bivector. The following proposition explains the relation between these two properties. 

Proposition 3.10. Let tt be a regular conformally Poisson bivector field on P with conformal factor 
If. Let € VL^{P) be as in Corollarv \3.Si Then any closed 3-form verifying ()3.23p for vr coincides with 

—dip Ail. on the leaves. 

f 

Proof. Since vr admits a conformal factor ip E C°°{P), then [tTjTt] = ^X^p A vr. On the other hand, 
if 0, is the 2-form on P associated to vr given by Corollary [321 then for 51,(72 S C°°{P) we have 
Q{Xg-^,Xg^) = TT{dgi,dg2). Thus, for 51, 52, 53 S C~(P) we obtain 

}-[7r,Tr]{dgi,dg2,dg3) = -X^ ATT{dgi,dg2,dg3) = --d(p A n{Xg^, Xg^, Xg.^). 
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Then we conclude that any closed 3-form (j) satisfying ()3.23p coincides with {^d^p A il) on the leaves. 

□ 

Let (P, vr) be a regular conformally Poisson manifold. A 2-form on P satisfying the conditions 
of Corollary 13.21 verifies that is conformally symplectic on each leaf O. However, Vt may not 
necessarily be conformally closed. 

3.3 Gauge transformations. 

In this section we will consider a natural action of the abelian group of 2-forms on P on the almost 
Dirac structures on P. This action is given by gauge transformations of almost Dirac structures by 
2-forms, and it was introduced in [33j. 

More precisely, consider an almost Dirac structure L in TP (BT* P . A gauge transformation by the 
2-form B is a map tb : L ^ TP ® T*P, given by Tij((X,a)) = {X,a + ixB) for (X, a) € L. The 
subbundle tb{L) of TP T*P given by 

TBiL) = {{X,a + ixB) : (X, a) G L} 

is an almost Dirac structure. If the 2-form B is closed and L is Dirac then tb{L) is again Dirac. Thus, 
the 3-form dB is what determines the integrability with respect to the Courant bracket. It is a direct 
computation to see that if L is a </)-twisted Dirac structure, then the gauge transformation of L by the 
2-form B \s {(f) — (ii?)-twisted Dirac (see e.g. [33]). 

If Li and L2 are almost Dirac structures on P and there exists a 2-form B on P such that tb{Li) = 
L2, then we say that Li and L2 are gauge equivalent or gauge related. 

Note that a gauge transformation does not modify the distribution pri{L). So, for a Dirac struc- 
ture L, the foliation associated to L will be the same as the one associated to tb{L). However, the 
presymplectic form on each leaf is modified by the pullback of B to the leaf. If L is a regular almost 
Dirac structure determined by the pair {F,Q,f) in the sense of item (i) of Proposition 13.11 a gauge 
transformation by the 2-form B corresponds to the operation: 

TB : {F,VLf) ^ [F.np - B\f). (3.24) 

Theorem 3.11. Any two regular almost Dirac structures Li and L2 are gauge related if and only if 
pri(Li) = pri{L2). 

Proof. It remains to prove the "only if" part of the statement. Let us denote F := pri{Li) = pri{L2) 
and let Qp and Qp be the 2-sections associated to Li and L2 respectively (Proposition 13. II (i)). Define 
the section Bp G r{/\'^{F*)) by Bp := n], - and let B G n'^{P) such that B\f = Bp (Prop. O 
(ii)). We claim that tbLi = L2. In fact, if {F,Qp) is the pair associated to the almost Dirac structure 
T^Li, then by equation ()3.24p . 

i^p = — Bp = ^"p. 

Since the sections associated to tbLi and L2 coincide, by Proposition 13. II (i) we conclude that r^Li = 
L2 which means that Li and L2 are gauge related. 

□ 
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We are especially interested in gauge transformations of almost Poisson structures. Consider the 
almost Poisson manifold {P, vr) and a 2-form B on P. Then, the gauge transformation of := 
graph{7r'^) is tb{Lt^) = {{X,a + ixB) € TP T*P : X = T:\a)} which does not necessarily 
correspond to the graph of a new bivector vr^. A necessary and sufficient condition for this to happen 
is that 

TB(L,)nrp = {0}, 

which is equivalent to the fact that the endomorphism Id + i?^ o tt" : T*P — ?> T*P is invertible [33] . 
Indeed, if such a bivector field vr^ exists, then, in view of ()3.3p . for any 1-form q on P we have 

TB ((vr''(a) , a)) = {tt^o) , a + = ((vr^)" (^a + , a + • 

Thus, vr^ is characterized by the condition 

(7r^)«(a + M(,)i?) =7rtt(a), (3.25) 

and we can write 

(^S)tt = 7r«o(Id-S^o^tt)-i. (3.26) 
In particular, if vr and tt^ are non-degenerate bivector fields, equation (j3.26p is equivalent to 

((vr^)tt)-i = (vr«)-i - 

Therefore, any two 2-forms on a manifold are gauge related. This is not necessarily the case with 
bivector fields. A necessary condition is that their characteristic distributions coincide. In view of 
Theorem 13.111 this condition is also sufficient if such distributions are regular. 

Although gauge related bivectors have the same characteristic distribution, their Schouten brackets 
[7r,7r] and [7r^,7r^] may not coincide. The following proposition makes this precise. 

Proposition 3.12 (|33j). // a (p-twisted Poisson bivector vr is gauge related with another bivector tt^ 
via the 2-form B, then vr^ is {(/} — dB)-twisted. That is, 

^[n^,n^] = (n^)i{cP-dB). 



4 Applications to Nonholonomic Systems and Hamiltonization 

In what follows, we will analyze the geometry of nonholonomic systems in the framework presented in 
the previous section. We introduce gauge transformations of the bracket describing the nonholonomic 
dynamics in order to study the process of Hamiltonization. 
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4.1 Nonholonomic systems 



A nonholonomic system consists of an n-dimensional configuration manifold Q with local coordinates 
q € [/ C M", a Lagrangian C : TQ — t- M of the form £(q, q) = ^^(q)(q, q) — ^(q), where ^ is a kinetic 
energy metric on Q and F : Q ^ M is a potential, and a regular non-integrable distribution V C TQ 
that describes the kinematic nonholonomic constraints. In coordinates, the distribution T> is defined 
by the equation 

e(q)q = 0, (4.27) 

where e(q) is a /c x n matrix of constant rank k where k < n\s the number of constraints. The entries 
of e(q) are the components of the M'^-valued constraint 1-form on Q, e := e(q) dq. 

The dynamics of the system are governed by the Lagrange-D'Alembert principle. This principle 
states that the forces of constraint annihilate any virtual displacement, so they perform no work during 
the motion. The equations of motion take the form 

( — 1 = M^e(q). (4.28) 



dt \di\) (?q 

Here /x : TQ — )• M'^ is an M'^-valued function whose entries are referred to as Lagrange multipliers. 
Under our assumptions, it is uniquely determined by the condition that the constraints (j4.27p are 
satisfied. The equations (j4.28p together with the constraints (j4.27p define a vector field on T> whose 
integral curves describe the motion of the nonholonomic system. A short calculation shows that along 
the flow of y^, the energy function Ec := ^ • q — is conserved. 

The above equations of motion can be written as a first order system on the cotangent bundle T*Q 
via the standard Legendre transform, Leg : TQ — >■ T*Q, that defines canonical coordinates (q, p) on 
T*Q by the rule Leg : (q, q) i— > (q, p = dC/dq). The Legendre transform is a global diffeomorphism 
by our assumption that Q is a metric. 

The Hamiltonian function, Ti : T*Q M, is defined in the usual way n ■= Eco Leg~^ The 
equations of motion ()4.28p are shown to be equivalent to 

and the constraint equations ()4.27p become 

6(q)^ = 0. (4.30) 

The above equation defines the constraint suhmanifold A4 = Leg(P) C T*Q. Since the Legendre 
transform is linear on the fibers, is a vector sub-bundle of T*Q that for each q & Q specifies an 
n — k vector subspace of T*Q. 

Equations ()4.29p together with ()4.30p define the vector field on A4, that describes the motion 
of our nonholonomic system in the Hamiltonian side and is the push forward of the vector field by 
the Legendre transform. The vector field X^i^ is defined uniquely in an intrinsic way by the equation 

L*nQ = i^^dn^ n^T*e), (4.31) 
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where ilq is the canonical symplectic form on T*Q, l : Ai ^ T*Q is the inclusion and t : T*Q Q is 
the canonical projection. The constraints ()4.30p and their derivatives are intrinsically written as 

X„i, G C := TM n T, (4.32) 

where T is the distribution on T*Q defined as T{T*Q) : {t*€,v) = 0}. Denote by the 

pull-back of flQ to A4, i.e. := i*^Q. The following proposition is of great importance for our setup 
of the equations of motion as an almost Hamiltonian system. 

Proposition 4.1 (\38\ [2]). The distribution C on Ad defined by ()4.32p is regular, non-integrable, and 
the point-wise restriction of to C, denoted by is non- degenerate. 

The non-integrability of C is a direct consequence of the non-integrability of V. One shows that the 
rank of C is 2(n — k) and that along Ai we have the symplectic decomposition 

Tm{T*Q) = C(BC^'^, (4.33) 

where C^*? denotes the symplectic orthogonal complement to C. 

Since r*e vanishes on C, by restricting (j4.3ip to C and denoting Tij^ := i*^. G C"^{Ai), the 
equations of motion can be written in the appealing format 

ix ,^c = {dnM)c, (4.34) 

where {d'HM)c is the point- wise restriction of dV-M to C. The above equation uniquely defines the vector 
field and is central in our treatment; with this in mind, we collect the data of the nonholonomic 
system in the triple {Ai , flc , T-Lm ) • 

Even though ()4.34p defines the vector field X^^ uniquely, and resembles a classical Hamiltonian 
system, notice that since the distribution C is non-integrable, then is a section in /\^ C* ^ Ai (not 
a 2-form). 

Let (A^,r2c)^x) be a nonholonomic system. For every / G C°°(A^), let Xf denote the unique 
vector field on Ai with values in C defined by the equation 

ix^^c = {df)c, (4.35) 

where (d/)c denotes the point-wise restriction of df to C. The vector field Xf defined by equation 
()4.35p is called the the (almost) Hamiltonian vector field associated to f. 

Since Qc is nondegenerate, by Corollary 13.41 there is a unique bivector field n^i-, on Ai associated 
to the pair (C,Oc)) that is 7rfjj(a) = —X if and only if ix^c = ct|c. On exact forms we have, for 
/ G C^{M), 

iXf^c = {df)c if and only if 7rL(d/) = -X/, (4.36) 

which is consistent with notation of Definition [5J The bracket {•, -jnh on functions on Ai associated to 
the bivector vTnh describes the dynamics in the sense that 

XM){m) = XHj^{f){m) = {f,%MU{^) for ah / G C°^{M). 

It follows from ()4.36p that the characteristic distribution of the bracket {■, -jnh is C Since C is non- 
integrable then {., -jjiij is an almost Poisson bracket that does not satisfy the Jacobi identity. 
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Remark 4.2. If the constraint distribution V were integrable, the same would be true for the distri- 
bution C. Let TV C 7W be a leaf of the corresponding (regular) foliation of A4 (i.e. Cx = T^M for all 
X G N). In view of ()4.33p the submanifold M is symplectic. Therefore, in this case, our construction 
of {■, -jnh coincides with the usual construction of the Dirac bracket on each leaf N of the foliation of 
Ad (see Oils]). Hence, in this case, the Jacobi identity holds. o 

Inspired by equation ()3.2ip and our data, it is natural to define the almost Dirac structure on 
M by 

L^^:={{X,a) eTM®T*M : X e C, ix^m\c = -a\c}. (4.37) 
as considered (up to a minus sign) in [361 EZ] i and subsequently in [2^ . 

Proposition 4.3. Let 7r„ij he the bivector field defined in ()4.36p and {■, the corresponding bracket. 
The following statements hold: 

(i) The almost Dirac structures ^Tr^^ := graph{Tr^i^) and given in (j4.37p coincide. 

{a) The almost Poisson bracket {•, -jnh coincides with the classical almost Poisson bracket for non- 
holonomic systems defined in |34l [25] . 

Proof, (i) It is immediate since both almost Dirac structure are defined by the same pair (C, ^c)- 

(a) The bracket on A4 for nonholonomic systems introduced in [341 [25] was shown in [lU] to be 
given by 

{f,g} = nM{VYj,VY-g), for /,5€C°°(A^), (4.38) 

where V : T_m(T*Q) — t- C is the projector associated to the symplectic decomposition (|4.33p . and Yj 
the free Hamiltonian vector field on the symplectic manifold {T*Q, 0,q) defined by iy_r2Q = df, where 

/ is an arbitrary smooth extension of / to T*Q. 

It is easy to check that along A4 one has i-pyj^c = idf)c, so VYj' coincides with the almost 

Hamiltonian vector field Xj defined by equation (j4.35p . Therefore, for any f,g(z C°°{A4), 

{f,g} = nM{VYf-,rYg) = nc{Xf,Xg) = -df{nUdg)) = {f,gU. 

□ 

The second item of the above proposition should not be surprising since the expression (j4.38p is the 
nonholonomic version of the Dirac bracket (see discussion in |23[ 110]). Hence, its description naturally 
falls in the ambit of almost Dirac structures as described above. As a consequence of the above 
proposition, we can equivalently describe our nonholonomic system with the triple {M.,-k,^i,,'Hm)- 

Definition 3. Let (A1, vTnh, ^A4) be a nonholonomic system. 

1. The bivector field vTnh on A4 given by (|4.36p is called the nonholonomic bivector field and the 
bracket {•, -jnh is called the nonholonomic bracket. 

2. We say that an almost Dirac structure L describes the dynamics of the nonholonomic system if 
the pair {—X^i-^dl-LM) £ r(L). 
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4.2 Gauge transformations of the nonholonomic bracket 

The main idea of using gauge transformations in our setting is that it opens the possibihty to modify 
the geometric structure on A4 that describes the dynamics. 

Consider the nonholonomic system (7W , TTn^, ) and continue to denote L^^,^ = graph{Tr^^i^). The 
gauge transformation of Tr^h associated to a 2-form B on A4 gives 

TBiL^^J = {{X, a + ixB) eTM® T*M : 7rL(a) = X}. (4.39) 

First of all we are interested in knowing when the pair {—X^Y,,dl-iM) is a section of tb{Lt^^^^). On 
the other hand, we would also like to know whether the almost Dirac structure tb{Lt^^^) corresponds 
to the graph of a bivector field or not. 

If the 2-form B on M verifies ^ i3 = 0, then from equation ()4.39p we see that the pair 
(— Xnh, (IT-Lm) belongs to r(rB(L7rj^^)). Moreover, in view of ()3.24p . the gauge transformation of n^i^ by 
the 2-form B has the form 

TBiL^^J = {{X, a)£TM(S T*M : X G C, ixi^M - B)\c = -a\c}. (4.40) 
Thus the equations of motion (|4.34p are equivalently written as 

ix^^iVLc - Be) = {dnM)c, 

where Be is the point- wise restriction of to C. 

Therefore, as a particular case of Corollary 13.41 we observe that if the section — Be is non- 
degenerate then the gauge transformation of vrn^ associated to the 2-form B is again a bivector field 
TTj^ . It follows from (|3.26p that the non-degeneracy of f]_A4 — S on C is equivalent to the invertibility of 
the endomorphism (Id — B^ o vrSh) onT*JV[ . 

Definition 4. Let (P, tt) be an almost Poisson manifold with a distinguished Hamiltonian function 
H G C°°{P). Given a 2-form B on P, the gauge transformation of vr associated to the 2-form B is said 
to be a dynamical gauge transformation if 

(i) ixn^ — 0> where Xh is the (almost) Hamiltonian vector field associated to H and 

(ii) TB{graph(Tr^)) corresponds to the graph of a new bivector vr^, i.e., the endomorphism (Id— P^'ott") 
on T*P is invertible. 

Of course we are interested in dynamical gauge transformations of the nonholonomic bracket where 
the distinguished Hamiltonian function is T-Im and the corresponding (almost) Hamiltonian vector field 
is Xnh- 

Note that if vr is regular, by equation ()3.24p . the gauge transformation defined by B is determined 
by the restriction Bp of B to F where F := 7r^{T*P). Then condition (i) of the above definition is 
equivalent to ix^Bp = 0. 
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Remark 4.4. The definition of an affine almost Poisson bracket for a nonholonomic system made 
in [21] corresponds to a dynamical gauge transformation of the nonholonomic bracket by a 2-form 
B = —L*i}() where ^Iq is a semi-basic form 2-form on T*Q. The proof is analogous to that of item (ii) 
of Proposition 14.31 In this case, the hypothesis that flo is semi-basic implies that the condition (ii) of 
Definition H] is satisfied (see Proposition 14.61 below) . o 

After this discussion, we observe that it is more appropriate to describe a nonholonomic system 
by the triple {M,^,7iM) where ^ is the family of bivector fields that are related to vr^h through a 
dynamical gauge transformation. Notice that C is the characteristic distribution of any bivector field 
in 5- Thus, the (almost) Hamiltonian vector fields defined by the corresponding brackets satisfy the 
nonholonomic constraints. It follows that the bivector fields in the family ^ are almost Poisson bivectors 
in a "strong" sense since the non-integrability of C prevents them from being twisted or conformally 
Poisson. Our interest in considering this big family of brackets relies on the outcome after reduction. 
In the presence of symmetries, the bracket that Hamiltonizes the reduced equations may arise as the 
reduction of a member of ^ that is not necessarily tt^^,. 

Since the distribution C is regular we observe 

Corollary 4.5. [of Theorem l3.ll] . All bivectors with characteristic distribution equal to C are gauge 
related (in particular, gauge related to the nonholonomic bivector tt^i-,). 

We finish this section by discussing some cases for which the second condition in Definition 2] is 
satisfied. Recall that M C T*Q is a vector bundle over Q. We have 

Proposition 4.6. If B is a semi-basic 2-form on Ai, then the gauge transformation of tt^y, associated 
to B corresponds again to a bivector field. 

Proof. The graph of vrJi, is an almost Dirac structure corresponding to the pair (C, ilc) in the sense of 
Proposition 13.11 (see Proposition 14.31 {i)). Thus, in view of ()3.24p . the gauge transformation of TTn^ is 
the almost Dirac structure corresponding to the pair (C, {^m ~ B)\c). It is shown in [21) that if S a 
is a semi-basic 2-form, then the point-wise restriction of {^m — B) to C is non-degenerate. Thus, by 
Corollarv 13.41 tb{L.,^^^) corresponds to the graph of a bivector. □ 

In fact this proposition is a special case of the following result: 

Proposition 4.7. Let P Q be a vector bundle and vr a regular almost Poisson bivector on P. If 
for all semi-basic 1-forms a on P the vector field 7r^(a) is vertical, then the gauge transformation of it 
associated to a semi-basic 2-form B corresponds again to a bivector. 

Proof. Consider local bundle coordinates (q, p) € f/ x y C M" x M™, on P such that q are local 
coordinates on the base manifold Q. Since 7r^((iq) is vertical and B is semi-basic we obtain 

B^oTT^dq) = 
B'on^idp) = 6(q,p)dq, 
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where 6(q, p) denotes the m x n matrix with entries baj{(i,p) = {B o ir^i^dpa), -^) , for a = l,...,m, 
and j = 1, n. Thus, the matrix representation of the endomorphism (Id — o vr^) on T*P is 

/ Idnxn -fe(q,p)^ \ 

xm J 

This matrix has full rank and hence (Id — B^ o vr") : T*P — )• T*P is invertible. □ 

Remark 4.8. It is interesting for future work to drop the condition (ii) in Definition H] that requires 
TB{graph(j:\^)) to define a bivector field. In this case, the family ^ consists of all the almost Dirac 
structures that are gauge related to L^^ and that describe the nonholonomic dynamics. In this sense, the 
Hamiltonization of the problem is achieved if the reduction of a member of 5^ is a Dirac structure (but 
not necessarily a Poisson structure). This approach requires the consideration of a general reduction 
scheme for almost Dirac structures. However, we are unaware of any examples of nonholonomic systems 
that justify the need of such a general framework. o 



4.3 Reduction by a group of symmetries 

We now add symmetries to the problem and perform the reduction. Our interest from the point of 
view of Hamiltonization is to find a bivector field in the family 5 whose reduction is either Poisson or 
conformally Poisson (see Section [4. 4p . 

Let G be a Lie group acting freely and properly on Q. We say that that G is a symmetry of the 
nonholonomic system if the lifted action to TQ is free and proper, and leaves the constraint distribution 
T> C TQ and the Lagrangian C : TQ M invariant. 

Denote by ^ : G x T*Q T*Q the cotangent hft of the action to T*Q. If G is a symmetry for 
our nonholonomic system, then ^ leaves both the constraint submanifold A4 and the Hamiltonian 
H : T*Q M invariant. We continue to denote by ^ the restricted action to A4. 

One can show that the tangent lift of ^ to TA4, preserves the distribution C and the section 17^. 
As a consequence, if G is a symmetry group for our nonholonomic system, then the action ^ preserves 
the standard nonholonomic bracket {•, - jni,- That is, for /, 5 € C°°{M), we have 

By freeness and properness of the action, the reduced space TZ := M/G is a smooth manifold and the 
orbit projection map p : A4 — t- 7^ is a surjective submersion. Notice that IZ inherits a vector bundle 
structure from M over the shape space Q/G. Moreover, TZ is equipped with the reduced nonholonomic 
bracket {•, ■}^^^ that is characterized by 

{/, 5}rcd o p{m) ■={fop,go /j}„h(m) for m€ M and f,g€ G°°{n). (4.41) 

The corresponding bivector field will be denoted vr^cd. The reduced nonholonomic bracket describes 
the reduced dynamics in the sense that the nonholonomic vector field X„h is /c-related to the (almost) 
Hamiltonian vector field Xy^^ = {■,T-lTi}rcd associated to the reduced Hamiltonian Hn defined by the 
condition T-lj\^ = H-ji o p. 
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Remark 4.9. The reduction of nonholonomic systems performed by Bates and Sniatycki in [2] shows 
that it is possible to define a 2-form a;,.;,d on TZ which is non-degenerate along a distribution C C TTZ. 
The definition of C is given by C := Tp{U) where C/ := C R (C n V)^m and V is the distribution on 
M. tangent to the orbits of G. In fact, the pair (C, i^^ed) is just the pair associated to the almost Dirac 
structure given by the bivector field vTrod in the sense of Corollary 13.41 o 

The analysis of the reduction of the bivector fields vr^^ in the family ^ that are related to tt^^, by a 
dynamical gauge B follows from: 

Proposition 4.10. Let (P, vr) he an almost Poisson manifold and B a 2-form on P such that the 
endomorphism (Id — B^ o n^) : T*P — > T*P is invertible. If the Lie group G, acting freely and properly 
on P, preserves the almost Poisson structure vr and leaves B invariant, then G preserves the bivector 
field TT^ obtained by the gauge transformation of tt associated to B. 

Proof. In view of equation (|3.26|) . we see that (vr^)^ is a composition of invariant maps and we conclude 
that TT^ is invariant as well. □ 

As a direct consequence of this proposition we have: 

Proposition 4.11. IfG is a symmetry group of the nonholonomic system (-M,7rnh,^x) and B is aG- 
invariant dynamical gauge on A4, then {•, -j^ is G -invariant. In particular there is a reduced bivector 
field TT^^^B on the reduced space TZ that determines a well defined bracket {•, ■}^^^b on TZ satisfying 

{/,5LdB op(m) = {/op,5op}fXm), forf,g E C7~(7^). 

Moreover, the reduced bracket {■,-}rcdB also describes the reduced dynamics in the sense that X-^^ = 

There is a good reason why we did not denote the reduced bivector field tt^^^b by vr^^^ and that is 
that in general the reduced bivector fields 7r,.;,d and 'k^^^b need not be gauged related. We shall see this 
explicitly with the analysis of our mechanical examples (Remark 15. 5|) . 

Remark 4.12. Since the almost Dirac structure L^^ given in (j4.37p is the graph of a bivector (see 
Proposition 14.31 (i)), then the reduction of L^^ as an almost Dirac structure is simply the reduction 
of the bivector TTnh in the classical way. The same observation is valid for the reduction of the almost 
Dirac structure tb[Lt^^^) defined in ()4.40p since tb{Lt^^-^) is also the graph of a bivector {{Q-m ~ is 
non-degenerate on C). o 

The G-Chaplygin case. 

If the reduced bivector fields vr^cd and tt^^^b happen to be everywhere non-degenerate, then they are 
gauge-related. This is the scenario that one finds after reduction of external symmetries of G-Chaplygin 
systems (see [TS1E7]). These systems are characterized by the property that the tangent space to the 
orbits of the symmetry group exactly complements the constraint distribution on the tangent space 
TQ of the configuration manifold Q. 
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In this case there exist unique non-degenerate 2-forms ^1.,.^^ and ri^^^-s on TZ satisfying 



Id, 



Id, 



where as usual, ftl^^{X) = —ix^rcd and ft^^^^g{X) = —ix^^^dB for ah X G X(TZ). In particular, the 
reduced equations can be written as: 

which gives different almost symplectic formulations of the reduced equations (compare with the results 
in [22]). The bivector fields vr^cd and tt^^^b are gauge related via the 2- form iJ^od := ^^red — ^red^ ^ 
and we recover the Chaplygin Hamiltonization studied in [16]. Moreover, the 2-forms ri^ed and ^,.^^3 
satisfy 

(p*aed)|c = and (p*acdB)lc = i^M - B)c, 

and thus B^^^ verifies that {p*B^^^)\c = Bq. In other words, we have the following commutative diagram 



L 



""nh 
P 



(4.42) 



Trod ^ -^'^rodS • 



We stress that one does not have such a commutative diagram in more general situations, see Remark 
[531 



4.4 Hamiltonization 

We continue with the notation of the previous sections and suppose that G is a symmetry group for 
our nonholonomic system. The solutions of the reduced equations on the reduced space TZ = Ad/G are 
the integral curves of the reduced vector field Xy^^ and preserve the reduced Hamiltonian l-L-ji- 

The issue of Hamiltonization in our context concerns answering the question of whether the reduced 
vector field X-^^^ on TZ is Hamiltonian. Our candidates for the Hamiltonian structure come from the 
reduction of the (invariant) bivector fields that belong to the family 5- 

It turns out that the above condition is too restrictive. We relax it by asking that the vector field 
X-^^ can be rescaled by a basic^ positive function (^9 : 7^ ^ M in such a way that the resulting vector 
field (fX-^i^ is Hamiltonian. 

In view of Proposition 14. 1 ll for any G-invariant bivector field vr^^ belonging to 5^, the rescaled vector 
field ipX-j-i^ satisfies 

{if7^,^,Bf{dTin) = -vX^^. (4.43) 
Hence, we are interested in finding a bivector field vt^.^^b satisfying 

[(^^,,dS,V5vr,,ds] =0, (4.44) 

that is, we want to find vr^.^^s conformally Poisson (see Definition [2]) . 

*We use the term basic with respect to the fibered structure of TZ inherited from M. That is tp — (p o t where 
r : 7?. — > Q/G is the bundle projection and : Q/G — >■ M. 
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Definition 5. If there exist an invariant bivector field vr^^ belonging to ^ and a strictly positive, 
basic function (p : TZ ^ M such that (j4.43p and (j4.44p hold, we say that the nonholonomic system 
is Hamiltonizable. Moreover, we say that the reduced equations are Hamiltonian in the new time r 
defined by dr = ^dt (see discussion below). 

In Section [5] we will extend the discussion of Section [2] and show that all generalized rolling systems 
are Hamiltonizable. The table (jl.ip in Section [1] shows how different scenarios of the Hamiltonization 
scheme described above are realized according to the rank of the matrix A. 



Time reparametrizations 

It is common in the literature to interpret the rescaling of the vector field X-^^ by the basic positive 
function as a nonlinear time reparametrization. One argues as follows, let c{t) G 7^ be a flow line of 
X-^^ (i.e. ^{t) = X-}i^{c{t))). Introduce the new time r by integrating the relation 

1 

ar = — , , dt. 

V?(c(t)) 

Since ip > 0, the correspondence between t and r is one-to-one and one can express t as a function of 
T. The curve c(r) := c(t(T)) is checked to be a flow line of (pX-^^ (i.e. ^(t) = (/9(c(r))X-^^(c(r))). 
This interpretation of the rescaling is quite subtle. The definition of r depends on the particular flow 
line c(t), so different initial conditions induce different reparametrizations. It is therefore not possible 
to interpret the time rescaling as a "global" operation. This contrasts with the natural procedure of 
multiplying the vector field X-^^ by the positive function tp. 

Remark 4.13. One might wonder why we only care about basic and not arbitrary functions ip : 
TZ M^. A detailed answer to this question would involve a careful study of the structure of the 
equations of motion that would gear us away from the main subject of this paper. We refer the reader 
to |17j where one can find a very good discussion on the Hamiltonization of G-Chaplygin systems. We 
simply mention that physically, the fact that ip is basic means that it is independent of the momentum 
variables and only depends on the (reduced) configuration variables. Thus, the time reparametrization 
changes the speed at which the trajectories are traversed depending on the position of the system but 
independently of the velocity itself. It is shown in [T7] how in order to obtain Darboux coordinates 
for the reparametrized system, one can keep the same (reduced) configuration variables but should 
rescale the momenta by ^. Since ip> is basic, the rescaled momenta continue to depend linearly on the 
velocities. 



Measure preservation 

Hamiltonization is strongly related to the existence of invariant measures. Suppose for simplicity that 
we are dealing with a G-Chaplygin system. As mentioned before, in this case the bivector field tt^^^b 
is everywhere non-degenerate and the reduced equations can be written as: 
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where ^^.^^b is the non-degenerate 2-form on TZ induced by tt^^^b . It fohows that the scaled vector field 
ipX-j-i^ satisfies 

Hence, Hamiltonization in this setting amounts to finding a positive function 93 such that the 2-form 
^^rcd^ is closed (which under our hypothesis, is of course equivalent to (|4.44p ). Suppose for a moment 
that this is the case so (7^, ir2j.^jjs) is a symplectic manifold. It follows from Liouville's theorem that 

the vector field ipX-^^^ preserves the symplectic volume (^f^red^ ) where m = 2 dim 7^. Therefore, 
the volume form (^)'"~^(f^rcds)"^ is preserved by the vector field Xy^^. 

The above argument shows that a Hamiltonizable G-Chapligyn system possesses an invariant mea- 
sure. One might wonder if the reciprocal statement is true, namely, if any G-Chaplygin system with 
an invariant measure is Hamiltonizable. The celebrated Chaplygin's reducing multiplier Theorem [12j 
demonstrates that the answer is positive if m = 2. For m > 2, a characterization of the systems for 
which this is true is an open problem. Interesting examples where this holds for arbitrary values of 
m have been found by Fedorov and Jovanovic in the study of the multidimensional Veselova problem 
|17) . See also the discussion in [16] where a candidate for the conformal factor ^ is given under the 
hypothesis that there exists a preserved measure, and JQ] where a set of coupled first order partial 
differential equations for the multiplier ^ are given. 

In the case where the nonholonomic system is Hamiltonizable but the corresponding bivector field 
TT^^^B is degenerate at some points in TZ, one can repeat the above argument to conclude that the 
reduced system preserves a measure on every leaf of the symplectic foliation of TZ corresponding to the 
Poisson bivector field ipir^^^B . However, this does not imply the existence of a smooth invariant measure 
on TZ (and this is the motivation for Fernandez, Mestdag, and Bloch [19] to talk about Poissonization) . 
An example of this situation is given by the reduction of the Chaplygin sleigh (see the discussion in 
|2U| ) that exhibits asymptotic dynamics that contravene the existence of a global invariant measure. 
The problem of the existence of a global invariant measure in this case is most naturally attacked by 
considering the modular class of the Poisson manifold (TZ, ^tt^^^b ) , see [3S] . 

5 Back to the Examples: Hamiltonization and Integrability 

In this section, we analyze the generalized rolling systems presented in Section [2] using the geometric 
framework that was developed in the previous sections. In subsection 15.11 we provide the geometric 
interpretation for the brackets {•, -iRankj and {•, -jRankj presented in Section [2j In 15.21 we consider the 
Hamiltonization and integrability of generalized rolling systems in detail and finally, in l5.3l we explicitly 
show that the brackets {•, -jRanki and {•, ■y^^^i^2 twisted Poisson. To our knowledge, this is the first 
time that the appearance of such structures is made explicit in the field of nonholonomic mechanics. 

5.1 The geometry of the rigid bodies with generahzed rolling constraints 

We begin by computing the nonholonomic bracket for the motion of a rigid body subject to generalized 
rolling constraints as introduced in Section [2j 
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Nonholonomic bracket via the non-degenerate 2-section 

Consider again, as in Section [21 the motion of a rigid body in space subject to a generalized rolling 
constraint as in ()2.3p . That is, the constraint relates the linear and the angular velocities of the body 
X = rAu = rAgft, where the matrix A satisfies any of the conditions of Definition [1] and cj, is the 
angular velocity written in space and body coordinates, respectively. 

Recall that the configuration space for the system is Q = S0(3) x R^. Denote by A (respectively, p) 
the left (respectively, right) Maurer-Cartan form on S0(3). Upon the identification of the Lie algebra 
so(3) with by the hat map ()2.4p we think of A and p as valued 1-forms on S0(3). For a tangent 
vector Vg € Tg S0(3) we have 

i^ = p{9){vg), n = X{g){vg), 

where cj (respectively, CI) denotes the angular velocity vector written in space (respectively, body) 
coordinates as discussed in Section [2j 

The Maurer-Cartan forms A and p are related by X{g) = g^^p{g) and satisfy the well-known 
Maurer-Cartan equations 

dp=[p,p], dX = -[X,X], 

where [•, •] is the commutator in the Lie algebra. For the rest of the section we will use three dimensional 
vector algebra notation in our calculations with differential forms and vector fields. In our convention, 
the scalar product of differential forms should always be interpreted as a wedge product (and is thus 
anti-commutative!). The Maurer-Cartan equations take the form 

dp = X p, dX = (5.45) 

where "x" denotes the standard vector product in M'^. 

The constraint distribution D, defined by the generalized rolling constraints, can be expressed in 
the terminology of subsection 14.11 as the annihilator of the M^-valued 1-form e on Q given by 

e = dx — rAp = dx — rAgX. 

We consider the (global) moving co-frame {A, dx} for T*Q that defines fiber coordinates (M, p) 
in the following sense. A co-vector aq £ T*Q is written uniquely as = M • A + p • dx, for a 
certain (M, p) € M'^ x M^. The Legendre transform Leg : TQ — )• T*Q associated to the kinetic energy 
Lagrangian ()2.6p is defined by the rule: 

M = m, p = mx. 

Physically, p is the linear momentum of the body while M is the angular momentum of the body about 
the center of mass written in body coordinates. 

In order to deal with the constraints, it is more convenient to work with the global moving co- 
frame {A, e} for T*Q. We denote by (K,u) the fiber coordinates defined by this co-frame. Putting 
K-A + u- e = M- A + p-(ix implies 

K = M + rg'^A^p, u = p. 
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Along the constraint submanifold A4 = Leg(P) we have p = mrAgfl so 



which is the expression for the kinetic momentum obtained in ()2.9p . Notice that is a vector bundle 
over Q and that K is a natural coordinate for the fibers of A4. In what follows we will use the 
components of g,x, and K, as redundant coordinates on Ai. 

Denote by = {X[,X^,X^) the moving frame of S0(3) that is dual to A = (Ai,A2,A3). The 
components of are the left invariant vector fields on SO (3) obtained by left extension of the canonical 
basis of M^. Along the points of the constraint subbundle A4, the non-integrable distribution C defined 
in (|4.32|) is given by 



C = span I X^ + r^M^ — , — I . (5.46) 

The canonical 2-form Qq on T*Q is given by 

Qq = -d{M • A + p • dx) = -d(K • A + p • e) 
= A • dK — K. ■ dX — p ■ de + € ■ dp, 

where "•" denotes the usual scalar product in R^. 

To compute de we use the identity dg = gX where " denotes the hat map (j2.4p . Using the 
Maurer-Cartan equations (|5.45p we get 

de = —rA{dg)X — rAgdX = —rAg{X x A) — rAgdX = rAgdX. 

Therefore, 

= A • dK - K • dA - p • {rAg dX) + e ■ dp 
= A • dK - (K + rg'^A^p) -dX + e-dp. 

Let i : M. ^ T*Q denote the inclusion. Since p = mrAgO, along A4, we have 

Qm ■■= i^*{^q) = a • dK - (K + mr'^g'^A^Agn) • dA + t*(e • dp). 

Since e vanishes along the non-integrable distribution C, we get the following expression for the restric- 
tion Qc of i^*{^q) to C: 

fic = A • dK - (K + mr'^g'^A^Agn) ■ dX. 
With this expression for fic we are ready to show: 

Proposition 5.1. The nonholonomic bracket {•,-}nh on A4 for the generalized rolling system is given 
in the redundant coordinates {gij,xii.,Ki}, i,j,k,l = 1,2,3, for Ai by 

{xi,Ki},,^, = r{Ag)ii, {gij,Ki}^^ = -e]i gik, {Ki,Kj}^^ = -e-j(K + mr'^{g^ A'^ Agii))i , 

with all other combinations equal to zero. In the above formulas the Einstein convention of sum over 
repeated indices holds and e^j denotes the alternating tensor, that equals if two indices are equal, it 
equals 1 if{i,j,k) is a cyclic permutation of (1,2,3), and it is equal to —1 otherwise. The entries of 
{g,x) € Q are denoted by {gij,Xk) and those o/K by Ki. 
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Proof. We will rely on the identity ()4.36p that characterizes the associated bivector field TTn^ in terms 
of r^c • Contraction of 0,c by the elements in the basis of C given in (j5.46p gives 

i(x^+raTAT-S-)^c = dK - (K + mr^ g'^ Agfl) x A , i^Q^ = -A, 

where we have again made use of the Maurer-Cartan equations (j5.45|) . It follows that 

\x.^+r9TAT£^iK+mr^gTATAaSi)x^)^C = dK, i_,,^gJ_J^C = rAgX = (ix|c. 

Therefore, according to (|4.36p we get 



7rL(dK) = -X^ - rg^A^^ + (K + mr^g^A^AgQ) x A, 



(5.47) 



In addition, for any canonical vector € M'^ we have d{g ej) = {g ej) x A, so 

d 

7rL(%"^eO) = {g~'^ei) x — . 

The proof follows using the above formulas and recalling that for any /, 5 S C°°{A4), we have {/, g}nh = 
-dfi-KUdg)). 

□ 

Finally, we state without a formal proof that the nonholonomic vector field X^h on Ai is given by 

X„h = fi-XL + rA5n- — + (Kxn).— . (5.48) 

ax aK 

The above expression can be shown by taking into account the equations (|2.1ip . the constraint (|2.3p . 
and the definition of fi, or, alternatively, by computing the the almost Hamiltonian vector field X-^j^ = 
—■K\^{d'HM) corresponding to the Hamiltonian Hm that coincides with expression ()2.10p . The latter 
approach requires one to write in terms of K and g as was done in Section [2] for the different values 
of the rank of A. 

The gauge transformation of the nonholonomic bracket 

We will now construct a gauge transformation of the nonholonomic bracket in the sense of subsection 
14.21 We are interested in describing the same dynamics so we look for a 2-form B that defines a dy- 
namical gauge transformation as introduced in Definitional In our case, the distinguished Hamiltonian 
is T-Lm that has as its associated almost Hamiltonian vector field. 

Following |21j . we consider the bi-invariant volume form v on S0(3) oriented and scaled such that 
v{Xi ,X^) = 1. We consider the natural extension of z/ as a 3-form on Q = S0(3) x M^. Denote 
by P G n^{T*Q) the 3-form given by = v*!/ where v : T*Q — > Q is the canonical projection. We can 
write = gA ■ (A X A). 
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Let -B be the 2-form on A4 given by 

B = mr\ix^^i*i^), 

where, as before, l : Ai ^ T*Q is the inclusion. Note that i? is a semi-basic 2-form on A4 that 
vanishes upon contraction with the nonholonomic (almost) Hamiltonian vector field X^^. Therefore, 
by Proposition I4.6t we can perform a dynamical gauge transformation of the nonholonomic bivector 
field TTnh by the 2-form B to obtain another bivector field vr^^ that also describes the dynamics of our 
problem. 

Using the Maurer-Cartan equations ()5.45p and the expression ()5.48p for X^ih, we obtain 

B = -mr'^fl ■ dX. (5.49) 

To compute the bivector field vr^^ associated to the gauge transformation we use equation (j3.25p . For 
an arbitrary one-form a on we have 

(vrfj" (a + y«^(^)i?) =7rL(a). 

Setting a equal to A and dx and using ()5.47p and ()5.49p we obtain 

(vrfj«(A) = ^, (n^Jidx) = vAg^. 
Similarly, putting a = dK and noticing that 

V (dK)^ = -ixL^ = mr^n X A, 

nh^ ' 

we deduce 

(vrfj«(dK) = -X- - r/A^|- + (K + mr^g^A^Ag - E)^) x A, 

where E denotes the 3x3 identity matrix. 
The above formulas imply 

Proposition 5.2. The gauged nonholonomic bracket {•, ■}^^^ on M, associated to the bivector field vr^, 
is given in the redundant coordinates {gij,Xk,Ki), i,j,k,l = 1,2,3, for Ai by 

{xi,Ki}^^ = r{Ag)u, {9ij,Ki}^^ = -e% ftfc, {i^i, i^ilf^ = + mr'^^g^ A^ Ag - E)n)i , 

with all other combinations equal to zero. 

Reduction of the symmetries 

Recall that the Lie group introduced in section 12.41 acts on the configuration space Q and that its 
lift to TQ leaves both the Lagrangian and the constraints invariant. From the discussion in section 
14.31 (and the regularity of the action) it follows that the reduced space TZ := M./H \s equipped with a 
reduced bracket {■, - jrcd determined by condition ()4.4ip and that describes the reduced dynamics. 

We are now ready to give the geometric interpretation of the bracket {•, -jRankj introduced in section 
[2] for the different values of the rank of A. 
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Theorem 5.3. The reduced bracket {•, -jrcd on TZ is precisely the restriction of the bracket {-j-lRankj 
(defined in section\^ to the Casimir level set ||7|| = 1, for the different values j = 0, 1, 2, 3, of the rank 
of A. 

Proof. Recall from section [2^ that the reduced space TZ can be identified with xM^ with redundant 
coordinates (')',K). Therefore, it makes sense to compare the two brackets on the Casimir level set 
[|7|| = 1 of the space (7,K) € x R^. 

Moreover, from the expression of the projection p : A4 —j' TZ given by (|2.18|) . and condition (|4.4ip . 
it follows that the reduced bracket of the (redundant) coordinate functions (7, K) can be computed 
using the formulas obtained in Proposition 15.11 (notice that 7 = ((731, (732, 533))- 

The proof is completed by considering the particular form of A for the different values of its rank 
given in Definition [U and by writing the bracket {/i,/2}rcd of arbitrary functions /i,/2 € C°°{7Z) in 
terms of the derivatives ^ and ^ using Leibniz rule. □ 

We now turn to the study of the reduction of the gauged nonholonomic bracket {•, -j^. First of 
all notice that the 2-form B that defines the gauge transformation is written in (|5.49|) in terms of left 
invariant objects on S0(3). Since the symmetry group H acts by left multiplication on the S0(3) 
factor of Q, it follows that B is invariant under the cotangent lifted action. Therefore, in accordance 
with Proposition 14. IH the gauged bracket {•, -j^ drops to TZ where it defines the bracket {•, -y^^^B that 
determines the dynamics. As usual, the corresponding bivector field on TZ will be denoted by ti^^^b. 

In analogy with Theorem 15.31 we have 

Theorem 5.4. The reduced bracket {•, -^^^^b on TZ is precisely the restriction of the bracket {■, - j^^nkj 
(defined in section\^ to the Casimir level set ||7|| = 1, for the different values j = 0, 1, 2, 3, of the rank 
of A. 

The proof is identical to that of Theorem 15.31 except that one uses the formulas obtained in Propo- 
sition 15.21 

According to the above theorems and the discussion in Section [21 the the properties of the reduced 
brackets {•, -jrcd and {•, - ^^^^b are summarized in the table p.ip presented in the introduction (Section 
[1]) according to the different values of the rank of A. 

Remark 5.5. Notice that the reduced bivector fields -k^^^ and tt^^^b are not gauge related. Indeed, 
from table (jl.ip one sees that for any value of the rank of A, only one of the two bivector fields tt^^^ 
or VTj.^jjS has an integrable characteristic distribution. It follows from Theorem 13.111 that there cannot 
exist a gauge transformation between their graphs. This situation is opposite to the one described by 
G-Chaplygin systems in diagram ()4.42p . 

o 

Remark 5.6. Recall from subsection 12.51 that as the rank of A increases, the constraint distribution is 
less integrable or "more nonholonomic" . Table (jl.ip seems to suggest that it is appropriate to perform 
a gauge transformation by the 2-form B when the nonholonomic effects are more important, while the 
reduction of the standard nonholonomic bracket works better for weaker nonholonomic effects. 

o 
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5.2 Hamiltonization and integrability of rigid bodies with generalized rolling con- 
straints 

According to the notion of Hamiltonization introduced in Section SHI (Definition ED, and the table (jl.ip . 
it immediately follows that the problem of the motion of a rigid body subject to a generalized rolling 
constraint is Hamiltonizable for any value of the rank of A. 

If the rank of A equals (respectively, 3) the reduced equations are Hamiltonian with respect to 
the bracket {•, ■}^^^ (respectively, {•, -j^^js). Recall that in both cases the reduced dynamics correspond 
to classical rigid body motion (with modified inertia tensor I + mr'^E if the rank of A equals 3). 

If the rank of A equals 1 or 2, the analysis of the Hamiltonization is a bit more delicate but it 
also follows directly from Definition [5] and the table (jl.ip . In the case rank ^ = 2, it follows that the 
reduced equations are Hamiltonian in the new time T2 defined by dT2 = -;^dt and with respect to the 
bracket (/52{-,-}redS where 

ipii'j) = Vl - mr2 (7 • (I + mr^E)-^'y). (5.50) 

Note that (p2 is a basic function on TZ corresponding to the restriction of (j2.15p to the level set ||7|| = 1. 

Analogously, if the rank of A equals 1, the reduced equations are Hamiltonian in the new time ti 
defined by dri = -^dt and with respect to the bracket -Ircd where 

V5i(7) = ^Jl + mr'^ il-l-^l)- (5.51) 
Integrability of the reduced equations 

In view of the Hamiltonization of the problem, the integrability of the reduced equations of motion 
(|2.1ip can be easily established using the celebrated Arnold-Liouville Theorem for classical Hamiltonian 
systems, see e.g. [1]. 

Indeed, for any value of the rank of A, the reduced equations are Hamiltonian on IZ (after a time 
reparametrization if rank A = 1,2). Independently of the rank of A, the symplectic leaves Oa of the 
foliation of IZ correspond to the level sets Ci(K, 7) = K • 7 = a and can be shown to be diffeomorphic 
to the tangent bundle T of the sphere (see the discussion in chapter 14 of [30] for the coadjoint orbits 
on 5e(3)*). 

Once the value of a is fixed, the reduced equations (j2.1ip can be seen as a two degree of freedom 
classical Hamiltonian system on Oa (again, after a time reparametrization if rank A = 1,2). These 
equations possess two independent integrals, the Hamiltonian T-Lu-, and F = K • K, whose joint level 
sets are compact in Oa- It follows from the Arnold-Liouville Theorem that these level sets are invariant 
two-tori and the dynamics are quasi-periodic on them (notice that the flow on the tori is rectilinear 
but not uniform if the rank of A is 1 or 2). 

The Arnold-Liouville Theorem also tells us that the reduced equations are integrable by quadratures 
(after the time reparametrization if the rank of ^ is 1 or 2). 

Finally, we state without proof that the reduced equations of motion ()2.1ip preserve the measure 
/x(7) o" A dKi A dK2 A dK^ where a is the area form of the sphere S^, and the basic density /i : — > M 
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is given by 

1 if rank j4 = 0, 3, 



1 

, ¥'2(7) 

where (pi,tp2 £ C°°(S^) are defined in ()5.5ip and ()5.50p respectively. 



4-^ if rank A = 1, 
if rank A = 2, 



5.3 Twisted Poisson structures for rigid bodies with generalized rolling constraints 

In Section I3.2[ we presented twisted Poisson structures which have been extensively studied in other 
contexts but not in mechanics. Now, we will show explicitly that twisted Poisson structures appear 
naturally in the study of nonholonomic systems. 

Rigid body with generalized rolling constraints of rank 2 

Here we show that the bracket {•, •}^^j,„ij2) in addition to being conformally Poisson, is twisted Poisson. 
Note that this cannot be the case for the other bracket {•, •}Rank2 that describes the dynamics since, as 
shown in Section [21 its characteristic distribution is not integrable. 

Recall from the discussion in 12.41 that {•, •}Rank2 should be considered as a bracket on the reduced 
space TZ = S^ xR^ with redundant coordinates (7,K). The characteristic distribution of the bracket 
is integrable and the leaves Oa of the foliation are the level sets Ci('y,K) = 7 ■ K = a. By regularity 
and integrability of the characteristic distribution, it follows from Corollary 13.71 that the bracket is 
(/>-twisted. The value of the 3-form (j) is given in the following. 



Theorem 5.7. The bracket {■,-}'^^n]^2 defined in { 2. 14\ ) (that in particular describes the reduced dy 



namics of the Chaplygin sphere for the appropriate choice of A), is a (j)-twisted Poisson bracket with 
(j) = —dB where 

= mr2(n-7)a, (5.52) 

and where a denotes the area form of the sphere ||7|| = 1. 

Proof. The idea of this proof is to show that the bracket is gauge related to a Poisson bracket via 
the 2-form —B. Thus, by Proposition 13.121 the bracket is (— (i;B)-twisted Poisson. More precisely, we 
will show that {-j-jH^nka is (— S)-gauge related with the bracket {•,-}Ranko defined in ()2.17p and that 
coincides with the Lie- Poisson bracket on se(3)*. 

According to Theorem [531 we denote the bivector field associated to the bracket {•, - Y^^^^^ by T^^cdB- 
The bivector field 'k^^^b and the 2-form B verify hypothesis of Proposition 14.71 and thus the gauge 
transformation of tt^^^s associated to B is again a bivector field that we will denote it by tt^^^ . Relying 
on equation ()3.25p and writing in the redundant coordinates (7,K) as 

B = ^?Tir^(ri • 7) 7 • (c?7 X d'y) 

one gets that 

(dK) = K X A + 7 X A, (vrl.)«(d7) = 7 x A, 
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and the proof is complete. 

□ 

The conformal factor and the 3-form (p. In accordance with Proposition 13. 10] since the bracket 
{•, •}^j,nk2 is both conformaUy Poisson and twisted Poisson, there is relationship between the conformal 
factor ip2 (given by (|5.50p ). and the twisting 3-form cp (defined in Theorem 15 .jp . 

We leave it to the reader to check that on the leaves Oa of the foliation of TZ corresponding to 
the bracket {•,'}Rank2i ^be 3-form cj) coincides with Tp := -^dip2 A where 2-form Q is given in the 
redundant coordinates (7, K) by 

= ^ (K - mr2(n • 7)7) • {d-y x d-f) - {dK x ^7) . 

This choice of Q satisfies the conditions of Corollary 13.21 for the graph of the bivector field vr^.^^s 
corresponding to {•, •y^^,^i^2 ® T*Tl. 



Rigid body with generalized rolling constraints of rank 1 

A completely analogous analysis can be performed if the rank of the matrix A equals one. This time 
it is the bracket {•, - jRanki that is both twisted and conformally Poisson. In analogy with Theorem 15.71 
we have 

Theorem 5.8. The bracket {•,-}R,anki defined in i2.16\) . is a (f)-twisted Poisson bracket with (j) = dB 
with B given by expression ()5.52p . 

The proof is the same to the Rank 2 case. The bracket {•, -jRanki is S-gauge related with the bracket 
{•1 •}Ranko defined in ()2.17p and that coincides with the Lie-Poisson bracket on se(3)*. 
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